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NOUVELLES MÉTHODE&O^ 

POUR LA DÉTJERMINATIO n( 

DES ORBITES DES COMÈTESe-^ 


XjE problème dont je m’occnperni dans ce Mémoire, consiste à 
déterminer l’orbite d’une comète d’après trois observations don- 
nées de sa longitude et de sa latitude. Depuis Newton , qui le 
premier adonné pour sa solution des constructions géométriques 
fort ingénieuses , ce problème a fait successivement l’objet des 
recherches d’un grand nombre de géomètres. 

Parmi ceux qui l’ont traité avec le plus de succès, on doit 
distinguer particulièrement Lambert, qui a donné de très-beaux 
théorèmes sur le mouvement des planètes et des comètes dans 
son ouvrage intitulé : Insigniores orbilœ cometarum proprictates, 
et dans les Mémoires de l’Académie de Berlin , année 1771. 

Dans ce môme recueil , année 1778 , La Grange a discuté les 
principales méthodes connues jusqu’alors , et après avoir fait 
connoltre les causes de leur imperfection , il a donné l’analyse 
complète du problème, fondée sur une belle théorie à laquelle 
il a ajouté ensuite de nouveaux développemens dans le volume 
de 1783. Cette méthode n’auroit sans doute rien laissé à desirer, 
si son illustre auteur en eût fait l’application à des exemples , 
ce qui l’auroit conduit lui-même à y apporter les modifications 
nécessaires pour en rendre l’usage facile dans la pratique. 

La Place , dans les Mémoires de l’Académie des Sciences de 
Paris , année 1780, a proposé une autre méthode' qui revient à 
supposer connue une portion infiniment petite de la trajectoire 
apparente de la comète , et à conclure de cette portion la gran- 
deur et la position de la trajectoire vraie. Par la supposition 
des quantités infiniment petites , les formules se simplifient et 
conduisent à une solution peu compliquée; mais la difficulté est 
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de déterminer avec précision les coefficiens différentiels du pre- 
mier et du second ordre , tant de la longitude que de la latitude. 

Il semble , au premier coup-d’œil , qu’en liant ensemble plu- 
sieurs observations par la méthode des interpolations , on en 
déduira les coefficiens dont il s’agit avec d’autant plus d’exac- 
titude qu’il y a plus d’observations combinées ; et cela auroit 
lieu en effet, si les observations étoient exemptes d erreur, ou 
si les lieux de la comète étoient calculés d’après une formule 
exacte j mais il en est autrement dans l’état réel des choses. 
Comme toutes les observations sont affectées de quelqu erreur , 
et que cette erreur ne suit aucune loi d’une observation à 1 autre, 
jl s’ensuit que plus on eombinerad’observations, et plus 1 erreur 
des coefficiens différentiels qui en sont déduits pourra devenir 


sensible. _ 

En effet, considérons trois observations de longitude a',a’,a , 
faites à des intervalles de temps que pour plus de simplicité 
nous supposerons égaux à l’unité. Ces longitudes répondront 
aux temps — 1 , 0 , +i, et la longitude pour un temps quel- 
conque t compris entre — i et + i , aura pour expression 
générale } 

* = o* + ï(«"'— o')t+i(a" — ia +a)f ; 
d’où l’on déduit les coefficiens différentiels 


dt * 


(<*) 


Considérons ensuite les cinq longitudes a, a!, a ", a!”, a”, qui 

répondent pareillement aux temps — a, — i, o, +i, +a, res- 
pectivement ; on en déduira la longitude au bout du temps t, 
x = a"+At + B <* + C* s + &c. , où l’on aura 


A ou -JT *= î («" — Ô — 75 °) 

at 

aB ou = j (®" — aa*+«0 — n (o"—- 3 a"+ o). 

Cil 


(*) 


Digitized by Google 


( V ) 

Supposons maintenant que l’erreur sur la longitude a' soit t'a'} 
alors , suivant les équations ( a ) , les erreurs des coëfficicns diffé- 
rentiels , dues à cette cause , seront 

„ ddx 




de 


— ta'. 


Mais par les équations (b) les erreurs de ces coëfficicns scroient 

ddx . 

-= T ta ; 


dx , 

t —— — — \tâ 
dt ’ 


de 


d’où l’on voit que dans le cas de cinq observations combinées , 
l’erreur des coëfficiens différentiels , due à la cause mentionnée , 
est plus grande dans le rapport de 4 à 3 , que celle qui a lieu 
dans le cas de trois observations. Elle augmenteroit encore si 
on combinoit ensemble plus de cinq observations. 

D’après ces réflexions, j’ai pensé que ce qu’il y avoit de mieux 
à faire dans le problème des comètes , étoit de partir des données 
immédiates de l’observation , et d’employer tous les moyens 
pour simplifier autant qu’il est possible, les formules et les 
équations qui servent à déterminer les élémens de l’orbite. C’est 
l’objet que je me suis proposé dans ce Mémoire. 

Je l’ai divisé en deux parties : la première comprend l’analyse 
générale du problème, avec deux applications détaillées de la 
solution qui en résulte aux comètes de 1781 et de 1769. 

Dans l’analyse je suppose , comme cela est indispensable, que 
les trois observations données ne comprennent pas un inter- 
valle de temps de plus de i 5 à ao jours , afin que les séries qui 
expriment les coordonnées , tant de l’orbite de la comète que de 
celle de la terre , soient suffisamment convergentes , et qu’on ne 
soit pas obligé d’employer les termes qui contiennent des puis- 
sances du temps supérieures à la troisième. 

Après avoir mis les équations générales du problème sous la 
forme la plus simple dont elles paroissent susceptibles , je déve- 
loppe en particulier, avec beaucoup d’étendue, le cas où les 
trois observations sont faites à des intervalles de temps égaux. 
Celte supposition, qui ne restreint guère la généralité du pro- 
blème, simplifie beaucoup les formules, et contribue même à 



( vj ) 

les rendre plus exactes , par la disparition de plusieurs termes 
dont il faudrait tenir compte , si les intervalles de temps entre 
les observations n’étoient pas égaux. 

Lorsqu’on ne fait aucune supposition sur la nature de l’orbite, 
le problème offre précisément autant d’équations que d’incon- 
nues ; mais alors il y a une circonstance dans laquelle les for- 
mules pourraient ne pas donner des résultats assez exacts. C’est 
lorsque l’orbite apparente est située à très peu près dans le plan 
d’un même grand cercle , et dans ce cas , la distance de la 
comète au soleil, lors do l’observation moyenne, diffère tou- 
jours très- peu de la distance de la terre au soleil. Si les méthodes 
analytiques ont peu de succès dans ce cas particulier , il en 
résulte au moins une connoissance utile sur la distance de la 
comète au soleil , laquelle servira toujours à diriger les promiers 
essais des calculateurs. D’ailleurs , en choisissant trois autres 
observations , à quelque distance des premières , on évitera 
l’inconvénient de cc cas particulier , à moins qu’il ne se ren- 
contre tout à-la-fois que la comète soit très-voisine du péri- 
hélie , et que la distance périhélie diffère très-peu de la distance 
de la terre au soleil ; circonstances qui, tout en faisant excep- 
tion , avanceraient beaucoup vers la connoissance de la véri- 
table orbite. • 

Mais comme on suppose communément , d’après les résultats 
de l’observation , que l’orbite est parabolique , celte condition 
donne une équation de plus que d’inconnues, et on a la faculté 
de choisir entre les diverses combinaisons des équations , celle 
qui doit conduire aux résultats les plus exacts. Je suis entré à 
ce sujet dans une discussion fort étendue. J’ai fait voir quelles 
sont les équations qui mèneraient, dans certains cas, à des 
résultats défectueux , et quelles sont celles sur lesquelles on peut 
établir la solution la moins sujette à être affectée des erreurs 
des observations. 

Les équations dont il s’agit donnent immédiatement par leur 
résolution les distances de la comète au soleil et à la terre ; il 
faut ensuite en conclure les élémens de l’orbite. Je donne pour 
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cet effet toutes les formules nécessaires, cl j'ajoute les moyens 
de recounoitre avec certitude si le mouvement est direct ou 
rétrograde ; si la comète marche vers le périhélie , ou si elle a 
déjà passé par ce point; si le nœud dont on a calculé la longi- 
tude est le nœud descendant ou le nœud ascendant ; de sorte 
qu’avec toutes ces directions les calculs pourront s’exécuter en 
quelque sorte mécaniquement, sans qu’on ait lieu de craindre 
de se tromper sur le signe d’aucune quuntité , ou sur la position 
d’aucun des points qu’il importe de déterminer. 

Après avoir traité complètement le cas où les intervalles 
de temps entre les observations sont égaux, il falloit au moins 
tenter de vaincre les difficultés d’analyse que présente le pro- 
blème considéré dans toute sa généralité. Je suis donc revenu 
sur les équations générales du problème , et après leur avoir- 
fait subir différentes réductions successives , j’ai trouvé assez 
heureusement qu’elles étoient susceptibles d’une solution pres- 
que aussi simple que celle du premier cas ; de sorte que sans 
aucune interpolation , on pourra appliquer immédiatement le 
calcul à trois observations données , quels que soient les inter- 
valles de temps qui les séparent, 

Venant ensuite aux applications qui sont relatives a la seconde 
comète de 1781 et à la comète de 176g , j’entre dans les détails 
qui peuvent jeter du jour sur toutes les applications en général. 
J’examine successivement les divers systèmes d’équations qu’on 
peut former , j’en donne la résolution ; je dévelqppe les calculs 
principaux avec beaucoup d’étendue, et en général avec plus 
de précision qu’il n’est nécessaire dans les applications ordi- 
naires : mais celles- çi étant données pour exemple, et devant 
faire juger du degré d’exactitude de la méthode, on ne pouvoir 
y mettre trop de soin. Les résultats, au surplus , eu sont très- 
satisfaisais , et on verra que par cette méthode il est possible de 
déduire de trois observations données d’une comète, dans un 
assez petit espace de temps, des valeurs fort approchées des 
vrais élémens de son orbite. 

Cependant la connoissuncc de lu véritable orbite ne scia 
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acquise d’une manière certaine, que lorsqu’on aura pu satisfaire 
à trois ou à un plus grand nombre d'observations , dont les 
époques soient assez éloignées entre elles. 

La seconde partie contient les méthodes nécessaires pour 
corriger , à l’aide de ces observations , les élémens connus par 
une première approximation. 

Je prends pour exemple trois observations de la comète 
de 1769 , faites à de grands intervalles de temps , et je fais voir 
comment , par l’emploi des corrections indéterminées ( que 
j’avois déjà indiqué dans les Mém. de l’Acad: des Sc. ann. 1787), 
on peut parvenir à déterminer la parabole qui satisfait le plus 
exactement possible aux observations. La méthode des correc- 
tions indéterminées a l’avantage de donner, par un seul calcul 
et sans aucun tâtonnement, le résultat de toutes les hypothèses 
peu différentes entre elles qu’on pourroit former sur la valeur 
des élémens. Il faut ensuite, lorsque toutes les conditions du 
problème sont exprimées convenablement , déterminer les 
coëfîiciens de manière à rendre les erreurs les plus petites qu’il 
est possible. 

Pour cet effet, la méthode qui me paroît la plus simple et 
la plus générale , consiste à rendre minimum la somme des 
quarrés des erreurs. On obtient ainsi autant d’équations qu’il y 
a de coëfficiens inconnus ; ce qui achève de déterminer tous les 
élémens de l’orbite. 

Comme la méthode dont je viens de parler , et que j’appelle 
Méthode des moindres quarrés, peut être d’une grande utilité 
dans toutes les questions de physique et d’astronomie où il 
s’agit de tirer de l’observation les résultats les plus exacts qu’elle 
peut offrir ; j’ai ajouté , dans une appendice , des détails parti- 
culiers sur cette méthode , et j’en ai donné l’application à la 
mesure de la méridienne de France , ce qui pourra servir de 
complément à ce que j’ai déjà publié sur cette matière. 
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PREMIÈRE PARTIE. 

Analyse du problème. 

I. Soit SAT (fig. 1) le plan de l’écliptique, S le soleil, AT 
l’orbite de la terre , S A le rayon dirigé vers le premier point 
d 'arias supposé fixe ; soient T et C les lieux de la terre et de 
la comète , en un même temps t. Du point C abaissez CK per- 
pendiculaire sur le plan de l’écliptique j menez KF, TE, 
perpendiculaires, et TG parallèle à SA; enfin joignez SC, TC. 

Nous appellerons*,^, *, les trois coordonnées SF, FK, KC 
de la comète, v son rayon vecteur SC , f sa distance à la terre 
TC; . * * ’ 

X, Y, les deux coordonnées SE, ET de la terre, V son 
rayon vecteur S T ; . 

a et Cia longitude et la latitude géocentrique de la comète , 
c’est-à-dire, les angles GTK, CTK. 

II. Si l’on prend pour unité la somme des masses du soleil 
et de la terre, et qu’on néglige par rapport à cette somme, la 
différence , d’ailleurs inconnue , entre la masse de la comète 
et celle de la terre; les équations différentielles du mouvement 
de la comète seront: 

ddx x ddy y ddz * 

~dF V ’ ~dF ~ ~ ’ df ~ S* 

et celles du mouvement de la terre , («) 

ddX _ X ddY _ Y 

di‘ ~ V* ’ dt * ~ V 5 ’ 

* n 

Ces équalious supposent en outre que le temps t est mesuré par 
le moyen mouvement du soleil réduit en parties du rayon , 


( 3 ) 

et qu’on fait abstraction de l’action réciproque entre la terre 
et la comète. 


♦ 


III. Cela posé , soient encore m , n , p , r, les valeurs respec- 
tives de x , y , z , r ; et M, N, R, celles de X, Y, V, lorsque 
t = o. Si on suppose qu’à compter de l’époque où t == o , les 
intervalles de temps ne soient pas trop considérables en deçà 
et au-delà de cette époque ; les coordonnées , tant de la comète 
que de la terre , seront déterminées avec une exactitude suffi- 
sante par les équations suivantes, où i’on a omis seulement 
les termes affectés de et des puissances supérieures de t. 

x ~ m + 

y — n + rit + A V + »''V 

z = p + p't + /> V + p "V • 

m 

X = M + M'< + M'/’ + M"V 
Y = N + N'/ y NV + N"V 

' . • 

Ce s valeurs doivent être substituées dans les équations (t),* 
afin da réduire au plus petit nombre les coëificiem indé- 
terminés , et comme on a 


, , , , f «i' + s’ + j»' 

v ~ x ^ * ( + af (mm’ + nn'+ pp') + &c. 


si on fait pour abréger 

k = mm + /ira' -f pp' , 

, K = MM'+NN', 


on aura 


v' -=. r* + ait + &c. , 
d’où résulte 




&c . , 


V* = R‘+ 3K/+&C.; 

JL - _L _ 

V’ “ R J R* + 


à 
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Jît la substitution faite dans l'équation différentielle en * 
donnera 

„ m mk m . * 

*1 TTt 5— ' " • « ^ ““ * ~ n »’ • I . 

a r 3 ’ If 6 r 3 t 

On obtiendra un semblable résultat des autres équations diffé- 
rentielles,, ce qui donnera , 

;* = m (, ~~ + 77 ) + (* ”67") 

' = n (v~7^ + J^) + *' G-S7-) 

/ ( 1 “^ + T7) +y G~6 t) *; w 

x=m G~^ + ^) + m ' G“ et 3 ) 

— G-^S)+-G-é) 

• , i . 1 . • „ t : ' . . • i\ » * . . . • 

IV. Au moyen de ces coordonnées , on peut déterminer la 

longitude * de la comète et sa latitude C , par les équations : 

. ' '4— Y tanglC x t 

tang a — rr , — .«= ÿ' 

0 x — X. cvs <e x — X 

De sorte qu’en général, pour un temps quelconque t qui n’ex- 
càik-ra pas certaines limites avant ou après l’époque choisie, 
on aura les deu^ équations : - j_ 

•‘(-■à+SM ■/ -(-£+£) 


- n ('-æ + æ)L 

-M-iff) ) 


= tang * 


( t' kt?\ 
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(4) 


( 3 ) 



:.r Y 

m 1 

-M* kt* \ 

1 -5- + . 1 




» \ 

ar a r / 

' ,Yi-— +— 



— M ( 

JL + üfl\ 

^ \ ar‘ + ar 1 / j 


tong-C 

aR‘ jü/ 

. t’ \ | 


cos « 


r’ \ 




»- «y 

1 6 r ) 




,/ 

t \ 




— M'{ 
1 ' 

t J 




6R/ 

V. Soit pour abréger : 


' 

* 



m' = 

M'+ P 

1 

n = N + r, 

fl' as 

N' + s' 

; (*) 

et de pins , 



“ 

1 

I 


a 

P R 5 



• * , 

h 

K 


(«) 

■ 

7 “ 

R‘ = f * 


kf 


Si on multiplie les équations (3) par i + — y — — - 5 ,et qu’on 

2r îr 

néglige toujours les t \ elles deviendront : 


t •+ » 


"( , + 37)1 


-;N 


= («ng » 


-îM (.f— fO 

-jM'.t 1 


( 6 ) 


F + f \ + 3rV cos* j 

On peut simplifier encore ces équations en omettant le facteur 

1 + -jp qui multiplie #»', //; car la snite du calcul prouvera 

que l’omission peut être faite sans qu’â en résulte d’erreur dans 
les termes qu’on doit conserver. 
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VI. Supposons maintenant que les trois observations don- 
nées de la comète répondent successivement aux temps : 

— 9, o, +9', 

en sorte que l’époque t — o réponde à l’observation moyenne. 
Soient 

les longitudes observées a" , a , a' , 

’ . et les latitudes b° , b , b'. 

Faisons de plus, pour abréger, 

tanga 0 — f°, tang a = f, tang a' = f ( 7 ) 

tan g b° tang b _ tang b' , 

cos a° ° ’ cos a cos a' ^ ’ 


Si dans les équations ( 6 ) on fait successivement t = o, t = —t , 
f = 2 i V t on aura les six équations suivantes': 

* • = /t* 


» — — ;N.f 1 9 /~tM«9* 

— -;Ni'9 J +-;N'-9’J ' [ 

r + IV — { N - 9'* 1 ,, U + IV— i M - » r ‘ 

+ iNC«' 3 — jN'-n 1 l + iM£4' 3 — JM-r* 

( 8 ) 

P “ W 


P - «/ = 


— iMf9’-t-iM'-9 J 


n . , V _ ■ ff* + — î M * 9'* 

P+ P~g[ 


VII. Après la substitution qui se fait immédiatement des 
valeurs de t ctdep, ces six équations se réduiront à quatre 
contenant les quatre inconnues g , g, »',//, sous forme linéaire. 
Il faut d’abord s’occuper de la résolution de ces quatre équa- 
tions. i 

Dans le cas présent, où les équations ( 8 ) ne sont exactes 



( fi ) 

qu’aux quantités près do l'ordre fi 4 , il faudra omettre, dans le 
calcul de l’élimination , tous les termes qui seroient de cet 
ordre ou d’un ordre supérieur, attention qui contribuera à 
simplifier les résultats. On pourroit aussi établir les formes 
des diverses inconnues n, p, p, en les représentant chacune 
par A + B8 + ÇS'-p 1) 9‘ + £ 8 fi' + F9'*+ &c.;on reconnoîtroit 
bientôt par la substitution , que rcs quatre inconnues ne con- 
tiennent aucun termo constant indépendant de fi et 9'; que 
la quantité p en particulier, ne contient pas même les termes du 
premier degré , et qu’dle doit être divisible par fifi', de sorte 
qu’on doit faire /a=A88' + B9*9'-t- CfiS'v On trouverapareille- 
ment que les trois autres inconnues doivent être de la forme 
A 8 + 139' + C8*+ D9 9'+ Efi'*, sans aller au-delà, afin de ne 
pas introduire dans les équations ( 8 ) des termes du quatrième 
ordre; et cette forme justifie l'omission que nous avons faite 

t‘ , • 

du facteur i + , dans les équations (6). On pourroit donc 

par les cociticiena indéterminés, effectuer la résolution des 
équations (8); mais on y parviendra aussi facilement par les 
moyens ordinaires, et voici le résultat de ce calcul présenté 
sous la forme la plus simple dont il paroit susceptible. 

VIII. Soit pour abréger, 

a = (f'g —fé) + (/v— yv ) + c f^-pg). (9) 

Cette quantité qui se présente d’abord dans le résultat de l’éli- 
mination mérite une attention particulière : elle peut se mettre 
sons la forme 

or-/) ig—g 0 ) -(/'— /°) {g -g), 

et alors on voit aisément qu’elle est de l’ordre 9fi' (9 + 5'), 
c’est-à-dire du troisième ordre, et cette extrême petitesse 
empêche le plus souvent qu’elle ne soit déterminée assez exac- 
tement par les données de l’observation. Sort de plus , 
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M (fg—fg) + N (g'— g‘) = B" 

M (f°g-fg°) + N (ff-g) = B 
M <Jg a -f°g) + N (g -tf) = B' 

M' ifg-fg) + *rc ï-.b) = fi* 

M'c f°g'-f'g°) + N'(g° — tf') = E 

+ N '(*-*»} = E*; 

1 

quantités entre lesquelles on a les relations 

. B° + B 4- B' = Ma • E^E + E^M'a ( 10 ) 
/°B° +/B +/'B' = Na /*E° + /E + /'E' = N'A 
g°B° + gB + g'B'= o • # D E’ > + #E+#'E'= o; 

les valeurs réduites de nos quatre inconnues seront : 


y , r^ ir B /< r 


») 


M = 


3 A 
3 A 

M 


6a 


88'E 


0*1 


8 ) 

(SB° — J' B') + ( 8E°— -8'E') 

DA 


+ + (iBT— |Mf) («•-»«' + l'p) 

a 

^ ( 8 /°B°— 8/'B') + ^ (8f E" — if E') 


3 A 


N 


+ ^«(8' — ?) + (VN'»— ;NO (8* — 8 8'+ 8' 1 ) 
p' =* ( •«* B °- V B ' ) + 8'g' E' ). 


a a 


IX. Avec ces valeurs, on calculera celles des six coefüciea* 
m , n , p, m, ri, p'i ce qui se fera par les équations 
m = f* + M n»' = f«' + M' 

»=/f.+ N n' — / •+■ N' (is) 

P — 8* P'—P' 

et ces six valeurs ne contiendront que les deux inconnues 
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( et a sous forme linéaire. Calculant ensuite la valeur de k par 
la formule 

k = mm + nri + pp , (i3) 

et substituant cette valeur ainsi que celle de » = — — — j, 

r 

dans l’équation 



on aura une première équation entre les quantités f et r, dans 
laquelle f ne montera qu’au sccdnd degré. 

Enfin on a de plus les deux équations 

m' + n' + p' = r' (*5). 

n'*+ p'*= * — ^ , (16) 

dans lesquelles il faudra faire de semblables substitutions pour 
obtenir les équations finales auxquelles se réduit la solution 
du problème. 

L’équation (16) appartient à une orbite elliptique dont le 
demi-grand axe = a ; mais comme on suppose ordinairement 

que l’orbite des comètes est parabolique , on pourra faire - = o , 

ce qui donnera une équation de plus que d’inconnues. Et cette 
circonstance pourra offrir diverses combinaisons dont on pro- 
fitera pour rendre moins difficile la détermination des incon- 
nues r et Ç. 

Nous donnerons ci-après le moyen de vaincre ces difficultés 
et de parvenir à des résultats dont la pratique puisse s’accom- 
moder. Mais avant tout, nous examinerons , avec tout le détail 
nécessaire, le cas où l’on a 9'= 3; c’est-à-dire où les trois 
observations sont faites à des intervalles de temps égaux. Cette 
condition apporte dans les calculs une grande simplification , 
et on peut toujours y satisfaire , en interpolant les observations 
faites à peu de distance les uues des autres. • 


■ 

v • 
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X. Soit donc fl' = fl , et les formules (ti) se réduiront aux 
suivantes : 

a# 3* 

* = TT® v . ’ 

a A * -• 

#•' = 7“ (B® — B') + (jM'a — 

2 A • 


» r = -(rB°-/'B'>+ (JIT-- 

p'= fi'). 

* • a ® . I' *, à * . * . -• 


Ôr j’observe qu’on peut omettre les termes affectés de fl* dans 
les valeurs de y et En effet, dans les équations (6) on a 

/fl* 

négligé , par rapport à y, la partie y-j , parce que le premier 

terme de / étant de la forme A 8 , cette partie serait devenue 

; et multipliée par fl , comme l’est / dans les équations (6) , 

elle serait montée au quatrième degré; ainsi elle a dû être exclue 
des équations (6) qui ne contiennent pas les termes affectés 
de t 4 ou 8 4 . D’un autre côté cependant , en considérant fl comme 
très-petit, Je coefficient A qui affecte fl doit être très-grand, 
afin que A fl ait la valeur finie qui convient à la quantité y; 
de Sorte que A9 3 désigne réellement une quantité de l’ordre 9*. 
Ainsi en conservant dans les formules précédentes les termes 
affectés de fl*, qui ne penvent être beaucoup augmentés par 
leurs coëfficiens , on n’ajoute pas à l’exactitude de ces formules 
dans lesquelles des quantités du même ordre ont été omises. 
Donc il faut supprimer entièrement ces termes, et on aura les 
valeurs très-simples . 1 ! 

«fl* _ 

y = — rn.fi » 1 

sa 

"'= C*7 ) 


i 


( io ) 


- (/ # B° — /'B') ' 

2 A 

/>’= ^ — jr'B' ) 

dans lesquelles l’ioconnue f ne se retrouve plus. 

XI. Avant d’aller plus loin, il convient de déterminer le# 
quantités M, N , M', N', qui dépendent du mouvement de la 
terre. Et d’abord , si on appelle A la longitude de la terre, ou 
celle du soleil augmentée de »8o°, au moment de la seconde 
observation , on aura 

M — R cos A , N = R sin A. 

Appelons encore pour un temps quelconque t : 

• l’excentricité de l’orbite terrestre , 
e la longitude héliocen trique de la terre , 
y la longitude de l’aphélie , . 

V ou R le rayon vecteur , 

on aura par les propriétés du mouvement elliptique et en 
regardant y comme constant, 


R * 

✓(*-«•) 
R* ' 


1 — t CCS {* — y) 


* sin (* — y) 
✓(*—*>• 


Or on a X = R cos * , Y = R ain * ; de-là on déduira les 
, , dX dY , 

valeurs de —r— et — r— , lesquelles en faisant t=o ou * =A, 
dt dt • ’ 

deviendront celles des coefficiens M' et N' : on aura donc 


— sin- A -t- 1 sin y 
✓ (»-•*) 


cos A — • cos y 
t/(‘ — f *) * 


XII- Si au lieu de la constante y on introduit l’anomalie 
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vraie f = A — y , ce qui donne y = A — « f ; on aura , en 
développant ces expressions jusqu’aux i* inclusivement , 

M ~ — sin \ ( i — i coS 1r -J- p « cos A sin t' 

N' = cos A ( l — • cos 1 r + r •’ ) — t sin Asinl'. 

Par un gemblahle développement on a 


R = » + • cos +■ — «* sin * ÿ 

= 1 — t CCS * + t* ; 

H 


( *9) 


(ao) 


f en 

( 31 ) 


De sorte qu’on peut écrire ainsi les valeurs de M' et N' , 

M' = — sin A ^ — • sin * cos A 

N' = + cos A ^ ^ — • sin ï sin A j 

d’où résulte 

M* + N" = -g i 

MM'+ NN'= —tsin *.R. 

On conserve dans ces expressions lequarré •* de l’excentricité, 
quoique ce terme soit très -petit , parce qu’il ne complique 
presque pas le calcul, et que si on l’oinettoit entièrement , la 
valeur de R se réduirait à 1 + t cos ÿ, et ne correspondrait 
plus assea exactement au logarithme donné par les éphémé- 
rides. On peut supposer pour les temps voisins de celui-ci, 
• = 0,01679, alors on a •*= 0.0003819 , et 

% ( * R ‘ ’ ) = l »g ^ — o. 000061a. 

A l’égard du terme • sin ’V, on le calculera en prenant dans 
les tables le lien do l’aphélie y qui est très-peu variable, d’où 
l’on conclura l’anomalie ’fr = A — y. Ou peut encore déter- 
miner la quantité isint avec une précision suffisante par la 
seule valeur connue du logarithme de R ; car tirant de ce 
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logarithme le nombre -^ = i + « , on aura « cos'ÿ = V — u , 
K 

et de -3 à 


tsinf = =fc */[(« + ** — b) («-*-•* + b)]. (aa) 

Quant au signe de cette quantité, il sera positif depuis le pre- 
mier juillet jusqu’au premier janvier , et négatif dans les six 
autres mois. 

v 1 - |* 

Enfin on peut encore observer que les quantités — ^ — et 

n. 


• sinf, qui entreront comme coëfRciens dans nos formules, 
sont liées entr’elles par cette relation : 



•• tin * ÿ 



XIII. Il s’agit maintenant de former les coëfficiens des équa- 
tions à résoudre, de la manière la plus simple et la plus com- 
mode pour la pratique. Four cela, il convient de remettre à la 
place de /, g,f°, (f , &c. , leurs valeurs on a, b ,.a° r b° t &c. 
données immédiates des observations. Cette substitution étant 
faite d’abord dans la quantité a, on prendra pour abréger, 

D = tang b' tin (a— d”) + tang b • tin (a — a) • 

‘ + tang b tin (a* — a'), (aî) 

I jj 

et on anra a = —y. Faisant ensuite les mêmes 

cos a cos a cos a 

substitutions et celles de M=RcoaA, N= R sin A, dans les 
valeurs de B", B, B', on aura 


— = ^ C ^ ~~ f tang b' sin (A — a’) — tang b * sin (A — o')] 

A 1/ 

• • » * • ’ » . 

[tangb’sin (A — a) — tang b sin {A. — a*)} («4j 

B* llc<wa° , . . , . ,, ,, . 

— = — — — [tang b sin (A — a) — tangb sin (A— a )J. 

A U <14. • 


Digitized by Google 


(» 5 ) . 


( i3 ) 

Soit encore pour abréger , 

C = tang b' sin (A — a") — tang 6* sin ( A — a') , 

R w 9 * 

ia première des équations (17) donnera /*= — — C cos a; or 

aD 

parles triangles rectangles TOC, TKG (fig. i), on aTK= p cos b, 
èt TG <n 1 /* = pcos b cos a ; donc en introduisant /> à là place 
de p, on aura 

R-**C - RjS’C /j_ t \ 

* aD cos b al) cos b Vr 5 Rv ( a ®) 

Substituant de même l’expression de /* en p, dans les valeurs 
de m, n,p , des équations (13), on aura * , 

. m — p cos a . cos b + R cos A , . 

n = p sin a cos b + R sin A (27) 

p —f sin b,. ... ... .» 

et l’équation m?+n‘+pP=si 1* deviendra . .1 

r*s=R*-K*Rrco« b cos (A— o) + p*. -, 

Celle équation seroit donnée immédiatement par le triangle 
CST; car en appelant c l’angle entre le soleil et la comète, 
les côtés compris sont R et p , et le côté opposé r, ce qui donne 
r* = R* — a R p cos c + p‘. 

Or l’angle c a pour mesure l’hypoténuse d*un triangle sphé- 
rique rectangle dont les côtés sont b et 180®— A + a, on a 
donc cos c ~ — cos b cos (A — a). * 

XIV. Les coëfüciens C et D se déduisent immédiatement des 
données du problème par les formules (s* 3 ) et (a&) ; si ensuite 
on lait 

R«*C 

(»*) 


cos c — — cos b cos (A — al, h = 

, Àt a ücosb* 

on aura les «eu*, équatious 


t 


»(Wi 

r'cs R * — p (aR cos c) + p‘, 


(* 9 ) 
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dans lesquelles il n’y a d’inconnues que les deux côtés r et f 
dn triangle 8CT. 

Si on élimine f de ces deux équations, on aura une équation 
en r du huitième degré , mais qui sera divisible par r — R , et 
s’abaissera immédiatement au septième. C’est le résultat auquel 
sont parvenus tous ceux qui bnt traité avec succès le problème 
des comètes ; et ce résultat est d’autant plus remarquable , que 
les équations (39) se déduisent assez facilement des données 
de l’observation, et qu’elles ont lieu sans supposer que l’orbite 
de la comète soit une parabole. 

XV. Comme les quantités r et c sont essentiellement posi- 

tives, il n’y aura d’admissibles, parmi les diverses solutions 
des équations (ag) , que celles qui donneront r et p positives. 
Or en construisant les deux courbes de genre hyperbolique 
qui représentent ces équations, et examinant les diverses inter, 
sections dont elles sont susceptibles , on parvient à cette con- 
clusion générale sur le nombre de solutions utiles que ces 
équations comportent : ... 

Si ih cos c est positif et plus grand, que R 4 , les équations (ag) 
admettront une solution et n’en admettront . qu’une ; dans les 
autres cas , ces équations auront deux solutions ou n’èn auront 
aucune. 

• • * .* 

XVI. A l’inspection de l’équation p = A ( on voit 

, ... ? » ’\r JA. 0 yr • 

que le signe A doit être le même que celui de Doné 

si h est positive , on aura r < R , .c’est-à-dire que la comète 
sera moins éloignée du soleil que n’est la terre. Au contraire, 
si A est négatif, on aura r>R. Dans le cas où l’on auroit 
D = o, A seroit infini et ou auroit r — R, c’est- à dire que la 
comète et la terre seroient à égale distance du soleil. Ce cas est, 
iTqprès la remarque de Lambert , celui où les trois lieux 
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upparens de la comète seroicut situés dans le plan d’un même 
grand cercle ; car pour que trois points dont 1rs longitudes sont 
a°,a, a, et les latitudes b a , b, b', respectivement , soient situés 
dans le plan d’un même grand cercle, il faut qu'ou ail l'équation 

o— tan g b' sin (a — a°) + long b’ sin (a' — à)+tang b sin (a* — a'). 

Donc lorsqu’on aura D = o, ou seulement D très-petit, on 
pourra en conclure, ou exactement, ou au moins par approxi- 
mation , 

r = R et f = a R cos c. * 

XVII. C’est nne circonstance peu favorable que celle où l’on 
tombe sur des valeurs de D très-petites ; car alors une erreur 
très-possible d’une ou de deux minutes sur quelqu’une des 
quantités données par l’observation, pourroit changer dans une 
proportion notable la valeur de D,ou même lui donner un 
signe contraire"* celui qu’elle doit avoir. Dans ce cas , il n’y a 
aucune conclusion certaine à tirer des équations ( ag ) , sinon 
qu’on a à-peu* près r = R et p = a R cos c; on se dispensera 
donc alors de chercher les racines de ccs équations , et on aura 
recoars aux autres équations qui seront données ci-après. 

XVIII. Dans les autres cas, la résolution des équations (ag) 
s’effectuera assez facilement par les fausses positions; niais pour 
cela il est bon de conuoilre d’avance les limites de r et de p. 

Lorsque h est positif on a r < R ; mais la seconde des équa- 
tions ( 39 ) donne r* = R* sin' c + (p — Reojc)'; donc on a 
en même temps r>Rst/»c. Ces deux limites auxquelles on 
peut joindre celles de p, savoir p > o et p<aRcoac, sont 
suffisantes pour diriger le calcul et rendre la résolution assez 
prompte. 

Si h est négatif et cos c positif, soit h = — » ; ayant déjà 
r>R, on doit avoir par conséquent p >aR cos c; mais de 

l’équation ^ — j , on tire p < ; donc on aura 


on 


- , aR i co» c i' \ ' , , , * ' 

r<.ÿ IR — ' + I* laut observer que si 

n’avoit pas i > 2 R 4 cos c, les deux 1 i ni i (es de p scroient incom- 
patibles, et il n’y auroit pas de solution. 

Enfin si h est négatif, ainsi que cos c; soit h rc — », et 

cos c = — l, on aura toujours p < , et l'autre limite sera 

simplement p > o : celles de r seront en même temps r > R 
2R»7 i‘\ 

' , '<< / ( R + -Rî- + ff> '■•••. . 

XIX. Il faut maintenant faire la substitution' des valeurs 
de mf, n'y p , dans J’équation (16) : or en vertu des for- 
mules (ia) , (17) et (»♦), on a 

, R» J f cos a° tang b sin ( A — o’) — cos a° tang b' sin (A — o)) 
+ aD j+cos a' tang b sin (A — 0°) — cosal tang b‘ sin(A~- •«)} 

- î * r 4 ' 

R« 9 f sin a,' tangb sin (A — a') — sin a? tang b' sin (Jk—o)\ 
aD j+sino' tang b sin (A — a") — tin a! tang tf sin (A— -a)f 

R»6 f tang b’ tangb sin (A — a') — tangb°tangb' sin{ A — a) 

- aD ( + tang b' tangb sin ( A — 0°) — tang b' tang b° sin(A — a) 

Appelons F, G , H les quantités qui multiplient — jj- dans les 


( 3 o) 


valeurs de m',n\p' respectivement; nous aurons en substituant 
aussi les valeurs de M' et N' : 


« A «, . r 

J ■ 


ro' = 

• 1 . 

— sin A ( 

s**\ 

■ 

1 — « sin T cos A + — =- F 

2U \ 

s. a ) 

yj 

co s A | 


R m fl 

i c crn "4' ctrï À I R- 

Tl 

L H ) 

1 » Slfl ” Sin fi y rv U 

r aD 


d’où résulte 
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a R*9/i— -f«*\ 

R 1 ~jj~ ^ g — J (F sin A — G co* A) 

R- fi 

~-ÿ-.iii»Ÿ(F WJ A + G sin A) 

R***9’ 

+ ~4D î ~ (F * + G ‘>- 

Donc si pour abréger encore lea expressions on fait 

r^,tmngb .in (À-«*) «m (A-o' ) - lang V .in (A -<f) «n (A — a) 

— 4» «û» (A — a') «ut (A— o) 

* * ^ 

Q= A»* -à-v 4< «« (A— a*) mn (A-o) 

— tangfrco. (A— a' ) «ûi (A-a) 

H = _ i (<,„* 4* /an* 4* «ta (A-a) + tang b tang 4* „* (A-a > 

, + long b tang b' .in (A — O*). 

Les trois quantités m', p', pourront être mises sous cette 

forme : * 

; * . . . . c3») 

et leur substitution dans l’équation (i6) donnera 

3 » 3 R«i/l-I,' V 

; â~P+'sôiŸQj 

+ “ÏD rCP ‘ + Û' +H *)’ («) 

XX. Si l’on suppose que l’orbite est parabolique, on aura 

3 
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- = o , et l'équation précédente devra être combinée avec 

? « . 

l’équation 


i. 

7 


• i. 

R* 1 


(34) 


d’où l’on voit que la résultante en r ne sera que du sixième 
degré. 

Celle équation aura toujours une racine réelle positive ; car 
en faisant u = x, - ~ y, si on construit les deux paraboles 
qui ont pour équations 

■ s'= .! ■« — ’ — à'x + B'x* 


R 

y* — x + 


R’ 


Iæ première sera une parabole ordinaire , ayant son axe 
dirigé dans le sens des ordonnées. Cette cOurbe est située toute 
entière du côté des y positives, et ne coupe point la ligne des 
abscisses, parce que le second membre de l’équation ( 33) formé 
de In somme de trois quarrés , ne saurùit devenir nul et a 
ses facteurs imaginaires. D’ailleurs lorsque x = o, l’ordonnée 


f • 1 . 

y *■ T, ou a- peu près y = f 

‘XV * * 


/• 


t La seconde courbe est une parabole cubique dont une branche 
s’étend à l’infini dans l’angle des x et y positives,' ët l’autre 
dans l’angle des i ety négatives. Lorsque x — o , d’or^oHnéç 


— ou à-peu-près . y == i ; ainsi dans çcltc courbe, les 

*V.O . i ‘I " * — ’ 

innées sont plus grandes vert l'origine des abscisses que 




ordonnées sont plus grandes vevs l'origine «les abscisses que 
celles de la première parabole. Mais comme à une grande dis- 
tance de l’origine c’cil le Contraire qui a lieu , il s’ensuit qu’il 
y aura nécessairement une intersection ; et il ne peut y en avoir 
qu’une y parce que la parabole’ ordinaire appose sa convexité à 
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la ligne des abscisses , tandis que la parabole cubique lui oppose 
sa concavité ( voyez fig, 4). 

Le système des équations (33) et (34) est donc préférable à 
celui des équations ( 39 ), non-sculcnient parce que le degré de 
l’équation finale est moindre, mais encore parce que la solu- 
lion de celles-là est toujours unique , et ne présente aucune 
nmkiguilé. 

Il Lut observer néanmoins que l’équation (33) est sujette 
au même inconvénient que la première des équations ( 29 ), 
dans le cas où D est très-petit, parce qu’alors les erreurs des 
observations ont une trop grande influence sur la valeur de D; 
et qu’ainsi la résolution des équations (33) et (34) pourroit 
conduire à des résultats défectueux. 


XXI. Pour obvier à cet inconvénient , il faut éviter tout 
emploi du coefficient D, et on en a heureusement la facilité, 
puisque le problème des comètes , dans la supposition d’une 
orbite parabolique , offre une équation de plus que d’inconnues. 

ùi 9 f cos b 

Par l’équation (a6), on a el celte valeur 

étant substituée dans l’équation (33 ) , on aura, en suivant tou- 
jours l’hypothèse parabolique, c’est-à-dire eu faisant - = o, 

a 



a t cos b ( 1 — { «* 

4C \ R P + ‘ 

(P*+Q*+ H*); 


,»V. Q ) \ 


(36) 


équation qui devra être combinée avec l’équation 

r* = R* — p ( a R cos c) + p* , (36) 

It . * ■ -* . J- t **-*#•’ * ■ 

afin d'en tirer les valeurs de r et f. 

Si on éliminoil r ou p de ccs deux équations , l’équation 
finale qe monteroit encore qn’au sixième degré. D’ailleurs les 
cjçtjEcicns de ces équations se.dcduiscut, aveu une exactitude 


( »o ) 

suffisante, des données de l’observation. Le coefficient C qui 
sert de diviseur, est celui qu’on connoitra avec le plus de pré- 
cision 3 car en regardant S comme une quantité du premier 
ordre, C est du premier ordre aussi , et les trois autres P, Q, H 
du second. Il paraît donc que la combinaison des équations (35) 
et (36) réunit tous les avantages qu’on peut deairer dans la 
solution du problème des comètes, et c’est à cette combinaison 
qn’il convient de s’arrêter définitivement 

XX IL On peut prouver encore que les équations (35) et (36) 
auront toujours une solution en nombres positifs. 

En effet , si on décrit les deux courbes qui répondent a ce* 
équations, en prenant p pour l’abscisse et r pour l’ordonnée; 
le lieu de l’équation (36) sera une hyperbole équilalère dont le 
centre est snr la ligne des abscisses à la distante. R cos c, et 
dont le demi-axe, dirigé suivant les ordonnées = R «u c(voytt 
fig.5). 

L’éqnatkm ( 34 ) , qu’on pourrait réduire k la forme 
y (x‘+ a') — a'b, appartient & nne courbe qui, comme la 
conchoïdc supérieure, a pour asymptote la ligne des abscisses , 
et est située toute .entière d’un même côté de cette ligne , sa 
plus grande ordonnée étant b. Lorsque f — o , l’ordonnée 
de l'hyperbole équilatère = R , celle de la seconde courbe , 

= j-t=iRà peu- près; donc puisque la première s’éloigne 

à l’infini de l’axe , tandis que l’autre s’en rapproche de plus en 
plus à mesure que l’abscisse est plus grande, il s’ensuit que les 
deux courbes auront nécessairement une intersection dans lo 
sens des abscisses positives ci n’en auront qu’une. Ainsi la 
résolution des équations (35) et (36) aura encore l’avantage 
d’étre toujours possible et de n’offrir aucune ambiguité. » 

Pour avoir les limites de r, il faut observer que les seconds 
membres de» équations (35) et (36) ne peuvent se réduire à 
zéro, et qu’ils ont tous deux un minimum. Le minimum «1er, 
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dans l’éq nation ( 36 ) , est R ain C. Soit Z le minimum du second 
membre de l’équation (35), et il est clair qu’on aura à la fois 

2 

r>Rsin c et /••< •=-. Au moyen de ccs limites, il sera tou- 

£j 

jours assez facile de résoudre les équations dont il s’agit par les 
fausses position». 

» jl- u * . r ^ v . » i t*j 

XXIII. L’application de nos formules exige que le temps 9, 
donné en jours et fractions de jour, temps moyen, soit converti 
en arc du moyen mouvement du soleil , cet arc étant lui-même 
réduit en parties du rayon i.Pour faire cette double réduction, 


a * 


i 


il suffira de multiplier 9 par jgg jggj g > c’est-à-dire qu’il fau- 
dra ajouter au logarithme de 9, compté en jours du temps 
moyen , le logarithme constant 8.a5558ai. 

En faisant abstraction des erreurs des observations, les coëffi- 
cicns qui entrent dans l’équation (35) ne doivent être censés 
exacts qu’aux quantités près de Pordrc 9* { d’où l’on voit 
qu’il faut choisir trois observations faites dans un intervalle 
de temps pen considérable, et qui n’excède pas quinze à vingt 
jours. Mais il ne faut pas non plus que les observations soient 
trop rapprochées, parce qu’a lors le mouvement géoccnlrique 
de la comète seroit trop petit, et les erreurs des observations 
auraient une trop grande influence sur les résultats. 


XXIV. Après avoir résolu lescqualions (35) et (36) qui 
feront connoîtrc les quantités r et p, il reste à déterminer les 
élcmens de l’orbite; mais comme ces élémens, tels qu’ils résultent 
d’une première approximation, ne peuvent pas être fort exacts, 
et qn’il y aura toujours lieu à les corriger par le moyen d’obser- 
vations plus éloignées entr’ellcs , on pourra se borner le plus 
souvent à chercher une valeur approchée de la distance péri- 
hélie , et de l’instant du passage au périhélie ; car avec la con- 
noissancc approchée de ces deux élémens , on est en état de 
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( aa ) 

procéder nu calcul nécessaire pour obtenir une orbite corrigée. 

Or Si on appelle « la distance périhélie , et qu’on calcule k 
par la formule 

k — mm + nn + pp'.i 


comme cette quantité est égale à 


rdr 

~d7' 


si elle est négative, 


on saura que la comète s’approche de son périhélie j cl si clic 
est positive, on saura qu’elle s’en éloigne. De plus, au moyen 
de k , on aura la distance périhélie 

n (3 7) ) 

Cette distance étant connue , on déterminera l’anomalie vraie 4 
de la comète au moment de la seconde observation , par la 
formule 

COs' j 4 = (58) 

. r 


Enfin l’anomalie étant connue, on cherchera dans la table 
des comètes le temps T qui répond à cette anomalie : faisant 
ensuite 

. . 1 

« \t t = ïi*.T, (3q) 

on aura le temps t employé par la comète à parcourir l’ano- 
malie 4- Ce temps étant ajouté à l’époque de l’observation 
moyenne, si la comète avance vers son périhélie, ou en étant 
retranché, si elle s’en éloigne , donnera l’instant du passage au 
périhélie. 

Tout se réduit donc à déterminer la quantité k, dont la 
grandeur et le signe sont également nécessaires à considérer. 
Or en substituant les valeurs données par les équations (3a), 
on trouve 

t 

> cos 


(f 

k f cas b (A— a) 1- 

+ (E. — f cos c ) Q — * «'» Ÿ )- 

|.,t» jm * j» ; V*V i - ' t.,*'. :! 


b „ (j — ï*‘)\ . . f cosb 

- p --îr-) +f ™ 4 --îc- H 

t» ; *. ’t 


(4o) 
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XXV. Mais si on vent déterminer à la fois tous les élémeus 
de l’orbite , il sera inutile de chercher partiellement la vateur 
de /t.et il faudra calculer les valeurs des six quantités m , «, 
p , m', ri, p\ par les équations (37) et (37). 

Ces valeurs étant connues,, on en déduira immédiatement 
celles des trois quantités mri — m'n, mp' — m'p , np — n'p, 

n* dp 

lia première étant égale a -^-7-.-^, fera connoître par son 

signe si le mouvement de la comète est direct ou s’il est rétro- 
grade. Dans le premier cas, mri — m'n doit être positif, et 
dans le second il sera négatif. 

' • . ° , J-l t 'K 

Appelons maintenant! l’inclinaison de Torbito , et S la lon- 
gitude de celui des deux nœuds que la comète devance dans 
l’ordre des signes. On aura', par les propriétés du mouvement 
parabolique, les trois équations : 

mri — m'n — =t cos I p/ïn 

•. mp' — m'p = =t sin I cou S /ïn (4r) 

np' — ri p = ± sin I sin S p/Tn. 1 . . , 

Dans ces formules qui doivent toujours donner I < go° et 
p/alT positif, il faut prendre les signes ambigus du second 
membre tous positivement, lorsque mri — m'n est positif , et 
tous négativement, lorsque mri — m'n est négatif. 

La longitude S en particulier sera déterminée par la formule 

( 4 a) 

Mais cette formule, qui convient également aux angles S et 
180°+ S , donne la position de la ligne des nœuds, cl n'indique 
pas l’un des nœuds plutôt que'fautre. O11 achèvera de déter- 
miner S en calculant l’inclinaison I par la formule 

\1L- . . 1 . n I VI ! • I ■•*”1 •',! ' • I • .'A .IV i f* 

tang I = 7— — — j -~- — — ; ( 4 $) 


, c np — np 

tang S = — ; —r . 

mp — m p 


.*» 


•**"6 — 7 / ; v c > 

(mn-m.n)cosb 

et prenant des deux valeurs de Scelle qui rend ian&l positive. 


( »4 ) 

Si cos S éloil très-petit, U sentit plus exact de déterminer 
tang I par la formule i 

tang I — — — St O^) 

° (/un — m n) sia £> 

et on lèveroit toujours de la même manière l’indétermination 
de S. 

Le nœud ainsi déterminé par la valeur de S , sera le nœud 
ascendant, si le mouvement est direct et la latitude boreale, ou 
si le mouvement est rétrograde et la latitude australe. Dans Us 
deux autres cas , ce sera le noeud descendant , tt il faudra 
ajouter 180“ à S pour avoir le lieu du nœud ascendant. 

Les angles S et I étant déterminés , oa aura la distance 
périhélie n par la formule 


( mn — jw»')* 

n a cos' I ’ 


on pourroit aussi la trouver indépenda mm ent de ces angles 
par la formule 

2 n = {mn — m n}' 4 . {mp‘ — m‘p)' + {np—np)', (45) 

qui résulte évidemment des équations ( 4 i); et il est aisé de 
•voir que cette dernière râleur s’accorde avec la formule (37). 
Car en général le produit (m' + n'+p‘) (m'* + n'* + p‘) se 
décompose en ces quatre quarrés : 

(mn'+Bn'+fp')‘+ (mn — mn)' + {mp— m’p)' + {np'—rip)'. 

Donc U somme des trois derniers quarrés = (m* + n* + pé) 

(m* + n ‘ + p‘) — (mm’+ nn’ +pp )* = r* . ^ —£• = a r — i I*= a n. 

* # . * S % 

XXVI- Avec la distance périhélie n et le rayon vecteur r, 
on calculera l’anomalie Vraie 4 de la comète au moment de la 
seconde observation, et de-là l’instant du passage au périhélie, 
ainsi qu’on l’a expliqué dans L’art. XXIV. 
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Enfin pour avoir le lieu ilu périhélie, on calculera la lon- 
gitude héliocentriqne ♦ de la comète par la formule 


mÊ 


1 4' 


tang 9 = 


ut’ 

i 


WM 


C«) 


ayant soin de choisir entre les deux angles ? et p -b i8o° qui 
conviennent également à cette tangente, celui qui rend sinp de 
même signe que ri, ou cos a de même signe qrtc m. 

♦ étant ainsi oonnu , on trouvera toujours que p — S est 
positif et pins petit que i8o°, conformément à la supposition 
faito que S est la longitude de celui des deux nœuds qui est 
moins avancé que la comète dans l’ordre des signes. Cela posé» 
soit r la distance de la comète au nœud dont la longitude est S, 
on calculera r par la formule 


tang (» — S) 


(4 7 ) 


tang r = 

co s I 

et le lieu du périhélie sur l’orbite sera 

« , ’ ! 

1 1 i S -f «• + 4 si la comète marche vers son périhélie , 
ou S + r — 4 si elle y à déjà passé. „ f 

_ ' ’ i _ __ " 

xxvîf. Pour éviter la peine de cTièrclier ailleurs la démons- 
tration des formules ( 4i), nouà croyons devoir l’ajouter ici. 

Des équations différentielles (i) on déduit immédiatement 
les intégrales suivantes , qui expriment que les aires décrit»* 
dans les diverses projections de l’orbite sont proportionnelles 

au temps ■ " '‘'rit » k ! 

xdy — ydx = C ’dt 
xdz — zdx — C'dt 
' ' ÿdz—zdy = C 'df, 

on en tire d’abord l’équation 

U > I •» , I ! Î 

C'a — C'y -f C"x = o , 


(48) 

. • 

i le iaq 


rtt na’b 

_ . _ ^ . . (Il r ( ^T. 

qui prouve que la trajectoire est située dans un mè«ic plan. 
En faisant x = o, cette équation doit donner pour la ligne 


«<*64 


v 

'tang S. Dp 


.., (taris léfyfcft^ès*;êt^?S*<i indéfini « dont 

l’origine soit sur l’axe des x ; les coordonnées du plan qui 

» ‘i » : («i _ .4 i 


On anra donc aussi 


• A 


r. 


de 

ds' 


CT 

c - 

x = — - 6 in s 

— coi <* 

S V -DP - • 

• . 

<r \ 



é - -4 ' ... T 

,■ ■ r i 

. 1. r 


-i. _ ■’* 


ob 


ton ou tang I : on à donc 

^T^gcosS + ÇsnsS^^; 

donc C° = C cos S tang T et C“ = C sin S tang I. 

Cela posé , les équations (48 ) rapporté» à l’époque où l’on a 
«•t/*») &c - deviennent , 1LC< . xf tô'jrfŸ 
*' *’■ ‘ ‘‘ mn'->É'a a» C' *1' * , .. ., 

mp — m’p = C' co* 8 tangï 

np' — »'p = C' «in S issjf I , .. t 

et il ne veste à déterminer que C'. Or dans l’équation (16) qui 
peut s’écrire ainsi : . , 

. ù . i£± d S+ d *' _ j ^ 1 \ 

r - • «&* T M a’ 

» v * 4 * r^cf# 11 "* 

le premier membre ,= — — — j multipliant de part èt 


d’autre par r*, et mettant à la place de ^ sa valeur k, on 

... . - 9 *:< 

déduit 

’ ' i*â't y r*-‘- ‘,'1 **'**'■ ..-i 

’ •* 1.; *.v'o* L 


on en 


i • 


d«* 
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Mai* l’aire j (xdy — ydx) est la projection sur le plan de 
l’écliptique, de l’aire {r'dv décrite dans le plan de l’orbite j 
donc puisque l’inclinaison mutuelle de ces deux plans est I , 


lu j”' i *i 

dt 


A* 


■ d y — ydx 
dt cos I 


CT 

N H. t»sp4'! 

^71 

\ 

»d igp%>.l 

-*>■ 

’F» • • * 1 y 


Donc 

o« 

MM ib ac*s Wsiy .1 JS 

Faisant A la fois r — n et * = o , le- second membre ac réduit à 

cos' I fa H — — J, Donc on a en général, 

-t> va/ 

iuav fv> fai] 

«9-fieui * 


C=± cos I 

■ 


... .UKj J t. iM 

«T' 

n — — i ii> i * t 


et en particolier , lorsque l’orbite est parabolique, m l. 

C ' = cfc cos I a IT ; ' • ‘ ' 


ce qui donne les formules (4r). 


XXVIH. La solation que nons venons de développer ] 
le cas où tes observations sont faites à des intervalles de temps 
égaux, paroît avoir tonte l'exactitude et la simplicité qo’on 
peut désirer dans une question dont la difficulté est générale» 
ment reconnue. Si par des circonstances peu favorables cette 
solution ne conduit pas 1 des résultats suffisamment exacts , ü 
faudra s’en prendre n l’im perfection des observations et non 
anx négligences que nous nous sommes permises dans la redno 
tion dps formates. En vain voudroit- on pousser plus loin 
l’exactitude de ces formules (ce qni ne paurroit se faire que 
par des calculs trop compliqués pour être de quelqu’utilité 
dans la pratique); les erreura des observations domineraient 
toujours et empêchcroient les résultats de passer un certain 
degré d’approximation. D'ailleurs, ou sait assez que trois obser- 
vations prises à peu de distance les une? des autres, ue peuvent 


pas déterminer nssex exactement l’orbite d*une Comète* ét qu’D* 
font totijburs recourir à des observations pins éloignées pour 
rettificr'les élémêns déterminés par une premier» approxima- 
tion. Il est donc inutile de pousser le scrupule trop loin sur cette' 
première approximation, et l’objet doit être suffisamment 
rempli par la méthode que nbus avons exposée'. 

Lorsque les trois observations données ne seront pas également 
distantes enlr’elles , il faudra, comme noua l’avons déjà dit, 
calculer par l'intcrpolalio'n un nouveau lieu qui soit avec deux 
des autre» dans la condition requise de l’égalité des «lislancesi 
Dans ce calcul d’interpolation ,.il ne faudra faire entrer que les 
trois observations données , ou en général les trois observa- 
tions connues les plus proches ilu,momcnt pour lequel on veut 
déterminer le nouveau lieu/ Tels sont les préceptes au moyen 
desquels on obtiendra tpujours une solution suffisamment appjq:, 
ebée dnns la pratique.^ , ^ } 

XXIX. Mais à considérer le problème soîis un point de ?ue* 
purement analytique , noos n’aurions rempli qu’irnpaf folle- 
ment notre objet, si nous ne gisions voir par quels moyens; 
an, peut obtenir une solution générale indépendante de la sup- ( 
position de; 8 ss 8. n <■ : rjiç,!, Ijj kj 

Dans «Ho vue, nous allons reprendre les équations (s ^ ) , et" 
dOu* observerons d’abord que les unîmes raison^pa,r desquelles , 
nous avons supi>riulé les. termes afféclés de 8' dans les forinplqaj 
de. l’art- ,X* doivent également faire omettre, (fon? les vaijjqr^ 
tle^'let le». termes affectés du facteur 8* — 3 5'+, fl Il nç reflje f 
alors, daus-Jes équations (| s) que deux., sortes deviennes , lc»j 
ups affectés du faotenr »— (8' — 3), les autres affectéf, du ( , 
facteur ♦>(9 — 8) ; or bienfait 

■ u « — ^(.S— .(49); 

leSecoml facteur ** (3‘ — 3) pourra être représenté pai^ (3’ — 
càr la différence de ce8 deut quantité*', qui e»t>£(*' — ♦ jt'f W 
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trouve de l’ordre des quantités qu’on a droit de négliger. 
Il suit de-là qu’on peut mettre les équations (u) sous celte 
forme: «* >• 

:y.- 


<•' - -*(&-'*>+ hr-'*\+ w-üh 

’ = X (IfE'—tJj E’ ) +;{ (l'_l)N 


* - <**-,»*** 


4 P 

», 


•*‘4 


*r 


,*U 

MàlJ 


Soit ('= x'+y et 0 = x — y,- ces équations deviendront : 

, «r 

tjj» 

** **. ' (B +:f^E)-‘ 1 f t ■ ■ 1 r • "t ■ .XJfX 

• • MÙin u.iuoi: jtr.aoi 4i . jy* , jpu,>L> , no 

f ÇX ' - * - - - - - 


'= y- [B‘ B' + j/(E--E.;]+ 


U 




_ÎT[B- + B'+i/(E- + E')] 

•j w vt 

- ,f - [ P B»-/' B' + * ÿ\f E ° — /' E')] + |y N 

' » ffcj Î WÎH -tw " ' 


21 A 

“[/’ B ‘+/ ( 





*»•'■'” e • , «ct>i,il u i rtiiint ;» es tus e-»!» Inwij'l A 

. ■ • + ff'B' + |y ^yE ‘4 g-'E'Oï ,d u:/lno,?vJ 

^ \¥ 9 Urtf%, , 

Lértroi* dernières expressions se simplifient encore au moyen 
des valeurs de B’ -f- B' ( f° B* +- y B', &c. données par les équa- 
tions (16» il en résulte V?!ti<sjPW|* 


iy(/"E*+/'E')] 


l ‘ = , ;»V < #4 .^f) • ,,, 


ï* r= 


1£ 
a A 


[B* — B' + jy ( E*— E* ) ] 


:.n . ■! 


+ LZ_(B+}/E) + H/M‘ ' 

a a 

[/» B‘ — / ' B' -h-^' (/• E* — / E") ) 

84 • , a 


. , M. , 

J>v 

• - - ✓ «■ , , • > * 4 

* XXX. Maintenant la substitution de» vateu» de fifif ' > *®' '• 
en a, i, a, Ace., donne, suivant les dénomination» employée» 
dans les arL XIII et XIX , ff — 

‘ ’ - = Rm oæ • ' ' 

D •• t ». ■> »>, 


* i. 

B*-»' F 

=* #• «T 


■s» a ' ' , - 


L 1 ( 4 * 4 ' 


A-~- 

r*-r »' _ R & ' 

,x0it-y\‘ -f 4 -“"S ' \ } — — 

£ï^£ -,• a.® ;v 

. pst.-* J», ■ s.rt i t ay * *’ - ' 

A l’égard des antre» quantités composées en E , comme celles-ci 
le sont en B, il faut observer que £ n’est .mire chose que B 
dans lequel on met M’ et N’ à la place de M et N. Or en 
négligeant l’excentricité « ( ce qui peut se tains sans inconvé- 
nient, parce que tous les tenues «fiertés de £ sont multipliés 
par la quantité très-petite y 1 ) , on a M’ — i l 1 eût A et N* — cos A , 
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ou encore, par la même raison, if = — R sin A , N'= R cos A. 

Observons de plus , que chaeune des quantités C , F, G, H , 
considérée comme fonction de A,e»tde la forme «.'cosA+Csin A, 
•' et C ne contenant point A; donc puisque pour passer des 
fonctions B aux , fonctions E, il suffit de changer ca» A en 

sin A et sin A en + cos A , il s’ensuit que si la fonction 

g itfnlr. f eTî qtnT. etT , ' 77 

— ou R cos a. — , par exemple, est représentée par «'cos A 
4- C sin A , la fonction )y' . — le sera par f ÿ\—**in A + Ccos A), 


et la sonnne des deux — (B + jy E) , le 

A 

mm ^^m 


sera par 


, « (cos A — $ y sin A) + f ( sin A 4- jycosA), 

ou ce qui revient au même, en négligeant les secondes puis- 
sances dey par ’* 

«' cos(A + }>') + C' sin (A + y y-), 
d’où il suit que celte fonctioa ne seça autre çhbsc. que Ta quan- 

tilé — ou Rcoso.rpr, dans laquelle au lieu do A on mettrait 
a JJ 

A+~y,- la partie jy' ou y(,®'— ®) devant être exprimée en 
degrés , minutes et secondés du moyen mouvement du soleil. 

a r ^ 1 1 ' 

XXXI. Un résultat semblable aura lieu pour les fonctions 

analogues de B et de E qui entrent dans nos équations. Soient 
donc 

C.., Fr, Gr, II. , 

ce que deviennent les quantités C, F, G, H, lorsqu'on y mat 
A + \y k la place de A , et nous aurons enün, après avoir mis 
P cos a cos b k la place de p* : 

p co» g _ R#r ' ?- 


5 , 

!L — — — g i f — — v, cos a + — — 

#trN P" < aüi ‘^V^ a 


Rx _ Ry’ Ry 

— Fi + — — Ci cos a -f cos A 


((-3a)) 

4 = d%^»ïh *?••«* IP 

«J. f a-D a D >„iu a-K i i iuU,, 

./. n' R#' R /■ _ I , • • IjLiaOU * 

T = Td H,+ 75T Cito ^ i - 

/. V-» » lU I fifW i » I; M njan-.p 1 X^MI K ' A 

s choses étant réduites à cel état de simplicité, on calculera 
trois coëfficiéns F»*, Ga , Ha , d’après les valeurs 


Les 

les 


i. tî> ‘<1 vAî; • t*> .alsirtTI Vf , -.nwrttt »:•» 

_ - cojs b ; ‘ .i, 1 'n. . . i. 

Fa = - , - (t Fi + .y Ci cos a -f y Dcos A) 

• * 9 • -i ; , . • -.i * -.ri,/ • 

„ cos b , , _ , „ 


Gî = -x^rpr— (*'Gi +.y'Cr sin a +y Dsi/iA) ( 5o ) 

98 Ci j ; — //i, >h!Wk il': 

. ,B» = 


•'T 


• I I . • J> O ) ' O 


- %£ («'», +yc«<fljvi)»' J . 

*.*. ' il n, • 

dans lesquelles il faut se rappeler qu'on a*'=£(8'+(), 
y = J ( 9' — 9 ) j et alprs on aura simplement , par l’élimination 
de l’inçonnue f , 

/ = Fa f , t = Gap » pV 1 À» f- (5i) 

J fl%rëtfîu'donc' !: "• ' " ” P rf t ). •» an - » uo - * * 

m' — M'+ Fa f . , ..V , 

»' 2 N ! + Gap v . . : ’ . (52) 

• Ha p ; 

et en substituant ces valeurs dans Péquition = -, 

i t ‘ 4 ' • ’kA H H * !i *• • T 

elle devient 


!.ji?s’ • > -V»; av '«si? H ü “Îj t* 

4' * 


- = 4" — i + a^M'Fa +N'6a) + p*(F*a + G*a + H*3), (53) 

r ' /<* ■ ■ . Ill*r i. ; '• -.••• 

•laqaeile devra être combinée avec l’équation ordinaire — i 

r* = R* — p ( aR coi e') '4- p\ *• ^ 11 •• •» * 

afin d’en tirer les valeurs de R et p. 

Il est visible que ces équations sont de même forme que 
celles qu’on a obtenues dans le cas de 8'i= 8 ; de sorte que la 
solution générale , au moins dans l’hypothèse d’une orbite 
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( 33 ) 

parabolique , n’offre guère plus de complication que celle qui 
a lieu dans ce cas particulier. D’ailleurs les équations précé- 
dentes une fois résolues, tout le reste s’achève comme on l’a 
expliqué ci-dessus. 

f,. . * '» 4^.' lit Al â I ; I -I «a 

Application à la seconde comète de 178t. 


XXXII. Proposons-nous de déterminer l'orbite de la seconde 
comète de 1781, d’après les trois observations suivantes, 
réduites pour chaque jour à 8 k >9' 4»’, temps moyen à Paris. 


Tempe de t’obe. 

Longitude. 

Latitude bor. 

Lieu du Q 

Log. K 

Novembre. 

i4 

'9 

*4 

D M S 

. D M S 

D M S 

9, 99/(86» 
9.994426 
9.994038 

o° 307 i4 |5 
a 3 06 5i 26 
d 3o6 43 ao 

4° 55 17 9 

4 39 14 48 
4* 3i 4 5s 

a3a 54 a 
*37 57 4 

a43 0 4 1 





Ces lieux sont tirés des Além. de l' yicad. des Scienc. année 
1780, pag. 67 ; et afin de rendre égaux les intervalles de temps 
entre les observations, on a interpolé les trois lieux des 19, aa 
et ?5 novembre, ce qui a donné le lieu du a4 tel que nous l’in- 
sérons ici. Ces sortes d’opérations se pratiquent fréquemment 
à l’égard des comètes, et elles sont utiles pour simplifier le 
calcul des élcmens. 

Si a°, a , ci sont trois longitudes observées aux temps 
| — m , t,t+n, respectivement, laJongitude 9 pour un temps 
quelconque t + z, intermédiaire entre ceux-là , se calcule par- 
la formule 

= (»+*) „■ + (£i) 

171/1 m [m+n) n {m+n) 

Il en est de même des latitudes. 


5 


( 34 ) 


XXXIII. A l'inspection tlu tableau précédent , on voit que 
le mouvement de la comète en longitude est très-petit'; cir- 
constance peu favorable pour l’application de notre méthode, 
et sur-tout pour la détermination du coëificient D. Cependant 
on peut observer que lorsque le mouvement en longitude est 
aussi lent que dans cet exemple , les observations ne sont pas 
sujettes à de si grandes erreurs , parce qu’ulors la comète est 
comparée pendant plusieurs jours à la même étoile, ce qui 
donne des différences assez exactes. 

Pour procédera l’application de notre méthode, on com- 
mencera par calculer les cocfficiens C et D au moyen des for- 
mules ( j 5) et (a 5 ), comme il suit : , > 




7-6ii 5546 7.58ai5gg 



jYemb. — o.oo 4 o 834 
— o.oo 38 aog 


+ 0.0077033 
— 0.0079093 


— 0.0079093 


D = — 0.0002060 


Lieu du © $37° 57' 4 " 
180 


A 5 7 * 5 7 ' 4 " 


'C . ** C * 


A — a°— uo° 4 a' 19' 

sin ....' 9.9710027 

tang b '.... g. 7801653 


• i 


A — d = j 1 1“ i4' 44' 

dm 9.9694327 

tang b'.... o.i5g3g2g 


9.7511680 
Piomb... + 0.563856 


o.ia 88*56 


— 1. 3453 ao 
+ 0.563856 


C = — 0.781464 
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(55) 

On remarqnern que si on eût augmenté a d’une minute , le 
coefficient C an roi t subi très-peu de changement ; mais le 
coefficient 1) se' serait réduit r — 0.0000237, valeur neuf fois 
moindre que l’autre. Une légère augmentation de plus dans la 
longitude a, au roi t fait passer D du négatif au positif; d’où 
l’on voit combien il est à craindre que la valeur de D ne 
soit pas donnée avec assez de précision par les observations. 
Dans des cas semblables , il conviendra d’éviter tout emploi 
de ce coefficient. Continuons cependant d’appliquer nos for- 
mules , comme si les observations données étoient exemptes 
d’erreur. '• 

Puisque C et D sont de même signe, il s’entait qn’on ar< R. 
c’est-à-dire qu’an 19 novembre, époque delà seconde obser- 
vation , la comète éloit moins éloignée qne la terre da soleil. 
Ensuite pour déterminer r ainsi que p, distance de la comète 
à la terre, on résoudra les équations (29) ; mais d’abord il faut 
calculer les coëfficiens A et 2 lt cos c, d’après les formules (28), 
où l’on fera fl = 5 j ', ayant soin d’ajouter au logarithme de fl 
le log. constant 8.2355821. 

A — a ss 111*5' 38" * 

COS. 9.556179 

COS h 9. 88898a 

COS C. 9.445161 

1R.1 o.2g5456 

(2R COS c)... 9.740617 

XXXIV. Cela posé , en substituant la valeur de R dans ks 
termes tout constans, les équations (29) deviennent 

p = A i.o3ga53^ 

1* ~ 0-974658 — p(aRcos c) + p*. 


C- 9**9 9 9°9 

fi. 096910 

* (6.471164 

R.... 9.994426 

1 : &..... 3.686i33 
1 : cos b.. o. 1 11018 

A i.2 5î56o 


( 5G ) 

Or sachant déjà que r est compris entre R et Rwi c, on trou- 
vera aisément par quelques fausses positions: 
log. r = g-ggo'tôa * 
log. f = 9.71S736. 

Les vraies valeurs de ces quantités calculées d’après les élé— 
mens connus de l’orbite (vol. cité pag. 71), sont 
log. r = 9.990071 
l °gf — 9- 709296. 

D’où l’on voit que par ce premier calcul la distance r est déter- 
minée , à moins d’un 1000*“, et la distance f à moins d’un iao 4 “; 
résultat beaucoup plus exact qu’on ne devoit s’y attendre , à 
cause de l’incertitude attachée à la valeur deD, et aussi parce 
que l’erreur des approximations étant de l’ordre 0', la solu- 
tion ne semble devoir être exacte qu’à environ un lèo 4 ™* 
lorsque 0 est de 5% qui sont le ia'“ des 58 jours pris pour unité 
(art. XXIII). . , . 

* Mais si on examine de plus près l’influence de cette seconde 
cause, on trouvera qu’elle n’est pas à beaucoup près aussi consi- 
dérable. En effet , la seconde des équations qu’on vient de 
résoudre , est une équation rigoureuse , qui ne peut être affectée 
que d’une très-legère erreur dans la valeur de cosC. Quant à 
la première équation où h est la seule quantité qui peut par- 
ticiper a la fois de l’erreur des observations et de celle de la 
méthode , il est aisé de voir , par les valeurs connues de r et 

, , . . . d p 1 dh dr 1 dh , 

de e , qu’on a a tres-peu-près — , et — = . — ; de 

r ’ , • p i5 h r lao h 

sorte qu’une erreur de par exemple, sur h, n’en donne 

qu’une de —7 sur p, et une de 7^ seulement sur r. Cela 

explique suffisamment l’exactitude que nous avons obtenue 

dans la solution précédente, sans avoir cependant une valeur 

de D fort approchée. 

XXXV. Venons maintenant aux équations indépendantes 
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de D et qui doivent donner des i^Mtats plus certains. Il faut 
d’abord calculer les coefficient P) Q, H, d’après les for- 
mules ( 3 i ). Voici le détail de ce calcul : 


tin (A— «*)... 9.9710017 
tin (A — a ' )... 9.9694327 

1 long b 0.2132187 

o.i5365/ji 

Nomb + 1 - 4 * 447*4 

— 1.7812457 

P = — 0.3567733 


•in (A— al... 9 9698779 «n (A— a’)... 9.9694327 
tin(A— a»)... 9.9710027 tin (A —a )... 9.9698779 
tangb' 9.7801653 tangb ° o.iSgSgag 

9.7110459 0.0987035 

— 0.526072g — 1.2551728 

— 1.2551728 1 : 1 

— 1.7812457 


■A — a , = no° 4*’ >9* 
A — a' = nt 14 44 

2A— a°— a'=22i 57 3 

tin 9.8250966 

langbt .... 9.9121887 

9.7372853 
— o.546 1 1 65 


• in [ A — a ).. 9.9698779 »in(A— a ).. 9.9698779 


®°*(A — a ').. 9.5591471 c*i(A — a 3 ).. 9.5484643 

tang l>° 0.1593929 long b' 9 7801653 

9.6884179. 9.2985075 

+ °>4 8 79979 + 0.1988417 

+ 0.1988417 

. + 0.C868 3g6 . 

— 0.5461 165 


Q= + 0.1407231 


«n(A— a').,. 9-9694327 1 «in(A— a*)... 9.9710037 «n(A— a)... 9.9698779 


tan g b 9.9131887 tang b 9.9121887 " tang b 3 ... 0,1593.929 

tangb» o. i5g3g2g long b '..... 9.7801*55 2 tang b '... 0.081 1953 

0.0410143 • 9.6633567 0.3104661 

4 1.0990420 4 0.4606Î48 - — 1 .62355 1 7 

4 10990420 4 ».55g6 7 68 

+ 1 .5596768 H = — 0.0638749 


Il faut ensuite substituer ces valeurs de P, Q, H, dans l’équa- 
tion ( 35 ). 
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( 38 ) 

( ir r^ p + * sin * q) 


Voici le calcul des coëfficiens dans lesquels la valeur de • sin ÿ 
a été déduite de log. R par la méthode de l’art. XII. 



9.892909 

P..../.. 

g.55a3g2 

Q q. 148366 

fl 

( : cos b... 

8.934552 
0,1 1 1018 

cos b 

«C " 

1 .o6i5ai 

I .o6l521 


8.933479 

*— T«* 

o.6i39t3 

o.oo55i3 

0.209887 
1 sinf.... 8.043643 

cos b 

i ,061 5a i 

R 

8.25353o 

*C “ 

H 

8.8o535o 

Nomb. — 

0.619426 

4.16319a 

4- 0-017928 

• - 

g.8G685i 

• 


— 4.163193 


’ 



— 4.t45a64 


(P'+Q* + H*) = 50.06817. 


On aura donc l’équation 

I = 0.512917 — ■ * p(4.i45a6) + p*(io.o34og) 
log .... 0.617053 1.001478, 

qu’il faut combiner avec l’équation déjà trouvée, 

r* = o . 97*658 — p ( 2 R cos c ) + f‘ • * 
log.... 9.740617, 

et on trouvera pour solution , 

log. r — 9-99°o5o 
log. p = 9 - 70 ^ 5 ^ ; 

valeurs qui ont toute l’exactitude qu’on peut désirer, puisqu’il 
n’y a que très-peu de différence entre ces valeurs et celles qu’on 
a données ci-dessus d’après le? vrais élémens de l’orbite. 
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( 3 g) 

XXXVI. Il faut maintenant déterminer ces élémens tels qu’ils 
résultent de nos trois observations : or par les calculs précé- 


dens on a 

co s b _ 

% 


... 0.209887 


TcT p 

0.613913 

• C Q 

g. 86685 i 
IC 3 

t 

9 ' 7 °g 58 a 

f 

... 9.70958* 

(•••• 9709582 

A coê b 

o. 3 Q 3 /,tj 5 

* co* b 

9 - 9 , 94®9 

/ 

p’... 9.576433 

* * • . . > 

—— p = 

#c 

2.106177 

Q — 

. »c V 

— 0.830748 


+ 

I H 

Ï* 

1^12775 
1 .09340a 

« sin 

Q— 

... 0.011057 
— 0.841805 

' ' . *1 


I)es deux nombres P' et Q' ainsi trouvés , il sera facile de 
déduire m! et ri par les formules ( 3a) , qui donnent 
m' — T? 1 sin A + Q' cos A 

ri = — P' cos A + Q' sin A . 

Enfin les trois autres quantités m, n, p se calculeront par les 
formules (27) , et on aura pour résultat les logarithmes sui- 
vans : 

m 9.881887 m 9-68 i 3 oi 

n g.7i54i8 ( — ri) o.ni83q 

p 9.510753 9.576433 

où l’on a indiqué que ri est la seule de ces six quantités qui 
soit négative. On conclura de-la 

*t .* •• ■ 

mri — m'n = — 1 . a34g55 o.ogj65i 

tpp — m'p — 0.131674 9.119499 

np ' — n'p = o.6i5t83 9.789004 

XXXVII. Puisque mri — m'n est négatif, le mouvement de 
la comète est rétrograde , et les équations ( 4i ) deviendront : 
m'n — mri — cos I 

— ( m p ' — m'p) — sin I cos 8 j/ïn 

— ( np ' — n'p) — smlsinSj/ân 


• (4o) 

Des deux dernières on tire 

mp'—m'p 9.119499 

np — n'p g. 789004 


long S o.66g5o5. 


Ainsi l’angle 8 est de 77* 55' 7', ou de 357* 55' 7*5 or il faut, 
par les équations précédentes , que tin S et cos S soient néga- 
tifs , puisque tin I ne peut l’être ni j/ân ; donc on a 
S = a5 f 55' 7", 

Connoissant S, on trouvera immédiatement I par la formula» 


qui donne 


tin S tang I = — 


( np' — n'p ) 
trin — mri ’ 


I = 26» 5g' 43* 5. 

Enfin on aura aussi la distance périhélie par la formule 

t ! 

— ^ m n " ii\ fz 11 * it • 1 . 

yan = — — — • ; dou Ion lire log. XI = g. 982^74, ou celle 

distance elle-même 


n = o.gGo44g. 

L’anomalie déduite de la formule cos* j 4 

4 «= is-ffée*. 


n 

= — . sera 
* * r ’ 


Avec celte anomalie, on trouve dans la table des comètes le 

1 

temps correspondant T = i o j • 9697 ; ensuite la formule <=n * T, 
donne /= io>.3a54. C’est l’intervalle de temps entre l’observa- 
tion du ig novembre et le passage de la comète au' périhélie. Pour 
savoir si la seconde époque précède la première ou en est pré- 
cédée , il faut calculer la quantité k , ou seulement en déter- 
miner le signe d’après la formule 

k = mm + nn + pp'. 

Or à l’inspection des nombres qui entrent dans cette valeur, 
on reconuoit bientôtque k est négatif, et qu’ainsi à l’époque du 
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( 4i ) « 

19 novembre, 8 k 39' 4i" temps moyen , la comète n’avoit pas 
encore atteint son périhélie. Donc en ajontant h cette date l’in- 
tervalle trouvé io* 3a5i ou 10' 7* 48' 35", on aura l’instant du 
passage au périhélie : ag novembre i6 k 18' 19", ou novembre 
39.679 4. 

XXXVIH. Il reste à trouver le lieu du périhélie , et d’abord 
il faut trouver la longitude héliocenlrique p de la comète , 

par la formule tangp = — qui donne p = 34° 16' 4a", ou 

771 

p = ai4° 16' 4a* ; la première est celle qu’il faut choisir, parce 
que sin ? et cos » doivent avoir les memes signes que n 
el m respectivement , lesquels sont positifs. Enfin la formule 

la ng,= t ° n S- ( *- S -l , donne 

CO S I 

9 34° >6' 4a" tang(f-S).... 9979379 

S — o 5 y 55 7 cos 1 9.949899 

p — S = i36 ai 35 tangr 0.039480 

t = i33« 3' a 5" 

et par la quantité • — 4 + S , on trouve ce .que les astronomes 
appellent le lieu du périhélie , 

i5°5i'46*. 

Rassemblant ces divers résultats , on aura les élémens de l’orbite 
de la seconde comète de 1781 , tels .qu’ils résultent des trois 
observations données ; nous les mettons ici en comparaison 
avec les élémens corrigés qu’on trouve dans les Mémoires de 
*781 , page 71. 


ÉLÉMENS APPROCHÉS. 

% - 

ÉLÉMENS conniGÉs. 
• 

Distance périhélie o . 960449 

Lieu dn nœud ascendant 77° 55 ' 7* 

Inclinaison ‘....*>6 5 g 44 * 

Lieu du périhélie .5 5 1 48 

Passage au périhélie: novembre agi i6 h 18' jg“ 

Sen» du mouvement Rètrof 

0.960995. 

77" aa' 55 " . 
a 7 1a 4 
16 3 7 

agi «a* 4** 46 * 

rradi. 


6 


esAk. 


( 42 ) 

XXXÏX. La différence de ce» deux systèmes d’élémcns est 
assez petite pour qu’on ait lieu de s’étonner qu’un premier essai 
de calcul , fondé sur trois observations faites dans l’intervalle 
de dix jours seulement, ait pu conduire à une connoissancc si 
approchée de la véritable orbite. Mais pour comparer encore 
mieux les deux systèmes , nous avons calculé dans chacun 
d’eux les lieux géocentriques de la comète aux époques des 
trois observations. Voici le résultat de ce calcul : 


ÉLÉMENS APPROCHÉS. 

jkpvqucs. 

Long, calculée. 

Long, obeervét. 

Dijpér. 

J.ahi 

calcul 

LalU- obêerv. 

Difftr. 


D 

M 

S 

u 

M 

S 

s 

D 

M 

S 

D 

M 

s 

S 

'4 

3o7 

■ 5 

46 

3 »7 

>4 

45 

+ 6 . 

55 

20 

55 

55 

'7 

9 

+ 226 

>9 

3oG 

6 i 

=7 

tofi 

5i 

a 6 

+ • 

3 9 

«4 

49 

3 9 

«4 

4» 

+ ' 


3o6 

4i 

59 

3o6 

4a 

20 

— ai 

3i 

3 

5oj3i 

4 

5a 

— fia 




• 

ÉLÉMENS 

C O R R I G*É a 





* 

Î07 

*7 

5 9 

307 

i4 

45 

+ '94 

55 

i5 

• 4 

55 

«7 

9 

— 1 15 

>9 

3o6 

53 

6 

3oG 

5i 

26 

+ 100 

39 

• 3 

53 

3 9 

>4 

48 

— 55 

*4 

3o6 

43 

s 8 

7o6 

4> 

ao 

+ 38 

3« 

6 

34 

3i 

4 

53 

+ 9 J 


Il résulte de celte comparaison , que nos élémms approchés 
représentent mieux les longitudes que les élémens corriges, et 
qu’ils n’ont du désavantage sur ceux-ci que dans la seule lati- 
tude du i4 novembre. Ce résultat inattendu prouve que notre 
méthode a toute l’exactitude nécessaire pour une première 
approximation ; maison devra toujours corriger, par des cal- 
culs faits sur des observations éloignées , les élémens donnés 
par des observations peu distantes entre elles ; car on ne réussit 
pas toujours à représenter par une même orbite parabolique 
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( 43 ) 

toutcs'les observations d’une comète; et en cas de différence 
sensible dans les résultats, on doit préférer l’orbite qui satis- 
fait le mieux aux observations éloignées. Enfin si la différence 
étoit plus considérable , il faudrait avoir recours à une orbite 
elliptique ou hyperbolique. 


Autre application à la comète de 1769. 

XL. Pour apprécier encore mieux l’exactitude de notre 
méthode, nous allons l’appliquer , non à des observations tou- 
jours susceptibles d’erreurs d’où peuvent résulter des compen- 
sations; mais à des lieux delaæomctedc 17(19, calculés d’après 
les éléraens connus. Les différences qu’on trouvera entre les élé- 
mens déduits de cgs observations fictives et les vrais élémens , 
seront des erreurs dues tout entières à la méthode, et servi- 
ront à fixer le degré d’approximation auquel elle peut atteindre. 


Epoques. 

Longitude. 

Latitude austr. 

Lieu du O 

Log. R. 

Septembre. 
8 à . 14 *- 
10 » ici. 

1 a ici. 

D M S 

D M S 

D M S 

0.0016648 

0.002424a 

o.ooai 838 

a“ 101 18 S 

a lia 5 i 3 i 
a 134 a6 47 

4 " aa 14 35 
b a 3 a 8 i 5 
4 ‘ a 3 48 36 

. 

1 G6 35 3 1 
168 3 a aa 
1 70 39 20 


Voici d’abord le calcul des cocfficiens D et C : 


a— a* -=ii* 33' a3* 

sin 9.3017509 

tangb'.. 9.64469^6 
8.g464465 

tfomb. d* o.o883g88 
+ 0.0831493 
d- 0.1705480 


d—a — 11° 35' 16' 

si/i g. 3039138 

tangb’.. 9. Gi 16907 
8-gi46o35 

d- o.o83i4ga 


0'— a°= 

sin g.5g44435 

tangb.. 9.6376970 
g.a3si4o6 

— 0.1706635 
d* o. 1706480 
T) — — 0.000 11 55 


( 44 ) 


A = 

348 ' 3 a' aa* 

r 


A — a' = 

247° i 4 ' i 4 " A — a 

’ = 924 ' 5 ' 35 ' 

• . . 

sin 

9.9647849 sin 

.... g. 84 a 5 oo 5 

— o. 4 o 68 g 33 
+o.a 8457 i 4 

tan g b'.... 

g. 6416956 long b°. 

.... 9.6116907 


9.6094805 

g.454i9ia 

C = — o.i 223 a 19 

Nomb... — 

o. 4 o 68 g 33 

+ 0.28457 14 

Puisque C et D sont de même sigde, il s’ensuit qu’on a r<R, 

conirae dans l’exemple déjà apporté j ensuite 

Eaisant 8 — a’ , 

on aura 


• 

fo. 3 oio 3 oo 

A — a 

r= 235 ' 4 o' 5 l* 



(6.471164a 

cos..... 


c "... 

9.0875041 

cos b.. 

9.9694938 

R 

o.ooa 4 a 4 a 

COS c... 


1 : D 

3 .g 374 i 8 o 

aR.... 


1 : cos b.. 

0.0376062 

( a R cos ( 

?)... 0.0170749 

h 

9.8378^67 


XLI. Au moyen de ces valeurs , on aura à résoudre les 
équations • 

’p = h Qy — 0.983394) 

1.011226 — p(2Rcos c) + p*} 

et comme on sait que r est compris entre R et.Rsine, on 
trouvera aisément la solution 

log. r — 9.946165 
log. p = 9.506414. 

Les vraies valeurs de ces logarithmes , calculés par les élé- 
mem exacts, sont 9.945573 et 9.513676 , d’où l’on voit que 
l’erreur est + 5 qa sur le premier logarithme, et — 7262 sur 
le secoud , ce qui répond à une erreur d’un 740“" sur le pre- 

• * 
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mier nombre , et d’an 60"** sur Je second. Si ces premières 
déterminations ne sont pas plus exactes , c’est que 3 h cote est 
plus grand que R 4 , mais en diffère peu ; de sorte que, suivant 
i’art. XV, les équations à résoudre sont très-près de la limite, 
passé laquelle elles seroient susceptibles de deux solutions. 
Alors les deux courbes qui se coupent sont très-près de se 
toucher , et leur intersection se détermine moins exactement. 
Pour faire cadrer entièrement la solution de nos deux équa- 
tions aVec les vraies valeurs de r et de p , il suffi roit de prendre 
n5o au lieu de nâ5 pour la valeur de D, changement qui 
répond à un dixième de seconde sur a° ou sur d. 

XLII. Il faut passer maintenant au calcul des coëlHciens 
P, Q, H, d’après les formules (3i) : 


tin (A— O*)... 9.964784g 
«in (A — a )... 9.8 iî5oo5 
a tang b 9.9387270 

9.7460124 

Nom b •+■ 0.5572017 

• — O.5710873 

P= — o.oi38856 

«n(A — a ).. 9.9169325 
co»(A— a').. 9*856a5i8 
tang b ° 9.6116907 

9.384875o 

Ncmb ... — 0.2425912 


(A— a°).. 9.9647849 
(A— a )... 9.9169325 
g t>‘ 9.6446956 

9.52641 3o 

— o. 336 o 57 i 

— o.a35o3oa 

— 0.5710873 

tin[A—a ).. 9.9169325 
co«(A — a 0 ).. 9.5876 1 74 

tang b ' g. 6446956 

9.1492455 

— 0.1410086 

— 0,2425g 12 

— 0.3835998 


tin (A— a ')... g.84a5oo5 
tin (A — a )... 9.9169325 
tang b ° 9.6116907 

9.3711237 

— o.a35o3oa 


A — a°=ra47° >4' '4' 
A — a ' =» 224 5 35 

2 A- O 0 — a'=lll J 9 4g 

«■'» 9.9691824 

tang b 9.6376970 

9.6068794 

q- 0.4044635 
— o.3835gg8 


Q =î + 0.0208637 


( 46 ) 

»i*(A — a)... 9.9169I2S ii/i(A — a')... g.84a5oo5 n'n(A — a 6 )... 9.9647849* 


tang b ° 9.6116907 tan g b g. 6876970 tan g b 9.6376970" 

2 tang b ' 9.9467256 tang 4 ». ... g.61 16907 tang b '... t .. 9.644*3986 

9-4743488 9.0918882 9.2471775 

JSomb -f- 0.2980910 — 0.1233629 — 0.1766760 

0.3002389 M O.I766760 

H = — 0.0021479 — 0.8002389 


Si l'on vent voir quel est le résultat qu’on lireroit des équa- 
tions (33) et (34), il faut multiplier les valeurs de P, Q, if , 
R9 

par -jj- dont le logarithme est 3 . 4 7 G 45 .fi; d’ailleurs par l’ai> 
ticle XII on a log. ( t sin V ) = 8 . ig6y643 , çe qui donne : 


P.... 

... 8.i4a5647 

Q 8.3ig3gi4 

H 7.33aoi4i 


2.4764543 

# a. 4764513 

3.4 7 64545 

ïV 

0.6190190 
... 9-99 7 5i46 

0.7958457 
isint... 8.ig6g643 

y.8o84684 

R 



0.6 16 5336 

8.9928100 

♦ 

AomA. . . 

— 4.135553 

+ 0.098358 




— 4.i35553 



t 

— 4-o3 7 ig5 = 

- coëf. 3e u. 


O 11 trouve de mémo , pour le coefficient de a.*, 

~(P*+Q*+H*)= i4. iga485. 

Substituant ces valeurs dans les équations (33) et (34), et 
faisant toujours - = o, on aura h résoudre les équations 
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I ss o. 494474 — »(a.oi85g8) + **(7.096243) 
log.... o.3o5o5o o.85ioa8 

^ =o.9833g4 + ». 

Or connoissant déjà les limites de r, on trouvera facilement la 
solution 

hg. r = g.g45435 
log. » = 9.676148 ; 

R 9 c 9 

et de la valeur de », on déduira celle de c = » . ~ — , 

D 2 cos b 

qui donne log. p = g.5i3ig5. 

Les vraies valeurs de log. r et log. p étant 9.g45573et9.5i3676, 
on voit que l’erreur est — i38 sur le premier logarithme , et 
— 48* sur le second, ce qui ne fuit qu’un 3oüo im * d’erreur 
sur r, et un 900*“* sur p . 


XLIII. Venant ensuite à l’équation (35)-, on pourra déduire 
ses cocfficiens de ceux de l’équation précédente , au moyen de 
2 D co s b 


la valeur » =?•- 


R9’C 


-, dont on a déjà calculé le coefficient 


pour passer delà valeur de » à celle de p. On aura donc de cette 
manière les équations suivantes en r et p : 

- = o.4 9 4434-r(C') + f *(C') 

r 

log.... o.468oo3 1.176934 

r* = 1.011226 — p(aRcosf) + f‘ # 
log.... 0.017075, 

dont la résolution donne 


log. r = 9.945619 
log. p — g.5i3ioa. 

ôn voit que l’erreur est encore diminuée sur Ie*premier loga- 
rithme , mais qu’elle est augmentée à-pcu-près d’autant sur le 
second. 

Comme le résultat de re^dernières équations est celui qui 


> 
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mérite le plus de confiance, lorsqu’on n’a encore aucune con- 
noissance des élémens de l’oçbite , nous continuerons , d’après 
ce résultat , le calcul nécessaire pour obtenir une première 
approximation vers ces élémens. 

XLIV. Il suffit le plus souvent de connoître la distance 
périhélie , parce qu’avec cette distance et le rayon vecteur r, 
on détermine facilement l’anomalie vraie au temps de la seconde 
observation , et l’instant du passage de la comète par le péri- 
hélie. Or ces premiers élémens étant connus d’une manière 
approchée , on a les moyens de les rectifier et de déterminer les 
autres élémens de l’orbite , en comparant les observations les 
plus exactes laites à de grands intervalles de temps. 

Si donc on veut se borner à la recherche de la distance péri- 
hélie, il faut faire usage des formules (37) et (4o); savoir, 
n = r — j *■ , et 

k = p cos b sin (A — 'a).P'+(R — p cos c; Q’ + p s in b. H' ; (56) 

formule où Von a fait pour abréger, 


p/_ Lf£ilp_( l — ï*‘) 


Q' = 


9 C 

p cos b 
~ ~6C~ 


tt 

Q — * sin 'f. 


Or dans cet exemple on trouve 

p cos b sin (A— a).... g.3ga5a83 
4 F « 7. 8g 836 3o 




<«7) 


7.3908913.... +o.ooig538 


pt — p cos c 9.9127514 

Q' 0.1754585 

• . " 

0.0983099.... — 1.3537470 
csinb.H' 8.3967865.... +o.oig8o55 

• -—1 .a3ig88 = k. 


» 
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Puisque l est négatif, il s’ensuit qu’au 10 septembre, époque 
de l’observation moyenne , la comète n’avoil pas encore alteiut 
son périhélie ; on en déduit 

n =-r — ji* = o.ia34o8. 

XLV. Pour avoir à la fois tous les éiémens de l’orbite , tels 
qu’ils peuvent résulter de nos trois observations , il faut cal- 
culer les valeurs des six coëfliciens m , n , p , m\ ri, p', d’après 
les formules 


m — R cos A 4- p cos b cos a 
n = R sin A + p cos b sin a 

p — p sin b 



P' sin A 4- Q'ros A 
— P 'cos A + Q’si/J A 
H' 


et on trouvera leurs logarithmes comme il suit , en observant 
que m est la seule de ccs quantités qui soit négative : 

ni g. gSpaaSg ( — m).... o.i66a465 

n 8. 8 7 88g56 ri g. 48 48 1^5 

p. ...... 9.1133928 p'.. g.i834g37 


13e— là on déduit les logarithmes des trois quantités suivantes , 
dont les deux premières sont positives et la troisième négative 

{mri — m'n) 9.5756968 

(rap' — m'p ) 9 . 5093008 

, — (np 1 — rip ) 8.4485i)?>7 


XL VI. De ce que {mri — m'n) est positive, on conclut 
d’abord que le mouvement de la comète est direct ; il faut donc 
prendre le signe supérieur dans les équations ( 4 1 ) , et on aura 

mri — m'n = cos I j/an. 
mp — m'p *ss sin I cos S t/’ïn 
np' — rip ' — sin I sin S [/~ïT\ 

A l’inspection des deux dernières , on voit que sin S doit être 

4 * *’ ‘ '* * ' * 
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négatif et cos S positif; de sorte que l’angle S calculé par la 
formule 

8 . 4i85g37 

tang S= n P~ n P. - 9.5092008 

mp — m'p 

tang S 8.9393929 

doit être compris entre 270° et 360” ; et comme la tangente, 
qui est négative, appartient également aux deux angles 175° 1' 
45". 2 et 355" 1' 45". 2, on doit s’arrêter nu dernier, 

S = 355" 1' 45". a. 

C’est le lieu du nœud descendant, suivant la règle de l’art. XXV, 
puisque le mouvement est direct et la latitude australe. 
L’inclinaison I se calculera par la formule 

mp — ni p 

tang I = ; j -!- — , 

[mn — m n) cos S 

qui donne I = 4o° 44 ' i4 '7. 

Enfin la distance périhélie étant calculée par l’équation 

( mn'— mn )* 

n = — —, 

2 cos' I 

on aura log. n = 9.0913596 , ou n = o. 1234126, valenr qui 
s’accorde suffisamment avec celle que nous avons déduite de la 
valeur de k. 

La distance n étant connue , on «ara l’anomal ie 4 par la for- 
mule cos' ï 4 = — -* d’où l’on tire 4 = i 36° 4' 3o".4. 
r 

Cette anomalie répond, dans là table des comètes, an temps 

T = 621'. 749, et ce nombre étant multiplié par n", on aura 
t = aS'.g56a, qui ajouté au temps de l’observation moyenne 
sept. io.5833, donnera sept. 37-53g3 ou octobre 7-53g3 pour 
l’instant du passage de la comète par le périhélie. 

La longitude hélioccntrique t est donnée par la formule 
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long P = — , où l’on doit observer de rendre le signe de sin p 

771 

conforme à celui de n , ou celui de cos p à celui de m. On aura 
ainsi p — 4° 58' 26" et p — S = 9 0 56' 4o".8. Ensuite la dis- 
tance <r de la comète au nœud dont la longitude est S , se calcule 
par la formule 


tan g r = 


tang (y — S) 
co s I 


» 


qui donne r — i3* 1' 5a". 8. Enfin la somme S + <r+4, puisque 
la comète marche vers' son périhélie (art. XXVI), donnera 
pour le lieu du périhélie, i44° 8' 8". 4. 

"Voici le résultat de tous ces calculs comparé aux vrais élé- 
mens de l’orbite : 

I. 


ÉLÉMENS APPROCHÉS. 

ÉLÉMENS VRAIS. 

Log. dût. périhélie.... 

9.091360 

•i 

9.090847 

Passage au périhélie Octobre 

7 . 53 9 3 


Octobre ... 7 . 53 io 

Lieu du nœud ascendant 

i 7 5 # 

45" 

i 7 S* 3- 4o" 

Inclinaison. 

40 44 

i5 

40 47 5 fî 

Lieu du périhélie. «... 

.44 8 

8 

144 n 3 a 

Sons du mouvement.. 


. Direct. 




OBSERVATION GÉNÉRALE. 

. ...... .. è * * il... .fUii 


XLVII. Les deux exemples que nous avons choisis pour 
l’application de notre méthode, sont aussi différens qu’il est 
possible. Dans le premier, le mouvement en longitude est très- 
lent , et le mouvement en latitude très-prompt. Dans le second, 
au contraire, Je mouvement en longitude est très-prompt, et 
le mouvement en latitude très-lent. Dans l’un , l’intervalle des 
observations est de 10 jours ; dans l’autre, il n’est que de 4 ; 
et cependant le succès de la méthode a été à-peu-près le même 
dans l’un et dans l’autre. On a obtenu les élémens de l’orbile 
très peu différens de ceux que donne la combinaison d’un plus 
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grauri nombre d’observations ou d’observations faites à de plus 
grands intervalle» de temps. 

Il y a lien de croire que le succès sera le même dans les 
diverses applications qu’on voudra faire de cette méthode , sauf 
peut-être quelques exceptions qui pourront dépendre ou de 
l’inexactitude des observations ou de quelques circonstances 
particulières peu favorables, telles qu’une très-petite inclinai- 
son de l’orbite, qui ne peruiettroit pas d’employer avec sûreté 
les latitudes ob>ervées. Le cas le plus défavorable de tous seroit 
celui où la vraie orbite différerait très-sen>iblcment d’une para- 
bole , et où en môme temps la distance de la comète au soleil , 
lors de la seconde observation, différerait très-peu de celte de 
la terre nu soleil ; car alors l’emploi des équations (39) pourra 
être défectueux à cause de l’incertitude sur la valeur de D, çt 
on 11e pourrait pas non plus faire usage de l'équation ( 35 ) y’ 

--f 

puisque - 11e serait pas assez petit pour être négligé. 

Au reste, nous avons réduit les formules autant qu’il a été 
possible pour la comtiiodilé des calculateurs : nous avons donné 
des règles sûres pour distinguer , dans le résultat des formules, 
le nœud ascendant du nœud descendant, les époques avant le 
périhélie des époques après le périhélie , et enfin pour prendre 
chaque chose avec le signe qui lui convient, sans avoir recours 
aux ligures qui sont cependant utiles pour diriger le calcul. 

- Mais on ne saura p >reco romancier aux calculateurs «l’ap-, 
porter le plus grand soin à bien déterminer les signes des différens 
termes qui entrent clans Je» formules, et principalement dans 
les valeurs des coeffieieus C, D, P,Q, II, d’après les règles 
connucsdes sinus, < usions et tangentes considérés dans le» divers 
quarts de la circonférence. L’attention , à cet égard, est d’aulaiit 
plus necessaire, que rien n’y peut suppléer, et que la moindre 
négligence peut faire augmenter considérablement un travail 
déjà long et fastidieux, ou donnera croire mal- à-propos que la 
méthode qu’on a suivie est défectueuse. 
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SECONDE PARTIE 


Méthode pour corriger les élémcns de l'orbite connus 
par une première approximation. 


XL VIII. IHous prendrons pour exemple la comète de 1769 
qui a été observée pendant un long espace de temps , avant et 
après son passage par le périhélie, cl nous choisirons les trois 
observations suivantes tirées de la Cométographie de Pingré, 
tome n, page 368 . 







Temps moyen 
à Paris. 

Longitude' 
de la Comète. 

Sa Latitude. 

Lieu du Q 

Log. BS 

Aoét. 14. 5235a 
Sept. . 1.5.by3y8 

Déc... 2.21413 

D M S 

D M S 

D M S 

0.00 5 2 440 

0 . 00 1 8 1 00 

<) 

3 <j 58 16 

1 4 o 3 g 17 
*76 4 | 20 

3 17 3 o A 

22 43 3 « A 

a 3 33 26 B 

l/|2 21 26 

173 3 i 3 o 
260 54 1 2 


Supposons que par les premiers essais on a trouvé le loga- 
rithme de la distance périhélie == 9.0900000, et l’instant du 
passage au périhélie le 7 octobre à 12 heures environ; nous 
prendrons pour ces deux élcmcns corrigés , 


• r >■■! 


log. n = g. 0900000 + « (10000) 

Passage au périhélie octobre 7.50 + t ( o . 35 ) 

La correction •» ( toooo) est exprimée en unités décimales dn 
septième ordre ; et la correction r (o.a 5 ) représente une frac- 
tion de jour qui iroit à { de jour ou 6 k , si on avoil t = j. Ces 
nombres sont ainsi choisis , parce qu’on suppose que lu véri- 
table époque du passage au périhélie, doit être entre 7.25 et 
7.75; et que le vrai log. de la distance périhélie est pareille- 
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ment compris entre 9,089 et 9,091 , ou du moins qu’il s’écarte 
peu de ces limites. 

XLIX. OU posé à l’aide de ces deux éléinem, on doit cal- 
culer les anomalies et les rayons vecteurs qui ont lieu à l’ins- 
tant de chaque observation. Voici le détail de ce calcul appli- 
qué à la seconde observation : on y verra les précautions à 
prendre pour tenir compte des corrections indéterminées. 

Passage au périhélie... ocl 1 "' 7.60, ou sept 1 " 37.50 + 1(0. a 5 ) 
Temps de l’observation sept 1 ” 15.69398 

Différence .. t == ai’, 80602 +1(0’. a 5 ) 

I<e logarithme de ai . 8060a est ï.33£576i ; pour nvoir la partie 
additionnelle due à la correction 1 (o. a 5 ), on pourroit prendre 
la différence logarithmique de 81806 à 81807, laquelle est 199, 
et multiplier cette différence par a 5 o , ce qui donneroit 49750 ; 
d’où résulte le logarithme entier de t = 1 .3385764+1 (49750). 
Mais pour plus d’exactitude, il faut prendre la différence qui 
répond à une unité plus élevée, moitié en dessus , moitié en 
dessous du nombre constant; par exemple, celle qui a lieu de 
21756 à ui 856 . Cette différence est 19916 : on la trouve facile- 
ment, par la disposition des tables, en remontant et descendant 
de cinq degrés dans la même colonne au-dessus et au-dessous 
de 91806. En vertu de cette différence, l’excès logarithmique 
qui réptyid à t(o.io) sera 1(19916); donc l’excès qui répond 
à t (o.s 5 ) sera 1(49790), peu différent de celui qu’on avoit 
trouvé immédiatement. On aura ainsi, 

log - t — 1 . 3385764 + 1 ( 49790 ). 

D’ailleurs on a 

ï log- II = 8 . 635 oooo + «• ( 1 5 ooo ). 

Retranchant le second du premier , il restera 

3.7035764 + 1(49790 ) — « (i 5 ooo). . - 
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C’est le log. du temps T employé à parcourir une égale ano- 
malie 4 dans la parabole dont la distance périhélie est i. 

Le logarithme 3.703576-'» répond au nombre 5o5.33a. Pour 
avoir la partie additionnelle, j’observe que l’unité de diffé- 
rence entre les nombres 5o4.83 et 5o5. 83, produit entre leurs 
logarithmes la différence 8699;. donc à proportion les diffé- 
rences logarithmiques 49790 et i5ooo , donneront entre les 
nombres les différences 5.796 et 1 -748. De sorte qu’on, aura 
T = 5o5>.33# 4- t(5> 796) — *(i* 746). 

Dans la table du mouvement des comètes , on trouve pour 
5o4 jours, l’anomalie i3a” 19' 69', et pour 3 jours de plus, 
hne différence de i'iy'’ ou 4' 4^o. Donc l’anomalie qui répond 
à la valeur de T est 

4= i3a # 23' 56" 8 + t ( 1 s' 894 ) — «(3' 885), 
où l’on voit que les corrections sont exprimées en minutes et 
millièmes de minute, ce qui est plus commode pour le calcul 
que des expressions réduites toutes en secondes , attendu qu’il 
faut l’apporter les différences à une unité au moins égale à la 
minute. 

Il faut chercher maintenant le logarithme-cosinus de l’angle 

14 = 66* 11' 28*4 + t(6' 447) — *(i'94a). 

• 

Ce logarithme, pour la partie connue, est 9.6060431 : pour 
avoir l’autre partie , je trouve dans la table qu’une minute de 
plus dans l’angle fait diminuer le tog.-co>inu » de 3865, et 
qu’une minute de moins le fait augmenter de 286a. Le milieu 
de ces différences est 3863 j- pour 1' ; donc pour 6 ^47 la diffé- 
rence sera i846i , et pour i'9'*a, elle sera 556i, On aura 
donc 

log. coi J 4 — 9- 6 o 6 o 43 i — r (18461) + *(556i) 
son double.... g. 3120863 — r (36gaa) + * (1 113a) 
log. n g .0900000 4- *(10000) 

Diff. ou log. r = 9.8779138 •+• r (36g3a) — * (usa) 
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c’est le logarithme du rayon vecteur calcule d’après la for» 
mule ( 38 ). 

L. Pour aller plus loin , il est bon de diriger le calcul par 
une figure qui représente à-peu-près la position relative de la 
terre T, du soleil S, de la comète C, et de sa projection K sur 
le plan de l'écliptique ( fig. a ). 

Longitude du © 173* 3 i' 3 o" 

Long, de la comète. • . , i 4 o 3 g 17 

Elongation KTS 3 a 5 a i 3 

Il faut en conclure l’angle CTS entre la comète et le soleil 
(angle qui a cté appelé c dans la première partie ) ; on le trouve . 
par la formule cos CTS = cos KTS cos CTK , où CT K 
= b = aa* 43 ' 34 ", 

cos KTS.... 9. 9247384 
cos b 9.9649015 

cos CTS.... 9.8891 a<)9 CTS = 39” i 3 ' a 1*9.’ 

Dans le triangle CTS, connoissant les deux côtés ST = R, 
CS = r, et l’un des angles opposés CTS , on calculera d’abord 

l’anglcTCSpar la formule s in TCS — — sin CTS , qui donne 

R 0.0018100 

sin CTS.... 9.8009486 

o.iaao86a — r ( 36 gaa ) + « (naa) 

sin TCS. . . . 9 .ga 48448 — r ( 36 gaa ) 4- «(lias) 

De là résultent deux valeurs de T CS (parce que'l’angle CTS 
n’est pas obtus ) , et il faudra se décider entre l’une ou l’autre , 
d’après la connoissance approchée de l’orbite- Dans ce cas , on 
doit prendre la plus grande , savoir 

TCS = iaa” 44 ' 38"8 + t( 45 ' 47 o) — *(»' 33 a); i 
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et comme on a déjà CTS = 39' i 3 ' ai" g , le troisième angle 
du même triangle 

CST = 18" 1' 59*3 — t ( 45 ' 470) + *(i' 38 a). 

Cela posé, la latitude héliocentrique CSK= * ^ se calculera 
sinbain CST 


par la formule sin a = 


CTS 


sin b g . 586 g 544 

si» CST.... 9.4907547 — t (176446) + «( 5363 ) 

1 : ain CTS.... O.iggo 5 i 4 

sin a 9.2767605 — t ( 176446) +■ «( 5363 ) 

De-là résulte a = 10” 54 ' 7' a — &c. ; mais les corrections de 
l’angle sont inutiles à calculer, et il suffit d’avoir celles du log. 
cosinus. On peut trouver celles-ci, en observant d’après les 
tables , qu’une différence de 65600 sur le log. sinus , répond à 
une de a 434 sur le log. cosinus : on peut aussi multiplier les 
corrections de log. sin a par tang' a , afin d’en déduire celles de 
log cos a avec un signe contraire. Ce dernier moyen , qui est le 
plus simple , donnera 

log vos a = g. 9920903 4- r ( 6546 ) — «(198). 
Maintenant l’angle de commutation K ST se calculera par la 
cos CST 

formule cas K ST = . 

COS A 

cos CST.... 9.9781236 + r (1871 1) — »( 56 g) 
r : cos a 0.0079097 — r (6546 ) + « (198) 

cos K ST... 9.9860333 ■+• 7 (1 ai 65 ) — «(371) 

Commutât. KST = i 4 * 27' 11' 3 — - t( 37 ' 36 s)+« (i' i 3 g) 
Longit. de la terre... 353 3 i 3 o . 

Longit hélioc. t = f 58 ' 4 i* 3 — r ( 37 ' 36 ï)+« (i' i 3 g) 

U. Par ces calculs, on a déterminé l’anomalie 4 > la latitude 
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héliocentrique a , et la longitude héliocentrique p de la comète 
pour le moment de l’observation du 1 5 septembre. On fera de 
semblables calculs pour les époques -données du i * août et 
du i décembre , et on aura les résultats suivons , dans lesquels 
on a distingué par l’accent ° ceux du i4 août , et par l’accent ' 
ceux du u décembre : < , 

[ 4"= 146” i3' 5a' + r ( 3^53 ) -—•(a'575) 

i t août. ' s= 35 7 i5 54 .4 + 1 (io'4^a) — *(o'5a3) 

(/off. sin *•’ = 8. 5367204 + V ( i8g54 ) — • (g45) 

iSaepierabref 4 = i3a" aa' 56*.8 + t ( i a'8.)4) — (3'885 ) 

avant < f sa. 7 58 4* -3 — T(37'36a) + ■»( I ' 1 3t^ ) 

le périhélie, (kg. sin*. =s g. 3767605 — 1(176446) + «(5363 ) 

2 décembre j 4' — l4®* 37’ 1 8* — t ( 3 3oo ) — 'b ( a’54o ) 

après < f' — 2Q7 3o 20. 3 + r ( 4'78o ) + « (o'a4* ) 

le périhélie, [ kg , sirt K > _ ' 9 . 7698879 + t( 3 7 3o ) + * ( 188 ) 

LII. Concevons maintenant une sphère concentrique au 
soleil , dont ta surface soit rencontrée par le plan de l’écliptique 
et par celui de l’orbite de la comète. Four mieux saisir la situa- 
tion respective des différons points, on pourra supposer que la 
circonféreuce de l’écliptique est développée sur un plan (fig. 3), 
et forme la ligne droite QgQ dans laquelle T représente le 
premier point d’or/es , T® l’ordre des signes , n le nœud ascen- 
dant , et u le nœud descendant de la comète. A chaque point 
de l’écliptique, comme K, on peut élever une perpendicu- 
laire KC qui représente la latitude de la comète, en même 
temps que T K est sa longitude ; et la suite .des points C for- 
mera la ligne sinueuse uC u C' il qui représente l’intersection 
du plan de l’orbite avec la surface sphérique concentrique au 
soleil. Dans cette sorte de projection , les longitudes et les lati- 
tudes sont représentées par des grandeurs proportionnelles ; 
mais les arcs du mouvement héliocentrique sont altérés dans 
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leurs proportions ; ils sont en général augmentés , puisque 
Parcdecourbe 0C0 ne représente qu’une longueur égale à u* 1 - 
Quoi qu’il en soit , cette figure est très -propre à diriger le cal- 
cul qui reste à faire pour la détermination de l’orbite. 

La distance Tg , supplément de la longitude du noeud ascen- 
dant O, est d’environ 5*. Au 14 août, la comète étoit à un 
point C“ très -voisin du noéud descendant, puisque la dis- 
tance T K’, complément à 3Go* de la longitude p*, n’est que 
d’envirou 3*. Au i5 septembre , la comète marchant vers le 
périhélie n, s’est trouvée en C avec la longitude T K d’envi- 
ron 8’. Arrivée au périhélie rf le 7 octobre , elle s’est trouvée 
distante d’environ 3i* du noeud ascendant Q. Continuant sa 
route, elle a passé le nœud, sa latitude est devenue boreale, 
et enfin le 3 décembre, elle est parvenue au point C' où l’on 
a rK' = 6a° 3o' environ , valeur de 36o* — r'. Si l’orbite de la 
comète est exactement parabolique, il ne lui reste à parcourir 
du a décembre 1769 jusqu’à l’infini, que l’arc C't , jtq étant 
égal à nu. Si cette orbite est elliptiqup , elle continuera son 
mouvementjusqu’à saprochaine apparition dans la partie uC’C, 
mais elle ne décrira qu’un angle très-petit dans un temps con- 
sidérable. 

* * * * * 

LIII. Après avoir pris cette idée du mouvement bcliocen- 
trique de la comète , venons à la détermination tant des coiiffî- 
ciens r et » que des divers élémens de l’orbite. -, 

Suivant la méthode ordinaire, on considère le triangle splié- 
riqne formé P ar l ea deux points C*,C, et par le pôle de l’éclip- 
tique. Dans ce triangle, on connoît les côtés de l’angle au pôle, 
qui sont les complémens des latitudes K*C*, KC ; avec l’angle 
compris mesuré par l’arc !QK= 3 60* — *"+#; on peut donc 
déterminer le troisième côté C’C qui doit être égal à la diffé- 
rence des anomalies 4° — 4 j et de-là résulte l’équation de 
condition 

cos { 4° — 4) = cos (#—?“) coj a cos\°+sin a sin a*. (55J 
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On en obtiendra une semblable de la comparaison des deux 
lieux C et C' qui donne 

cos ( 4’ + 4 ) = cos (9 — e' ) cos a cos a'— sin a sût a' ; (56) 

et par le moyen de ces deux équations , on pourra déterminer 
t et ». 

Cette méthode paroit conduire assez directement au but ; 
mais elle a l’inconvénient de supposer que les trois points 
C”, C, C' sont exactement dans un même plan , et elle n’offre 
pas les moyens de corriger ou au moins de diminuer l’erreur 
qui auroit lieu si cette condition n’étoit pas remplie. 

On n’exprimeroit point la condition essentielle de l’unité 
du plan , en formant par la comparaison des points C et C' une 
troisième équation semblable aux équations (55) et (56). En 
effet, cette troisième équation seroit nécessairement une suite 
des deux autres; car il faudroit que les deux plans C'C’, C’C 
lissent en C’ un angle bien différent de 180 °, pour que l’arc CC' 
mené directement de C en C', fût sensiblement plus petit que 
la somme des arcs C° C f , C” C. Voici donc le moyen qu’on pour- 
roit employer. 

LIV. On reconnoit que trois points C‘, C, C’, sont dans le 
plan d’un même grand cercle , si les longitudes de ces points 
t, t', et leur» latitudes a", a, a' satisfont à l’équation 

sin ton g a '+$in (*'—») tang tf+sin (?*—»' ) long a = o , (5y) 

ou ce qui revient au même, si entre les anomalies 4 °, 4 , 4 ', 
et les latitudes a', a, a', on a l’équation 

«»( 4 — 4 ”) sût a' + sin ( 4 '— 4 ) *in a °+sin ( 4 °— 4 ' ) sin a =5 0 . (58) 

Si donc on substitue les valenrs données dans l’article LT, 
en observant de prendre a' et 4 ' négativement , on^bura une 
troisième équation , laquelle devra être combinée avec les 
équations (55) et (56), afin d’en tirer des valeurs de r 
et de», qui y satisfassent aussi bien qu’il est possible. Je dis 
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aussi bien qu’il est possible, parce qu’iiyant trois équations 
pour déterminer deux inconnues, on ne doit pas prétendre 
satisfaire exactement à ces trois équations , mais seulement 
chercher à rendre leurs premiers membres très-petits. 

LV. Connoissant r et •», toutes les quantités contenues dans 
le tableau de l’art. LI deviendront entièrement déterminées ; 
on calculera ensuite l’inclinaison I par la formule 
T sin {9 — S°) 

1 = r ;îü(r — 4 ) cotKcoa » » (%) 

où l’on aura soin de prendre 4* — 4 (qui représente la dis- 
tance CC°) de même signe que 9 — 9°, parce que l'inclinaison 
est toujours supposée moindre que 90”. Lorsque I sera déter- 
miné , on aura la longitude du nœud S par la formule 

sin ( r — S ) = tau g a cos I. ( 60 ) 

Cette équation , dans laquelle a doit toujours être regardé 
comme positif, donnera pour l’angle 9 — S deux valeurs com- 
prises entre o° et 180 0 ; on en tirera donc deux valeurs de S. 
L’équation sin (9 0 — S ) = tang a’ cos I , {jui a lieu semblable- 
ment, donnera de même deux valeurs pour S. On aura enfin la 
troisième équation sin( S — i) =* tang*.' cos I ( où l’on a changé 
Je signe du premier membre, afin de prendre positivement la 
latitude a', qui est de signe contraire aux deux antres a, a*) , 
d’où résulteront encore deux valeurs de S. Or, dans chaque 
système, il est nécessaire que la même valeur de S se retrouve, 
ou exactement ou au moins , à une petite différence près ; celte 
différence venant de ce qu’on n’a pas pu satisfaire exactement 
aux trois équations qui déterminent t et *. On prendra donc un 
milieu entre ces trois valeurs de S peu différentes entre elles , 
et on aura la longitude d’un nœud auquel on appliquera la 
règle de l’art. XXV, pour savoir si c’est le nœud descendant ou 
le nœud ascendant. 

Par les valeurs connues de l’anomalie 4 et de la distance 
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périhélie n, on trouvera aisément l'instant du passage an péri- 
hélie ; le signe positif ou négatif de s — p n , indiquera si le mou.- 
veruent est direct ou rétrograde ; enfin le lieu du périhélie se 
déterminera par les formules de l’art. XXVI. Ainsi la connois- 
sance de tous les élémens de l’orbite se déduit aisément de la 
résolution des trois équations (55), (66) et (67) on (68), 
qui donnent les coëfiiciens ret*. 

Nous ne faisons qu’indiquer oette méthode sans Rappliquer 
à notre exemple, parce que, quoique sa marche soit assez 
naturelle, cependant elle a l’inconvénient de ne pas réduire 
les erreurs a une même échelle. Deux erreurs égales dans deux 
des trois équations que l'on considère pourroient répondre à 
des erreurs très-inégales sur les longitudes ou sur les latitudes ; 
et c’est ce qu’il faut éviter. Nous allons donc exposer une autre 
méthode qui paroît devoir mieux remplir notre but. 

LVI. Supposons que par les formules précédentes ou pur 
une première connoissance approchée des élémens de l’orbite, 
011 ait trouvé I = 4o° 5 o', S — 355 ° 3 ' (c’est le lieu du nœud 
descendant), et aQod r = 16" 48'. 

Appelons pour abréger * la quantité f — S , ou le côté oK 
du triangle sphérique rectangle tjKC , on aura dans ce triangle 

les deux équations 

ain \ = ain tain r (61 ) 

tong a = cos I long o. (6a) 

Nous prendrons pour valeurs corrigées de I et de r les deux 
expressions 

I = 4o° 5o'.+ z (6000) 

<r = 16° 48' + « (ô'ooo); 

- et par leur substitution dans l’équation (61) . on aura 

ain 1 9-8i5i854 + z(73oo) 

ain r 9.4609466 + 11(20920) 

. tin a 9. 97G4310 +.4(7300) + «(30920). 
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Comparant ce iog. sinus à celui de l’art. Ll, on a la différente 
— 3 ag 5 + * (7800) + u ( 30920) + r (176446 )— «( 5363 ) , 
d’où résulte sur la latitude a l’erreur 

E' = — o' 5 oo + z (i'i i3) + ij ( 3 ' 189) + 1 (36 900) — « (o' 8 i 8 ). 
Calculant de même les deux formules sin k° = sin I tin <r , 
si» a' = sj» I tin <r‘ , d’après les valeurs r* =a o — 4“ + 4, 
«' ~ 36 o”— » 4'— .4 — <r, on trouvera deux valeurs de log. tin a", 
et log. tin a', lesquelles comparées aux valeurs semblables déjà 
trouvées (art. LI ) , donneront sur les angles a”, a' les erreurs 
K“ = o'ia 5 +• z(o'igt) + »( 3 'afi 4 ) + r (5*657 )-~*(o'» 3 o) 
E"* = — o'ia 3 + *( 4 'ao 8 ) — «( >'761 ) — T( 5 ‘ 5 ag) + «(a'i 54 ). 
Maintenant le calcul de l’équation tang k=cos I tang o, donnera 

cot 1 9.8788748 — z ( 546 o) 

tango 9.4798887 + jj(aa 83 o) 

tang 9.3587635 — z ( 5460) + u (aa 83 o) 

x=3i“5a' 3 " 4 — z (o'g 38 ) + u(3'9a3) 
d’ailleurs s = 7 ° 58 ' 4 i -3 — T( 37 ' 36 a) + «(i'i 3 g) 

donc-Son 1*0 — 4 ° 53 ' 22''i-z(o'g38)+«(3'g23}+T(37'362)-«( 1 ' 1 3 <j) 

Les deux équations tang k°=cos! tango* et tang z'op cos I tango' 
donneront semblablement ; . ,.v> 

t 0 = 4 ° 58 ' 7"5 — z (o' 1 68) + «(3787)— t ( 3 ' 33 j ) — «(o' 46 g) 
r u = 5 ’ 4 ' 38 "5 + z(rg 6 i)+jj( 5 ' 78 g)+ t(6 J 38 j — «(7'68o). 
La seconde valeur de Tg étant retranchée surcessiveraent de 
la première et de la troisième , on aura deux nouvelles diffé- 
rences ou erreurs , 

E" = — V767 — z (0*770) -f u (o' 1 36 ) + r (4o'fig 3 ) — « (0*670) 

E’ = 6 ' 35 o *t- z (a' 139) + « (a'002) + r (g'ôSg ) — « (7'a 11). 

LVII. Toutes les conditions du problème étant ainsi expri- 
mées, si on vouloit anéantir à la fois les cinq erreurs E , E', 
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E", E", E', on aurait une équation de plus que d’inconnues, 
ainsi que le comporte la nature de la question. Mais comme ce 
ne serait que par un hasard singulier que les valeurs des incon- 
nues z, u, t, *, tirées de quatre de ces éqnalions, satisferaient 
à la cinquième , et que même dans ce cas , la méthode que 
nous de votls suivre conduirait au vrai résultat, nous tâcherons 
seulement de diminuer les errears de manière à avoir la para- 
bole qui satisfait le plus exactement qu’il est possible aux trois 
observations données. 

Nous commencerons donc par déterminer r et « par la con- 
dition que les erreurs E" et E" soient nulles ; ce qui donnera 
r = o. t344 4- x (o.oa43) 4- a(o.ooi3) 

* = i .o6o5 + * (0.3^78) + 0(0-2793). 
Substituant ces valeurs dans les trois quantités E' , E", E”', on 
aura 

E' = a'.s48 + *(i'- 499) + oCa'.ggG) 

E' a= o'-oo6 + z (o'.oRg) + «(3 .040) 

E'" ~ r.-fiS ■+■ s (4\ 780) — u( 1'. 167). 

On voit déjà qu’en faisant a et u égaux à zéro , toutes les 
erreurs se réduiront à deux; l’une E d’environ a' j, l’autre E''' 
de i'4i; erreurs très-tolérablcs dans la théorie des comètes. 
Mais il est possible de les atténuer encore en cherchant le 
minimum de la somme des quarrés des quantités E', E', E"'. 

LVIIJ. Pour cela , il faut multiplier tous les ternies de E' 
par 1*499 coefficient de x, tous ceux de E" par o.o5q , tous 
ceux de E'" par 4.780, et ajouter les trois produits. Opérant 
ensuite de la même manière par rapport aux coëfficiens de u, 
on formera les deux équations , 

o = 30.148 + *(35.099) — 0(0.908) 
o =ss 5.098 — z ( 0.908 ) + u (1 g.58o) ; 
d’où l’on lire 

z — — o.4(44 
u = — 0.9796. > 


O 
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Ati moyen de ces valeurs, les trois erreurs' E', E J , E" , 
deviennent 

E' = o' 789 , E” = — - o' 868 , E"=— o'a36j 

de sorte que la plus grande ne va pas à une minute. 

Substituant ces valeurs do-i et u dans celles de t et de •», 
on aura * > , . . 

T = o . 1 2,(0 

« =» o.8466; 

d’où résulte le log. de la distance périhélie — 0.0908466 , et 
l’époque du passage au périhélie : octobre 7,53 10. 

Des valeurs trouvées pour z et u , on déduit l’inclinaison 
I = 4o’ 47' 55" 7, et la distance <r ou uC = 16" 46' 36“i. 
Ajoutant à eC l’anomalie 4 = i 32° sa' 56*8 + T(ia'8g4) — 
* (3' 885) = i3a° si' i5"4 , et retranchant la somme de 180“, on 
aura la distance du périhélie au noeud no = 3o° 5a' 8' 5. Enfin 
par l’une des expressions de Tg on trouve Tu = 4 0 56' tg"6; ce 
qui donne la longitude du noeud ascendant 175° 3' 4o"4. 

Si de celte longitude on retranche la distance n n du péri- 
hélie au nœud , il restera ce que les astronomes appellent le 
lieu du périhélie = i44° u' Si* 9- 

LIX. On peut donc établir ainsi les élémens de l’orbite de 
la comète de 176g, en tant qu’ils résultent des trois observa- 
tions données : 

Log. dist périhélie * g. 0908466 

Lieu du nœud ascendant 175° 3' 40* 

Inclinaison de l’orbite 4o 47 56 

Lieu du périhélie i44 11 3a 

Passage au périhélie octobre 7 . 53io 

Sens du mouvement direct. 

Si d’après ces élémens, on calcule la longitude et la latitude 
géocentriques de la comète à l’instant de chaque observation, 
on trouvera les résultats suivans : 


9 
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Temps moyen 

Longitude 

Longitude 

Viff. 

/attitude 

Latitude 

Diff. 

à Paris. 

calculée. 

obiervèe. 

calculée. 

observée'. 

Aoùl. 1 ;.5î3 3a 

3g°58' a a' 

3 9 ° 58' 16 " 

+ 6* 

3‘i5'4i 

3 * 1 y » 3 ' 

- 89 " 

Sept. ID.by3q8 

140 îg 14 

i4o 3g 17 - 

— 3" 

aa 44 5t 

aa 43 34 

4-77 

Dér. a.ai4t3 

>76 4i aa 

376 4 > >0 

+*" 

a3 33 16 

i3 33 a5 

— 9 


où l’on voit que les erreurs sont insensibles sur la longitude, 
et ne vont qu’à l'j sur la latitude. 

Les élémens trouvés par Pingré , d’après les mômes observa- 
tions (Cométographie, tome n , page 58i), représentent bien 
les observations extrêmes du 14 août et du 2 décembre; mais 
sur l’observation du i5 septembre, ils donnent une erreur de 
+ a' 12* sur la longitude, et une de + 2' 37" sur la latitude. 

LX. Au reste , si d’après les élémens que nous avons trouvés, 
ou d’après ceux de Pingré, on calculoit les lieux de la comète 
pour d’autres éppques où elle a été observée, on trouveroit 
des erreurs de 12 ou 1 5 minutes et môme plus, tant sur la 
longitude que sur la latitude. C’est par cette raison que Pingré 
lui- même trouve différentes orbites eu calculant différentes 
observations; tantôt il trouve l’inclinaison de 4o” 48' 29", 
tantôt il la trouve de 4o° 4a' 16”, et ainsi des autres élémens. 

Il n’y a rien à conclure de-là contre l’exactitude, des calculs 
d’après lesquels nous venons de déterminer l’orbite de la comète 
de 1769; ces calculs donnent , entre toutes les paraboles pos- 
sibles , l’une de celles qui satisfont le mieux aux trois obser- 
vations données. Mais il est possible , et il paroit même très- 
probable, d’après ces différences, que la vraie orbite de la comète 
de 1769 n’est point une parubole. Il faudrait discuter de nou- 
veau les observations de celte comète, choisir les plus exactes 
après les avoir corrigées de l’aberration et de la parallaxe , et 
alors on pourroit essayer de délcrmiuer l’ellipse ou l’bypcr- 
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bole qui est'la Traie trajectoire. Le calcul ne serait pas beau- 
coup plus compliqué que les précédens; il y entrerait l’élément 

3c plus -, eltroisobservationssuffiroient toujours pour résoudre 

le problème ; mais il serait bon d’en (aire entrer au moins 
quatre en comparaison. 

C ,* « ' ' ♦ 

Observations sur la méthode précédente. 

LXI. Pour avoir une idée nette de la méthode que nous 
proposons , et dont nous avons détaillé les calculs dans un 
exemple, il faut observer qu’elle est composée de quatre par- 
ties distinctes. 

Première partie. Après avoir représenté le logarithme de la 
distance périhélie , et l’instant du passage au périhélie par deux 
expressions contenant chacune une correction indéterminée , 
on calculé , d’après ces valeurs , l’anomalie 4 et le rayon vec- 
teur r-pour le moment dè chaque observation. 

Seconde partie. Connoissant tout à-la-fois le rayon vecteur, 
la longitude et la latitude de la co/nète données par l’observa- 
tion , la longitude du soleil et sa distance à la terre données 
parles tables, on calcule la longitude et la latitude héiiocen- 
triques de la^oihète; et pour cet effet , il est bon de construire 
une figure qui représente à peu -près les positions respectives 
des différens points , et qui aidera le plus souvent à Axer l’indé- 
termination que présente dans certains cas le triangle CTS 
dans lequel on commît deux côtés , et l’angle opposé à l’un 
d'eux. . " . 

Le Résultat de ces deux premières opérations, dont on peut 
former an petit tableau , donnera , pour l’instant de chaque 
observation , l’anomalie 4 , la longitude héliocentrique s , et le 
log. sinus de la latitude héliocentrique a. Il faut ensuite procé- 
der à la détermination du plan de l’orbite et de la position du 
périhélie. ■ 
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Troisième parité. Ayant représenté i’écliptiquc par une ligne 
droite, et l’orbite de la comète, vue du soleil, par une ligne 
sinueuse; on marquera sur celle-ci les différcns lieux de la 
comète , d'après les longitudes et latitudes hélioccntriques déjà 
calculées. On y marquera également le lieu du périhélie qu’on 
doit déjà connoitre d’une manière approchée par l’une des ano- 
malies. On donnera ensuite des valeurs en partie connues, en 
partie indéterminées, à l’inclinaison I, et à la distance <r de 
la comète à l’un des nœuds lors de l’observation moyenne. 
Au moyen de ces valeurs, on calculera la latitude a, et la 
quantité «(qui n’est autre chose que la distance * réduite à l’éclip- 
tique) par les formules sia a = sinlsin <?, tan g * = cos 1 lange. 
La comparaison du log. sia a, avec celui qu’on a déjà trouvé 
dans la seconde partie , donne une première erreur E' sur 
l’angle a ; on en trouve une pareille sur les deux antres lati- 
tudes, ce qui fait trois erreurs E', E", E". 

Quatrième partie. Des trois valeurs de * correspondantes 
aux trois observations, on déduit trois expressions de la lon- 
gitude du nœud, qui égalées entre elles fournissent deux équa- 
tions entre les quatre coéflioicns inconnus. Tirant de ces équa- 
tions les valeurs de deux des coëfliciens , et les substituant dans 
l’expression des trois quantités E', E", E'", il_ ne restera plus 
qu’à faire en sorte que la somme des quarrés de ces trois quan- 
tités soit un minimum. Dc-là résultent deux équations qui 
achèvent de déterminer tous les cocfhciens, et par suite tous 
les élémens de l’orbite. 

LXII. On pourra procéder de même, si on veut combiner 
ensemble plus de trois observations j mais alors il conviendra 
d’établir deux systèmes d’erreurs. Le premier comprendra 
toutes les erreurs E', E', E'", &c. relatives aux latitudes; le 
second comprendra toutes les erreurs de la position du nœud ; 
si par les diverses valeurs de * calculées aux époques des diffé- 
rentes observations , on trouve , pour la longitude. du noeud, 
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les diverses valeurs S', S', S"”, &c. ; et qu’on appelle S' + e' la 
vraie longitude S du nœud , on aura le* erreurs successives 
é = e\ e" — e' + S' — S", e *' = e' + S' — S'", &c. , en nombre 
égal à celui des observations , et e sera une indéterminée nou- 
velle à joindre anx quatre déjà employées. Dans le système des 
erreurs^', e" , e", &c., on déterminera e' avec deux dcsaulies 
inconnues, de manière que la somme des quarrés de ces erreurs 
soit un minimum ; on aura ainsi l’expression des trois incon- 
nues en fonction des denx autMB . et la substitution de ces 
valeurs étant faite dans les quantités E', E", E"', &c. du pre- 
mier système, il n’y restera plus que deux inconnues. On 
déterminera enfin ces deux inconnues, par la condition que la 
somme des quarrés des erreurs E', E*, E", &c. soit un mini- 
mum , et on en déduira.les élémens de l’orbite. 

Il seroit plus direct de donner des valeurs en partie connues, 
en partie indéterminées , aux cinq élémens de l’orbite , et 
ensuite de calcnler pour le moment de chaque observation, la 
longitude et la latitude géoccntriques de la comète. On curn- 
pareroil les longitudes et les Jatiludcs calculées avec les lon- 
gitudes et les latitudes observées, cequi donneroit pour chaque 
observation deux différences on erreurs, et on égaleroit ensuite 
à un minimum la somme des quarrés de toutes ces erreurs. Mais 
le calcul conduit de cette manière seroit plus long que celui que 
nous avons indiqué , et il n’en résulterait qu’un avantage assea 
petit pour la diminution des erreurs, à cause du peu de varia- 
tion que subissent les quantités lorsqu’elles sont voisines du 
minimum. 

Observations sur le calcul des corrections indéterminées. 

LXIII. Les tables de logarithmes contiennent tout ce qui 
est nécessaire pour le calcul des corrections indéterminées, 
ainsi qu’on l’a vu dans les exemples précédens. Cependant il 
est bon de remarquer qu’on peut calculer dans tous les cas les 
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coéfficiens des corrections , sans reconrir an* différences tics 
tables qui sont toujours affectées de quelqu’crreur. Voici quel- 
ques préceptes à ce Sujet. 

Si l’on demande le logarithme du nombre a + Ax, clans 
lequel bx est une correction toujours supposée très-petite par 

rapport À a, le logarithme cherché sera log.a + x 

mêlant le module o.434ag, &c. Mais comme les logarithmes 
ont ordinairement sept décitsales, pour que la cocrection soit 
exprimée en unités décimales du septième ordre , il faudra, 
au lieu de m, mettre M — 434*945 , dont ielog. est 6 . 6377843 , 

et on aura log. (a-fAx) = a + x 



bm\ 
a ) ’ 


Réciproquement, ayant log. y — log. a + x( A) , on en tire 


le nombre y = a -+- x 


© 


Si bx représente un nombre de minutes serrant de correc- 

• . % 

, | .. * , |. | _ I O' » f/i 

tion a l are a\ et qu on appelle n le nombre constant 7 — , 

• r 

( r' étant le nombre de minutes comprises dans le rayon), dont 
le logarithme est 3. ioi5io4 , on aura 

log sm (o+Ax) = log si n a + x(nbcot.o ) 
log cos ( a+bx ) = logeas a — x(nb tanga) 

logtang(a + bx ) = log tanga + x (— — — ^ 

00 \sin a cos ar 

log cot. + Ax) = log col. a — x (— — ^ 

\ sm a cos a / 

formules où les quantités qui multiplient x sont réduites en 
unités décimales du septième ordre , de même que les loga- 
rithmes. Réciproquement si on a 

/ c \ * 

log tin y = log tin a + * (e) , il en résulte y = a + x I — long a J 
log cot y = log cot a + x (c) , il en résulte y — a — x rot. a ^ 
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log tangy = log tan g a + x (c) , 
log coi. y = log cot. a -j- x (c) , 


il en résulte y 
il en résulte y ~ a 



êin a cos a 
s in a cos a 


) 

) 


la correction ajoutée à la valeul* de a étant exprimée en 
minutes. 

Il est souvent utile de trouver le log. cosinus, lorsqu’on a 
le log. sinus, ou vice versâ, sans être obligé de chercher l’arc. 
Or par les formules précédentes si l’on a • 


log tin jr— log tin a +*(c)> il en résulte log cot y — logeât a — x (c long* a) 
log cot y = kg cot a-fjr(c) , il eu résulte log tin y — log tin a — x (c eut .' 1 a) 

A\cc ces formules et autres semblables , qn’on pent employer 
• suivant les circonstances , on s’habituera aisément au calcul 
des corrections indéterminées qui peut être fort utile dans toutes 
sortes d'approximations. 
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APPENDICE. 


Sur la Méthode des moindres q narrés. 

+■ 

Dans la plupart îles questions où il s’agit de tirer des mesures 
données par l’observation , les résultats les plus exacts qu’elles 
peuvent offrir, on est presque toujours conduit à un système 
d’équations de la forme 

E s= o + ij + cy + ft + &c. 

dans lesquelles a , b , c , f, &c. sont des coëfficieps* connus , 
qui varient d’une équation à l’autre , et x,y, z, crc. sont çjcs 
inconnues qu’il faut déterminer par la condition que la valeur 
de Esc réduise, pour chaque équation , à une quantité ou nulle 
ou très-petite. 

Si l’on a autant d’équations que d’inconnues x ,y , z, &c. , 
il n’y a aucune difficulté pour la détermination de ces incon- 
nues , et on peut rendre les erreurs E absolument nulles. Mais 
le plus souvent, le nombre des équations est supérieur à 
celui des inconnues, et il est impossible d’anéantir tontes les 
erreurs. 

Dans cette circonstance , qui est celle de la plupart des pro- 
blèmes physiques et astronomiques , où l’on cherche à déter- 
miner quelques élémens importans, il entre nécessairement de 
l’arbitraire dans la distribution des erreurs , et on ne doit pas 
s’attendre que toutes les hypothèses conduiront exactement aux 
mêmes résultats ; mais il faut sur-tout faire en sorte que les 
erreurs extrêmes , sans avoir égard à leurs signes , soient ren- 
fermées dans les limites les plus étroites qu’il est possible. 

De tous les principes qu’on peut proposer pour cet objet , 
je pense qu’il n’en est pas de plus général , de plus exact, ni 
d’une application plus facile que celui dont nous avons fait 
usage dans les recherches précédentes , et qui consiste à rendre 
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minimum la somme des quarres des erreurs. Par ce moyen , il 
s’établit entré les erreurs une sorte d’équilibre qui empêchant 
les extrêmes de prévaloir , est très-propre à Caire connoître 
l’état du système le plus -proche de la vérité. 

La somme desquarrés des erreurs E*+ E'*+ E®*+ &c. étant 

» (a + bx 4- cy + f z + &c.)* 

+ ( a + b' x -t- c'y + f z + &c.)* 

+ (a* 4- b" x + c'y + f'x + &c.)* 

+ Sec. j • • * 

si l’on cherche son minimum., en faisant Varier x seule, on 
aura l’équation 

, o =fab + xfb‘ +yfbc + zfbf + &c. , 

dans laquelle par fab on entend la somme des produits sem- 
blables ab + a'b'+ a‘b"+ îkc. ; par fb * la somme des quarres 
des coëfficiens de*, savoir i* + À” + b"’ + &c. , ainsi- de suite. 

lie minimum, par rapport ky, donnera semblablement 

o —fac + xfb c +yfc' + xffe + &C. , 
et le minimum par rapport à x , 

o =faf+ xfbf + yfcf+ zff' + &c., 
où l’on voit que les mêmes coëfficiens fbc,fbf, &c. sont com- 
muns à deux équations, ce qui contribue à faciliter le calcul. 

En général , pour former l’équation du minimum par rap- 
port à l’une de » inconnues , il faut multiplier tous les termes 
de chaque équation proposée par le coefficient de l’inconnue dans 
cette équation , pris avec son signe , et faire une somme de tous 
ces produits. . 

On obtiendra de cette manière autant d’équations du mini- 
mum , qu’il y a d’inconnues , et il faudra résoudre ces équa- 
tions par les méthodes ordinaires. Mais <on aura soin d’abréger 
tous les calculs, tant des multiplications que de la résolution, 
en n’adinettant dans chaque opération que le nombre de chiffres 

io 
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entiers ou décimaux que peut exiger le degré d’approximation 
dont la question est susceptible. 

Si par an hasard singulier, il étoit possible de satisfaire à 
toutes les équations en rendant toutes les erreurs nullcs, on 
obtiendrait également ce résultat par les équations du mini- 
mum ; car si apres avoir trouvé les valeurs de x, y, x , &c. 
qui rendent nulles E, E', &c., on fait varier x, y, s, &c. de 
Sx , Sy, Sz , &c., il est évident que E’ qui étoit zéro deviendra 
par cette variation ( aîx 4- bSy + cSz, &c.)\ Il en sera de 
même de E'*, E'* , &c. D’où l’on voit que la somme des quarrés 
des erreurs aura pour variation une quantité du second ordre 
par rapport à Sx, Sy, &c. ; ce qui s’accorde avec la naluredu 
minimum. 


Si après avoir déterminé toutes les inconnues x,y, z, &.O., 
on substitue leurs valeurs dans* las équations proposées , ou 
counollra les diverses erreurs E, E', E", &c. auxquelles ce 
système donne lieu , et qui ne peuvent être réduites sans aug- 
menter la somme de leurs quarrés. Si'parmi ces erreurs il s*cn 
trouve que l’on juge trop grandes pour être admises, alors on 
rejettera les équations qui ont produit ceB erreur*, comme 


venant d’expériences trop défectueuses , et on déterminera les 
inconnues par le moyen des équations restantes, qui alors don- 
neront des erreurs beaucoup moindres. Et il est a observerqu’on 
ne sera pas obligé alors de recommencer tous les calculs j car 
comme les équations du minimum se forment par l’addition des 
produits faits dans chacune des équations proposées, il suffira 
d’écarter de l’addition les produits donnés par les équations qui 
auront conduit à des erreurs trop considérables. 

La règle par laquelle on prend le milieu entre les résultats 
de différentes observations , n’est qu’une conséquence très- 
simple de notre méthode générale , que nous appellerons 
Méthode des moindres quarrés. 

En effet , si l’expérience a donne diverses valeurs a’, a ", a'", &c. 
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pour une certaine quantité a , Insomnie des quarrés des erreurs 
sera (a' — i)’ + (a" — *)’ ■+■ (a'” — xf + &c. , et en égalant 
cette somme à un minimum , on a 

o = (a* — *) + (a* — *) + (a"' — x) + &c. ; 


,, . , . a + a*+ a'^ + ^c. 

d ou rcsullp x = | 


n étant le nombre des 


observations. 

Pareillement, si pour déterminer la position d’un point dans 
l’espace , on a trouvé, par une première expérience , les coor- 
données o', b', c ; par une seconde , les coordonnées a", b", c\ &c. 
ainsi de suite; soient x ,y, z , les véritables coordonnées de ce 
point : alors l’erreur de la première expérience sera la distance 
du point(o c) au point (x,.y,.i); Je quarré de cette 
distance est 

(«'-*)*+ (*'-*)* + («>'— O*; 

et la somme des quarrés semblables étant égalée à un minimum, 

. , . , fa fb fc 

on en lire trois équations qui donnent x = — , y = — , x — — , 

Tl Tl IX 

n étant le nombre des points donnés par l’expérience. Ces for- 
mules sont les mêmes par lesquelles on trouverait le centre 
de gravité commun de plusieurs masses égales , situées dans les 
points donnés ; d’où l’on voit que le centre de gravité d’un 
corps quelconque jouit de cette propriété générale. 

Si on divise la masse d'un corps en molécules égales et assez 
petites pour être considérées comme des points , la somme des 
quarrés des distances des molécules au centre de gravité sera 
un minimum. 

On voit donc qnc la méthode des moindres quarrés fait 
counoitre , en quelque sorte , le centre autour duquel viennent 
se ranger tous les résultats fournis par l’expérience, de manière 
à s’en écarter le moins qu’il est possible. L’application que nous 
allons faire de celle méthode à la mesure de la méridienne , 
achèvera de mettre dans tout son jaar sa simplicité et sa 
fécondité. 



application à la mesure- des degrés du méridien. 

Supposant que le méridien est une ellipse dont les axes sont 
dans le rapport de 1 à 1 + «t, si on désigne par D la longueur du 
45'" degrc, et par S celle de l’are compris entre les deux lati- 
tudes L et L', on aura par les formules connues , et en expri- 
mant L' — L en degrés : 

O 

S = D(L'— L) — i«D.— si/x(L'— L)cos(L'+L){ 

• V I 

d’où résulte 

L — L = + f«.— tin (L'— L )cot (L'+,L). 

Comme le 45*”' degré est d’environ a85oo modules, égaux 

chacun à deux toises , on peut faire = ■ 1 ■ - , C étant une 

D a85oo 

fraction très-petite , et on aura 

L '- L “ » 

équation qui pour chaque arc dont on connoît la longueur, avec 
la latitude d.e scs extrémités , donnera une relation entre « et C. 

Voici maintenant les longueurs des différons arcs de la méri- 
dienne de France et les latitudes des parallèles qui les séparent, 
telles qu’elles résultent de l’opération exécutée par les célèbres 
astronomes Delatnbre et Méchain. 


Lieu 

de T observation. 


Panthéon a 


Evaux, 


Carcassonne 


Monljouy 


Digitized by Google 


Sa latitude. 

Arcs compris 
exprimés en modules. 

L' — L 1 / + L 

5i* a’ io'5o 
48 5o 49.75 
46 10 4 a. 5o 

43 ta 54.4o 

• 

4 ■ » « 44 80 

DP 62472.59 

PE 76145.74 
£C 844*4-55 
CM 52749 

2°i 1' ao" 75)99° 53' 0" 1 
a 40 7 .a 5 g5 1 3a 
2 57 ^-îolüg a 3 37 
1 5i 9.60:84 3 , 39 
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Nous avons donc quatre arcs dont les mesures étant substi* 
tuées successivement dans l’équation ( a ) , fourniront quatre 
équations entre « et £. Mais comme ces quatre équations ne 
peuvent pas être satisfaites toutes à-la-fois , nous supposerons 
qu’elles ont lieu , en attribuant une certaine erreur à la latitude 
de chaque lieu , et nous appellerons E' , E", &c. les corrections 
addilives aux latitudes de Dunkerque , du Panthéon , &c. Ces 
erreurs n’entrent que dans le premier membre de chaque équa- 
tion : elles sont trop petites pour affecter le terme multiplié 
par a. dans Iq second membre. Voici donc les équations qui 
résultent des quatre arcs mesurés dans l’opération de la méri- 
dienne , 

E' — E" = 0.002923 + £ (2. 192) — « (o.563) 

E" — E'"= o.oo3ioo + £(2.672) — « (o.35i ) 

E"'— E ,r = — 0.001096 + £(2.962) 4- a. (0.047) 

. E” — E T = — 0.001808 + f(i.85i) + a(o.a63) 

Comme il importe de considérer les erreurs séparément, on 
regardera comme une nouvelle inconnue l’erreur E"', par 
exemple , et on aura les cinq équations : 

E' = E" + o.oo6oa3 + £(4.864) — a (0.914) 

E' = E'" + o.oo3ioo + £(^.672) — a (o.35i) 

e"'= e" . . y,,, *<*> 

E''= E"' + O.OOlOyli — £(2.962) r— a (0.047) . , . , 

E' = E"-f 0.002904 — £(4.8i3) — a (o.3io) 

Il faut maintenant faire en sorte que la somme des quarrés de 
ces cinq erreurs soit an minimum, et d’abord cette condition 
exprimée par- rapport à 'l'inconnue E ,( dont tous les coëfficiens 
sont 1 , donne par l’addition de toutes ccs équations : i • * 

o = 5E" + o.oi3i23 — £(o. a3g) — «t ( 1 .622 ) 
donc E"'= — o.ooa6a5 + C(o.o48) ■+■ a(o.3»4). 

Substituant cette valeur dans les cinq équations (4), on aura 
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K' = 0.003398 + f(4-9«a) — «(0.090) 

E" = 0.000^75 -+■ £ (□ .720) — • « (0.037.) 

E"= — o.ooa 6 i 5 ■+• f (0.048) -t- « (o. 3 a 4 ) (c) 

E” = — o.ooi 5 ag — £(9.914) •+* « (0.977) 

E' es O.POO979-— £(4.765) + « (O.0l4) 

Pour exprimer ensuite la condition du minimum par rapport 
à £, it faut multiplier la première équation par 4,91a coeffi- 
cient de £ ; la seconde, par 3.720; la troisième, par o.o48; 
la quatrième, par — 2.914; la cinquième , par — 4.765, et 
égaler à zéro la somme de tous les produits. On opérera sem- 
blablement par rapport à « , et on aura les deux équations : 

► 

0= o. oaoçj83 + £(62.726) — *(3.83o) 
o = — 0.003287 — £( 3.83o) + * (o.53i) 

d’où l’on tire « = 0.00675, et £= 0.0000778 , donc. 


l’aplatissement a = 
et le 45 **” degré D = - 


i48 
2K5oo 
i + f 


28497.78. 


L’aplatissement déterminé par la longueur du pendule et par 
quelques .phénomènes astronomiques, n’est que de -jy-„ , et le 
45 *“* degré tel qu’on Fa déduit de la comparaison des mesures 
faites en France avec les mesures faites au Pérou , estde 28604. 1 o. 
C’est sur ce dernier résultat qu’est fondée la détermination défi- 
nitive du mètre : il devroitêlre diminué d’environ un 45 oo*“*, 
si on s’en tenoit aux seules mesure» exécutées en France ; mais 
l’aplatissement -pÿj trop peu d’accord avec celui que l’on con-r 
noît par d’autres phénomènes , ne permet pas d'adopter ce 
dernier parti. 

Les valeurs trouvées pour « et £, déterminent l’ellipse qui 
satisfait aussi exactement qu’il est possible aux mesures de l’arc 
du méridien compris entre Dunkerque et Barcelonne. Celte 
ellipse est beaucoup plus aplatie que celle qui convient à U 
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figure générale du globe; elle suppose dans les latitudes obser- 
vées des erreurs que l’on déterminera en substituant les valeurs 
trouvées pour « et C, dans les expressions de E', E", &c. : on 
trouvera , en réduisant ces erreurs en secondes , 

E'= — o" 7 3 , E'=i* 83 , E"= — i' 55, E"=o'’4a, E’=o'o3. 
La plus grande de toutes ne monte pas à a*, et la moyenne, 
sans égard aux signes , n’est que de o"9i. 

Si au lieu de chercher les deux quantités * et C qui con- 
viennent au minimum absolu , on commence par faire la quan- 
tité a égale à l’aplatissement connu j— j , les équations (ç) 
deviendront 

E' = o.ooi554 + C( 4.911 ) 

E" = o.ooo3gi +C(2.7io) 

E"= — 0.001611 + C ( 0.048 ) 

E" = — o.ooo663 — C( 3.914) 

E* = o.ooo3a3 — 0(4.765) 
et on aura pour l’équation du minimum, o— ô. 009010+ C (61.716), 
d’où résulte C = — 0.0001 436 ; donc 

le 45*"* degré = a85oo (1 — C) = 28504.09 ; 
ce qui s’accorde suffisamment avec la détermination adoptée ; 
mais alors les erreurs E', E*, &c. exprimées en secondes, 
deviennent 

f;=f 3 *o 6, E" sa o'Vro , «"=—5*83, E"=— o"88, E*=3’6i.' 

Ces erreurs sont plus fortes qu’elles n’étoient dans le cas du 
minimum absolu ; la plus grande tombe sur la latitude d’Evaux, 
et la moindre, qui est même entièrement nulle , sur celle du 
Panthéon . • “ 

Au reste, les anomnlies'dnns les latitudes, qui sans au&iti 
doute ne doivent point être attribuées aux observations , tien- 
nent vraisemblablement à des attractions locales qui agissént 
irrégulièrement sur le fil à plomb. Il suffît, pour cela, d’un défaut 
d’homogénéité dans les couches qui avoisinent le point où l’on 
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observe la latitude ; et In même cause qui rapproche le zénith 
apparent du midi ou du nord , peut aussi le détourner de quel- 
ques secondes vers l’est ou vers l’ouest ; ce qui explique les 
inégalités qu’on a aussi observées dans les azimuths. 

Il résulte de l’existence bien constatée de ces anomalies , que 
la longueur dea arcs du méridien est moins propre que celle du 
pendule j à la détermination d’une mesure universelle; et il 
n’est pas étonnant que des observateurs, d’ailleurs très-exacts , 
ne sc soient pas accordés dans les mesures qu’ils ont prises des 
degrés du méridien , puisqu’à raison des attractions locales , les 
latitudes de deux lieux également éloignés de l’équateur , pou r- 
roient différer entre elles de plusieurs secondes. 


Paris , le 1 5 ventôse an j3. 
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SUPPLÉMENT 

AUX NOUVELLES MÉTHODES 

POUR LA DÉTERMINATION 

DES ORBITES DES COMÈTES. 


La Comète découverte par M. Bouvard, le ao octobre der- 
nier, m’ayant fourni l’occasion d’appliquer les méthodes expo- 
sées dans le Mémoire précédent , je me suis apperçu que dans 
certaines combinaisons des observations, les élémens de l’or- 
bile, calculés d'après les équations (35) et (36), n'étoient 
pas suffisamment exacts. J'avois pensé d’abord qu’il ponvoit 
8’êlre glissé quelque erreur dans la réduction des observations ; 
mais eu y réfléchissant davantage, je n’ai pas tardé à recon- 
noitre que le peu de succès de la méthode étoit dû à la 
petitesse du coefficient C , et par suite , à celle des coeffi- 
ciens P, Q, H, avec lesquels se forme l'équation (35) ; il 
est évident en effet, que si, par les circonstances de la ques- 
tion, le coefficient C est ou nul ou trop petit pour Être déter- 
miné assez exactement par les observations , alors l’équation (35) 
cesse d’être suffisamment approchée , et sa combinaison avec 
l’équation (56) ne donne plus que des résultats incertains. 

On peut remarquer que le coefficient C est nul , lorsque 
les lieux apparens de la Comète, dans la première et dans la 
troisième observations , sont situés dans tin même grand cercle 
avec le lieu du soleil, au moment de la seconde observation; 
et c’est ce qui a eu lieu à fort peu près dans les observations 
dont nous avons parlé. 

Puisqu’alors on ne peut obtenir la vraie solution par lez 
formules que nous avons indiquées, il faut reprendre l’analyse 
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(*) 

générale du problème , et examiner ce que deviennent , dan# 
Je cas de C très-petit , les deux équations entra fi, u et a, d'où 
l’on avoit déduit le rapport de fi k fi, et qui ne sont exactes 
qu’aux quantités près de l’ordre fi*. On trouve que dans ce cas 
ces deux équations sont très-pou différentes l’une de l’antre ; 
desorte que par l’élimination de m , les coeffieiens de fi et fi 
devenant très-petits, les fermes restans sont confondus avec les 
erreurs de l’ordre ê*, et il n’y a plus de conclusion certaine à 
en tirer. 

Mais quoique les deux équations entre fi , fi, et u ne permet- 
tent pas d’éliminer t» ; néanmoins on peut toujours en déduire 
la vraie solution , en combinant l’une de ces équations ou 
une sorte de milieu pris entr’elles, avec les autres équations 
du problème; caria réduction accidentelle de deux équations à 
une seule , laisse encore subsister un nombre d’équations suffi- 
sant, puisqu’en général le problème des comètes, dans l’hypo- 
thèse parabolique, contient une équation de plus que d'inconnues. 

Nous nous sommes proposé , dans ce supplément , le double 
but, et d’obvier ainsi à la difficulté qui n’avoit pas été prévue 
dans le Mémoire précédent , et de simplifier encore assez 
considérablement la solution générale du problème.', c’est à 
quoi nous croyons être parvenus par différons moyens / et 
surtout en évitant l’emploi de coeffieiens trop petits. 

Par cette nouvelle méthode , nous avons déterminé de deux 
manières différentes, les éléraens approchés de l’orbite de la 
première Comète de i8o5. Nous avons cherché ensuite à corriger 
ces élémens, pour avoir une orbite qui représentât, aussi exac- 
tement qu’il est possible , les observations données. 

Cette recherche nous a donné lieu de perfectionner, à plusieurs 
égards, la méthode qui avoit été indiquée pour le même objet 
dans le Mémoire précédent , et surtout d'employer sons une 
forme plus commode , les corrections indéterminées qui peuvent 
être d’une grande utilité dans diverses sortes d’approximations. 
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ANALYSE DU PROBLÈME. 


x. V-iOJfSBRVANT les mômes dénominations que dans le 
Mémoire précédent (art. I, II...VI), et faisant de plus 


.m'— M'=V4*, «*— N*=/, rn’—W^fS, /T— N*=»“, etc., 
on aura en général ; 

x — etc. 
y— Y= » + »'/+»'/* -f- »V-f- etc. ; 
et l’équation — Y =(x-~-X) tang a , dans laquelle cm sub- 
stituera successivement les valeurs /= — 8, /= o, t—tf, donnera 
les trois suivantes : 

» — ►'8-f-r'ô* — s'fl’-f-efc .=_/■• (fi — /a'0+/u'â* — M^+etc.) 

('»') 

, 4 . »' 8 ' 4-»V4-r'6' J 4.etc.=/’(|t4- At '6'-fM'â ,, +^V J +etc. ) 

Les quantités 8, U ,f' yf , étant connues, on peut faire pour 
abréger 


i^!=r » 

a * 


£±£=f, 

=F, 


a 


(*0 


ce qui donnera 

6'=t-J-t', S=r-r',/'=F+F'T,/’=F-Fr. 

Substituant ces valeurs dans la première et la troisième des 
équations (a ") , et omettant les quantités du quatrième ordre 
en r et ■/ , on en tire les deux suivantes : 

» + ►' ’'+'•('■+>'•) 1 _ f FCM + /*V+/»-(T»4-,'*)+t‘V(3T*4t , ‘)] 
+ f — \ +r7-C^+^V+ # i*(T»43.'*)] 

a ,V 4-»«(T*+3r'M -f F 'C#‘-f#‘'7 , +j>**t^+’'*)D 


( 4 ) 

a. Pour donner une forme encore plus simple à ces équa- 
tions , faisons 

v— Fp=y, v'—F/x'=ÿ, y" — F >*'=/, y’—Ffi—y’ , 
elles deviendront : 

>+>V + »' <»•+✓*) l_ rT . c / , V v 
+j,V J T* 7 + a /“ T ) 

>' + +>• (T*+ 3r'*)=r (f< +,*" ( T* + ï'*) ) 

d’où résulte, en éliminant y\ 

i'*) + aj,V(r*— t , *)=F'(v*— v") (ft'+f»V)— FV#*. 

Or on a y-j-F''r / /i=r — F/u+F't7*=(/— -F-f-FV)^ soit done 

la quantité connue 

• ' \ (cj 

et l’équation précédente, divisée par t*— t'*, donnera 

F V + >' + 2>*r' = F' (/*'+/*' t'). 

Celte équation ne devant être exacte qu'aux quantités près 
de l’ordre r* , puisque l’équation primitive ne contient point: 
les termes de l’ordre r 4 , il faudra pareillement supprimer les 
quantités de l’ordre r*, dans l’équation en y'-, desorte que nous 
n'aurons plus à considérer que les équations réduites, 

FV+/+ 2 > v=F' o*'-h«v) 


Y + a*' t' = F(m+J»V). 


(<0 


3. Maintenant on 

a ( art. 

III) 


- 
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ce qui donne 
H-' 
/S 


(5) 


= M *> 


= — T ïNû», 




~ 6?‘ + — s N '» H- f N £ , 


et il en résulte 


~ ^ - i ( N'-FM') 4- ? (N-FM). 

Substituant ces valeurs dans les équations (a f ), remettant en- 
suite au lieu de y' sa valeur / - F/x' , et au lieu de y sa 
valeur t — Ffi, ou (/— F)/*, on simplement — F'r>, parce- 
que/— F + t'F est du second ordre, d’après i’équation (c') 
on aura les deux équations suivantes : 

F"n— F (fi— ‘-£— iMr'.) = (^— (y) (N— PM) V 

<,, .. +t<n-— pm-)| 

4- Fjir 4- FV/*'af- T'a (N — FM). 

Ces deux équations jointes à l’équation >=>, expriment les- 
conditions relatives aux trois longitudes de la Comète. 

4- De même l’équation des latitudes ‘ 

<= Sr f (— X); 

où l’on fera successivement r=_ fl, t=o, r=fl', f ournira des 
résultats entièrement semhlables à ceux qu’on vient d’obtenir;, 
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c'est-à-dire qu’en faiaant pour abréger , 



on aura le» trois équations p = jft , 

GV-C' M/-)= ("-</) ( — GM) — £ GM' 

(V) 

p' = G/ + G> 4. G'; y — GM/». 

A ces si* équations il faut joindre l’équation donnée par le 
triangle S CT , savoir: 

7 - = m’ + «• + p* = ( M -f- m)* 4- C N -f 4-3>‘ , 

ou 

r* = R* + N/) + (i 4-/* 4-*’) , 

laquelle revient à l’équation (56), et l’équation du mouve- 
ment parabolique 

4* *'* 4- p\ 


5. H est visible que les équations (<d) et (g") se simplifient 
beaucoup en faisant x':=o, c’est-à-dire en supposant que le» 
intervalles de temps entre les observations sont égaux. Ce cas, 
auquel on peut ramener tou* les autres , à l’aide de l’inter- 
polation, est celai qu’il convient de discuter avec le plus de 
soin. 

Soit donc t'ato, ou r=ô's=S, on aura d’abord entre/*, 
p,' et a, les deux équations : 

Ppc — F/*' = i a» (N — FM ) 

vo 

G>— G>'=— >GM, 

lesquelles peuvent donner lieu à diverses combinaisons qui in- 
flueront diversement sur l’exactitude de la solution. 
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Si ces équations étoient rigoureuses, de quelque manière 
qu’on les combinât, on parviendrait toujours à une solution 
exacte. Mais ces équations ne sont qu’approcbées ; d’une part 
elles supposent l’omission des termes de l’ordre fl* qui y entre- 
raient, si on avoit conservé ceux de l'ordre 6* dans les voleurs de 
X >X> 1 > d’autre part, comme les observations sont sujettes à er- 
reur, les coefficicns F et G sont affectés de cette erreur j les coeffi- 
eiens F’ et G' le sont encore davantage à cause de la division 
par fl, qui a Jieu dans les valeurs: 



enfin les coefficiens F", G', déduits des formules 


F' = 


/ F 

t* » 



» 


sont ceux sur lesquels les erreurs des observations ont le plus 
d’influence. 


6. Maintenant si on élimine p' des équations (A') pour avoir 
l’expression de p. en «, savoir i 




/(N— FM)G'+GMF'\ 
\ PG'— F’C" 


il est érident que cette combinaison est sujette à être inexacte, 
tant pareeque le dénominateur F'C» — F'G' contient les eoei- 
ficiens les plus imparfaits, que pareeque ce dénominateur 
pouvant se réduire à une quantité très-petite ainsi que Je nu- 
mérateur , les termes négligés pourraient devenir comparables 
au numérateur restant |»[(N — FM.)G'-f-GMF']. Cette com- 
binaison revient à celle des équations (ag) ; elle est la plus 
facile pour la détermination des quantités p et r; mais elle 
est la plus sujette à erreur, et il convient généralement de 
n’en pas faire usage, non plus que des équations (35) et (34), 
qui sont sujettes au même inconvénient. 


( 8 ) 

7 . Si, en second 'lien, on chasse « des équations (A'), 
pour en tirer le rapport de y à n, on obtient 


, /-C*(N— FM) + F*CM\ 

M P ' \G' (N — FM) FOM/‘ 


CO 


Cette combinaison conduit à l’équation (35), et elle est de 
beaucoup préférable à l’autre, puisque le dénominateur est 
composé de quantités qui en général sont déterminées d'une 
manière plus précise. C’est pourquoi nous avons eu raison de 
présenter les équations (35) et (36) comme devant conduire au 
résultat le plus exact. 

Cependant si dans l’expression précédente le dénominateur 
G'(N — FMJ+F'GIVl , qui répond au coefficient C, devient 
très-petit, le numérateur doit devenir aussi très-petit, et les 
résultats de la formule cessent d’être exacts par les raisons 
que nous avons déjà données. 

Le cas dont il s’agit est celui où les deux équations (A') dif- 
fèrent très-peu entr’elles, les rapports entre les coefficicns de 
// et eu étant presqu’égaux dans les deux équations. Alors on 
ne peut plus éliminer avec sûreté l’une des inconnues de ces 
deux équations, et il faut conserver dans l’expression de y les 
deux inconnues n et a>\ ce qui n’à d’ailleurs aucun incon- 
vénient, puisque le problème des combles offre une équation 
de plus que d’inconnues. 

8 . Si donc l’une des équations (A') donne y = hfi -f- la , 
et que l’autre donne y=hy~\-tai, les coefficiens A' et l 
étant peu différens de h et l, on pourra choisir entre ces deux 
valeurs celle qui est déduite de l’équation où y a le plus grand 
coefficient ; mais il vaut encore mieux tirer des deux équa- 
tions un résultat moyen dont l’erreur soit la plus petite pos- 
sible. Pour cela , supposons qu’on ait les erreurs 


F. = FV — F> + i a. (N— FJV1) 

E'= Gy — Gy — > . GM ; 

et 
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et cherchons une Valeur de p! qui rende la somme des qnarrés 
de ces erreurs égale à un minimum , nous trouverons 

_ (F'F' +G*G> + i »f (N-FM)F'-MGG') 
r F*» + G'* 

Cette valeur, que nous représenterons ainsi 

p' = Hp-j-Lu, ( i ’) 

est celle qu’il convient d’adopter. 

9* Cela posé , voici la marche à suivre pour parvenir 
aussi simplement qu’il est possible, à la solution du problème. 
On donnera une valeur arbitraire à p pour servir de première 
hypothèse , et on aura soin de prendre son signe conforme à 
celui de cosa cosbj ( car on a p = p cos a cos b , et p est tou- 
jours une quantité positive ) ; on calculera la distance r par 
l’équation du triangle SCT qui est 


'-=h-+K*- S£ S =:! )+ 


— Ü .. 

cos* a cos* o* 


<*) 


r étant connu , on calculera m par la formule 

(0 

ensuite on substituera les valeurs de p et et dans l’équation (/*) 
ce qui donnera la valeur de p\ 

Au moyen des valeurs de p' et p, on calculera les quan- 
tités m‘, n', p , par les formules 

m' = M' -f- p! 

»'=N'+F^ + F> (b0 

P =G/t'-f- G' p. 

Enfin substituant ces valeurs dans l’équation du mouvement 
parabolique 


s 
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2 = m'* + «'•+//*, (n') 

on en tirera une seconde valeur de r qu'on comparera à la 
première, afin d’avoir leur différence, qui sera l’erreur delà 
première hypothèse. 

Lorsqu’on aura fait deux hypothèses qui donneront deux 
petites différences ou erreurs sur la valeur de r, il sera fa- 
cile d’en conclure, par les règles connues , les vraies valeurs 
de /x et de r, et par suite celles de m' , n', p , Ou calculera 
enfin m, n , p, par les valeurs 

» = N +fp. (/>') 

P=gP) 

et à l’aide des formules <4a), (45), (44), (58), (5g), <4®), (47), 
on déduira aisément de toutes ces données les élément de 
l'orbite. 

10 . Cette méthode devient plus simple, si au lieu de la 
formule (r*) on peut employer avec sûreté la formule (f),qui 
donne directement le rapport de p à p-, mais pour cela il 
faut que son dénominateur ne soit pas trop petit. Alors cette 
méthode coïncide entièrement avec celle que nous avons don- 
née dans le Mémoire précédent, laquelle se réduit à la réso- 
lution des équations (55) et (3fi). Mais la nouvelle forme sous 
laquelle nous venons de la présenter, semble préférable à la 
première, en ce qne le» ooefficiens sur lesquels on opère se 
déduisent plus facilement des données de l’observation , et ne 
sont point sujets à être d’une petitesse qui rende leur calcul 
difficile ou inexact. Nous rappellerons ici, pour l’usage de ces 
nouvelles formules , qu’au logarithme de S , compté en jours 
du temps moyen, on doit ajouter le logarithme constant 
8. a555âar. Il est <bou que B n'esoède pas 5 à 6 jour», c'est- 
à-dire qu’il D’y ait pas plus de io à ia jours d'intervalle entre 
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la première et la troisième observation ", afin que ô* n’excède 
pas 737, et que l’approximation soit plus exacte. 

Quant aux ooefficiens F, F', F', G, G', G', ils se calcu- 
leront immédiatement par les formules 


f=£±C, 

C = £±£, G'=«ÿ, 



W) 


où les quantités f, g,f, çf , etc. sont les mêmes que dans 
les équations (7). Enfin les quantités M, N, M', N', relatives 
à l’orbite de la terre, devront toujours être déterminées par 
les formules de l’art. XL 

11. Observons encore que , dans la vue de simplifier les 
calculs autant qu’il est possible , on peut prendre à volonté 
l’une des quatre longitudes cC, a , a', A , pourvu que les 
trois autres conservent la même différence avec celle-là, 
et qu’on tienne compte de ce changement dans le résultat , 
en rétablissant l’origine des longitudes à sa vraie place. On 
peut donc faire, par exemple, A = 90°, ce qui simplifie le 
calcul des quantités M, N, M', N', dépendantes du lieu de la 
terre ; mais de toutes les hypothèses de ce genre , celle qui 
convient le mieux à la nature du problème, consiste à faire 
a=î=o, ce qui donneray = o, cos a = 1. Alors cosa‘ et cos a 
différeront peu de l’unité •, les quantités f° et f ne seront pas 
susceptibles de devenir trop grandes, non plus que g", g et g', à 
moins que la latitude de la Comète ne soit très-près de 90* dans 
l’une des observations -, et les équations du problème formeront 
deux systèmes presqu’indépendans , l’un pour le mouvement 
en longitude , l’autre pour le mouvement en latitude : en même 
temps la quantité F devenant de l’ordre 0 *, on pourroit faire 
F =0 , tant dans les équations (A') que dans les équations (/«') j 
ce qui les simplifieroit d’autant. 

ia. Après avoir traité avec toute l’étendue convenable le 



( 


( ) 

cas des observations équidistantes , revenons à la solution géné- 
rale du problème. Nous observerons d’abord que, quoique les 
intervallesS et O'soient inégaux, cependant on doit supposer qu’ils 
le sont peu , ou que leur différence est petite par rapport à l’un 
d’eux. Car si ces deux intervalles de temps étoient fort inégaux j 
comme il est nécessaire pour la convergence des séries que cha- 
cun des temps 9, 0', n’excède pas 5 à 6 jours, il y auroit donc 
deux observations fort rapprochées entr’elles, par exemple deux 
observations qui seraient éloignées d’un jour seulement ou de 
moins d’un jour. Concevons alors un nouveau Heu déterminé 
par interpolation entre ces deux lieux voisins ; toute orbite qui 
passera par le lieu moyen et par le lieu le plus éloigné , satisfera 
à-peu-près aux deux observations voisines l’une de l’autre ; de- 
sorte que le problème deviendrait indéterminé , ou qu’au moins 
on pourrait changer assez considérablement les élémens sans 
cesser de satisfaire aux observations. Il faut donc modifier en 
ce sens ce que nous avons dit sur la solution générale du 
problème ( pag. 7 du Mém. préced. ) , et supposer r' assez petit 
par rapport à r. 

i3. Cela posé, en examinant les équations (<?'). et ( g ), 
on y remarque une inconnue £ qui n’affecte à la vérité que 
les termes de l’ordre t', mais qui ne laisserait pas de compliquer 
beaucoup le calcul, si elle n’étoit pas éliminée. Effectuant donc 
l’élimination, et observant que MN' — M'N = 1 — £<*, ou sim- 
plement = 1 , on aura 


, fi ( MGF* — MFG’-t-NG'> — f «t'G , a 

*= «(F'ti— FG'j-TW (r') 


formule qui sera suffisamment exacte tant que le dénominateur 
M(F'G — FG')-|-NG', ne sera pas trop petit. Si le cas con- 
traire avoit lieu , on serait obligé de conserver £ dans les deux 
équations qui déterminent ft , et qui se réduiraient à une seule : 
mais comme £ n’affecte que des termes très-petits, on négli- 
gerait ces termes dans une première approximation , et on les 
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rétablirait dans nne seconde , en calculant la valeur de £ par 
la formule ( 5 ). Le calcul deviendrait ainsi beaucoup plus com- 
pliqué , et il seroit préférable , en pareil cas , d’interpoler les 
lieux observés pour en avoir trois équidistans auxquels la mé- 
thode précédente seroit applicable. 

14. Faisant donc abstraction de ce cas particulier, qui 
répond à celui où C est très-petit, la marche à suivre pour 
parvenir k la solution sera à-peu-près la même que ci-dessus. 
On déterminera les coefficiens F, G, F', etc. par les formules 



on fera une première hypothèse sur la valeur d e p., d’où l’on 
déduira les valeurs de r et », comme ci-dessus. Substituant 
ces valeurs dans l’équation (r'), on connoîtra /*'; ensuite on 
calculera m', n', pf , par les équations 

= + (O 

n' = N' 4-(F+FV)a'+F>4-t'«(N — FM), 
p' = (G + GV)/t'-fG>— • fs.GM. 

Et la substitution étant faite dans l’équation (n'), on aura une 
seconde valeur de r, dont la comparaison avec la première 
donnera l’erreur de la jpremière hypothèse. Lorsqu’on aura fai* 
ainsi deux hypothèses qui donneront de petites erreurs , on en 
déduira bientôt la vraie solution (*). 


(*) Cette solution s'accorderait avec celle qu'on a donnée dans le Mém. 
précéd. , n° XXXI , sous une forme différente , s'il ne s’étoit pas glissé dans 

I* 

celle-ci une erreur provenant de la suppression du facteur 1 4- dans les 
équations (8). Ce facteur ne peut être supprimé que lorsqu'on a V = t , on 

r 1 = o : dans les autres cas il produit dans la valeur de y.' le terme TÉrdonta 
faut tenir compte. 


( i4 ) 

• ! ' ' , ' .... 

Application « la première Comète de i8b5. 

i5. Nous allons calculer l’orbite de cette Comète , au moyen 
des observations suivantes que M. Bouvard a eu la complaisance 
de nous communiquer , et dont les époques sont fixées en temps 
moyen à l'Observatoire de Paris. Nous avons joint à ces ob- 
servations la longitude du soleil et le logarithme de sa distance 
calculés pour les mêmes époques. 


Époques. 

Longitude. 

Latitude bor. 

fjongitude du @ 

Log R 

! 

Octobre. 

D 

M 

S 

D 

M 

s 

D 

M 

s 


2 . . 6q5aa 

l 60 

l 9 

ao 

23 

Sa 

9 

208 

ao 

36 

9 997539 

22.68488 

.63 

ao 

53 

aa 

5 9 

53 

209 

*9 

5o 

9- 9974 » > 

30.68670 

■ 83 

48 

3a 

.5 

37 

ai 

317 

>9 

48 

9 998498 

3i . 6q5a3 
Novemb. 

■ 85 

5 4 

53 

>4 

42 

5i 

218 

ao 

25 

9 -998384 

3.7*088 

> 9 * 

4 S 

t5 

13 

a 

*9 

aai 

àa 

a6 

9.996053 


La lacune qu’on remarque dans ces observations entre le 
33 et le 3o octobre , exige qu’on calcule par interpolation un 
lieu assez éloigné des lieux observés. Ainsi voulant faire usage 
des observations des 3o octobre et 3 novembre , qui sont à un 
intervalle de 4 jours , il faut calculer le Heu équidistant du 
a6 octobre. Ce calcul peut se faire soit par les lieux des 3i, 
aa et 3o octobre , soit par ceux des 33 , 3o et 3» ; et les résultats 
étant peu différons, on prendra le milieu entr’eux. Voici donc 
les trois observations par lesquelles nous devons déterminer 
l’orbite. 
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Époques. 

Longitude. 

Latitude bor. 

Lieu du 

£^)gR. 

Octobre. 
a6.65:j5a 
3c . 68670 
34.72088 

D M S 

D M S 

D M S 

9-99 6 49S 

174 90 46 

1 83 48 5s 

>9> . 4S i5 

'9 s3 14 
i5 37 si 
1 a a 29 

317 i 9 48 


16. Avant d'appliquer la nouvelle méthode exposée dans 
ce Supplément , nous allons d’abord donner les résultats de 
l’autre méthode *, et comme il convient de ne s'occuper que 
des équations (35) et (36) qui sont les plus exactes, voici, 
d’après les données précédentes, les logarithmes des coefficiens 
et des quantités qui servent à les fermer : 

aRcwc 0.202179, y..... 3.336634, 

P 6 . 54 o 455 , 

(— Q) 7.096148, 

( — H) 6.3299274 

Cela posé, les équations à résoudre seront 

t* = 0.983993 — p(2Rcorc) ~ bp ‘ , 

i = o. 508101 — p£L'] +p*[L'], 

L coeff. . . . 9-856883. . . . 0.608574 ; 

1 : . . . 

On en tire la solution 

/û^^=^.79p5*3, 

/pg-p = 9.791894^- ' 


1— 

___ 


o.oo3443. 


tsirrir 8.172008, 

C 7.8o568i, 
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et de là les logarithmes des coefficient 

( — rri) o. aa3753, m 9.387074, 

(—ri).... 9.775141 , n 9.749670, 

( — /)■••• 9.450455, p 9.333137, 

(mri — m'n) 9.916578, 

(mjf — m'p) 9.349515, 

— (nff — rip) 8.793158. 

D’après ces valeurs on voit que la quantité mri — m'n est posi- 
tive et mm' ~\-nri +pp' négative; d’où il suit (XXIV et XXV ) 
que le mouvement de la comète est direct et qu’au 3 o octobre 
elle n’avoit pas encore atteint son périhélie. On trouve ensuite 
pour les éiémens : 

Log. dist. périb 9. 565 1 58 . 

Pass. au périh novembre 18.2097. 

Inclinaison i 5 ° 4 3 ' 6\ 

Lieu du nœud asc 344 a 9 * 5 . 

Lieu du périh. sur l’orbite i 5 o 8 14. 

Mouvement . direct. 

Si avec ces éiémens on calcule le lieu de la Comète aux 
époques de la. première et de la troisième observations , ou 
trouvera les résultats suivans : 

36 oct. lloDgit... 174° ia' 38*1 latit. . . . 19* 38 

3 nov. I 191 54 37 I 11 59 a6 

d’où l’on voit qu’il y a environ 8' d’erreur sur chacune des 
longitudes , et a ou 3 ' seulement sur les latitudes ; ce qui est 
une approximation passable , eu égard à la petitesse des coef- 
ficiens C , P, Q , H. Venons maintenant à la seconde méthode. 

ty. En se conformant à la remarque du n* n, on retran- 
chera la longitude i 83 * 4 &' 3 a* de chacune des deux autres et 
. de 
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de celle de la terre , et on aura les longitudes fictive* : 


a* = 55 o* 5 a' i 4 ' » 

a =s= o o o , vf 5= ai 5 ° 5 i' i6* , 
a' = 7 57 43 * 

Les quantités A*, A , 4 ', restant d’ailleurs les mômes, ainsique 
log R , on trouvera les logarithmes des coefficiens et des quan- 
tités auxiliaires, comme il suit: 


8 8.841357a 

(—/*)•• 9.2218700 

/' 9.1457046 

g * g.55a5756 

g 9.446581a 

g g. 333 aai 5 


( — FJ. .8.137388 

F.. o. 544 n 8 

F . .0.444613 

G. .9.456478 
(•—G').. 0.007995 
(— G")> .0. *26578 


(— M).. 9. 91 7497 
(—N).. 9. 788637 
M'. .9.755135 
(—N'). .9.93018a 


Ces valeurs étant substituées dans les équations (ù'J , on 
en déduira deux valeurs de /*', 

p! =/t(i.a6o36) + »(o.iai5i) 
ft,' = ft. ( 1 . 3 1 3g6) •+• m (o . 1 16 1 3). 

Il est visible qu’on ne pourroit éliminer avec sûreté a» da 
ces deux équations, parceque le rapport de y! à n se dédui» 
roit de coefficien» trop petits. C’e st le cas où C est très-petit, 
et où les deux équations se réduisent presqn’à une seule. Pre- 
nant donc un milieu entr’elles (*), on aura 

(*) Il anroit été plu* exact de se serrir de la formai* (i') qui donne 
/»' = /*( 1 .08375) +« (o. iao 5 S) , 

ou de choisir entre les deux valeurs de // , celle qui est déterminée par de 
plus grands coefficiens , c'est-à-dire la première. Mais le calcul étant déjà 
fait d'après la valeur moyenne , on n'a pas cru devoir le recommencer sur nnc 
antre valeur de y! , d'autant que les résultats, qui ne peuvent guères différer 
l'un de l’autre, ont toujours besoin d'être corrigés ultérieurement. 

3 


( i8 ) 

fl ' = /* (i .28716) -f- » (o. 1 1882) l 
0.10963a 9.074890 / 

équation qu’il faudra combiner avec les suivantes : 

r* =0.983993 — fi [a*] -f-M’O*] 1 

0.318537 0.032696 | 

«=^-— 1 .02450 

m'=M'-hfi' ' * ' ' 

n'= N'-f F/*' + ï> 

P' = Cfi' ■+■ G '/i 
* = m'* •+■ »'• + t?'*. 

Or d’après la marche que nous avons indiquée (art. 9}, on 
trouve aisément la solution suivante: 

P- 9 - 77 3 943 

r 9.790963 

m' . ,. .0.234428 ( — m). .. .9.566934 

.9.667360 9. 738637 

(■— p"). . . .9.44^853 p ... .9.320534 

(mn’—m’n) 9 • 9*9939 

— (jnpf — rn'p) g . 345483 

( np ' — n'p) 8 . 872884. 

De là résultent les élémens suivans, dans lesquels on a 
ajouté iSS* 48' 3 a* au lieu du nœud, et à celui du périhélie , 
pour revenir des longitudes fictives aux longitudes vraies. 

Log. dist. périh 9.571698 

Passage au périh.... . .... Nov. 18.18747 

Inclinaison i 5 ° 38 ' 24* 

Lieu du nœud asc 345 6 5 i 

Lieu du périh. sur i’orbite 149 o 28 


Digitized by Google 


( «9 ) 

i8. Ces élémens ne diffèrent pas beaucoup de ceux que 
nous avons trouvés par la première méthode : pour juger 
de leur exactitude, nous avons calculé les lieux qui en ré- 
sultent aux temps de la première et de la troisième observations. 


a6 oct. longit. 174* 17' 54* latit, 19° a3' 4** 
3 nov. 191 49 10 ia 1 5a. 


L'erreur est donc d’environ 3' sur chaque longitude, et 
d’une demi-minute sur chaque latitude. C’est toute l'exactitude 
qu’on peut attendre d’une première approximation, d’autant 
que le lieu du a6 octobre n’étant pas un lieu observé , mais 
bien un lieu conclu , il peut s’écarter de la vérité de quelques 
minutes, et même s’en écarter dans un sens qui ne permet- 
troit pas de faire passer exactement une orbite parabolique par 
les trois positions données. On ne calcule pas le lieu du 3o oc- 
tobre , pareeque ces solutions satisfont toujours exactement à 
l’observation moyenne. 

Pour juger encore mieux des deux systèmes, nous avons 
calculé dans l’un et dans l’autre le lieu de la Comète à l’ins- 
tant de l’observation du aa octobre. Voici le résultat: 

I er Système long. \6a* £6' 5i* , lat, a3* a' 4»' 

a® Système i63 6 aa a3 o a 

Observation. . . .~.‘.TTr*~. . . «65 -a* 53 , aa 5g 53 


D’où l’on voit que les latitudes s’accordent assez bien , mais 
qu’il y a a4' d’erreur sur la longitude dans le premier sys- 
tème, et 14' i dans le second. Cette différence de 14'î peut 
paroitre considérable ; mais on verra qu’il faut changer assez 
peu les élémens pour corriger cette erreur-, du reste l’avan- 
tage de la comparaison reste toujours à la seconde solution. 

*9. Pour donner une seconde application de notre mé- 
thode , nous allons calculer l’orbite de la Comète d’après les 

. . ..If Uv/ - rr . • I • 1 
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trois lieux suivans , dont l’un (celui du 28 octobre) a été 
conclu par interpolation de ceux des ai, aa, 3 oet 3 i octobre. 


Epoques. 

Longitude. 

Latitude bor. 

Lieu du G ■ 

Log. R. 

Octobre. 

D 

M 

5 

D 

M 

s 



aa.68488 

1 63 

20 

53 

aa 

5 9 

53 



28 . 70388 

»79 

*7 

36 

■7 


5 a 

ai 5 * ao' 3 g* 

9 - 99 s 723 

34.72088 

* 9 * 

46 

i 5 

13 

a 

a 9 




Réduisant à zéro la longitude du 28 octobre > on aura les 
longitudes fictives: 

û*= 344° 3 ' 17* 

a = o o o A ss aifi* 3 ' 5 *) 

a = 12 28 5 g 

«onservant d’ailleurs les latitudes b *, b, b', telles qu’elles sont 
données, et faisant fl = 6^,018 , on trouvera les logarithmes des 
coefficiens comme il suit : 


fl. .g.oi5o34o 

(-/•). .9.4559286 

y' 5 44 94 74 

g°. .9.6448505 

s . . 9.4964562 
gf . .9.3393954 


(— F).. 8 . 5 o 8 o 54 
F* . . o . 588940 
F*. .0.477986 
G . .9.51844^ 
( — G') . . o . o 5 a 1 37 
(—G*).. 0.181779 


(— M). .9.904401 
(—N). .9.766472 
M'. .9.781829 
(~N'). .9.906122 


Par la substitution de ces valeurs les équations (A') donnent 

/fc' = A* 0 ' 4 ,o 7 ) + < 1 * (0.12294) 

A*' = A* (« .3276) + u ( 0 . 12454 ) ; 

d’où l’on voit qu’on ne peut éliminer a>. C’est donc toujours 
le cas de C très-petit, et il faut prendre une sorte de milieu 


* 
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entre ces deux Valeurs. Ce milieu, déduit de la formule (i*), 
donne 

ï = fi [ a' ] + « [ b' ] I 
0.099456 9.094403 ( 

Combinant cette équation avec les deux suivantes : 

« = jj -- 1 .033895 

r % = o.g 85 oa 4 — p- [ a' ] ■+- /u.* [ b* ] 1 
o. 30545 1 0.040754 j 

et avec les autres équations en m', ri, p', on trouvera la 
solution 

, N 

^..9.754545 j ni . .0.318973 I ( — m). .9.869514 


r.. 9. 8 i 5595 71'. .9.74346a ( — n).. 9.76647a 

pi. .o.oi 55 i 6 ( — p'). .9.453818 p.. g.a 5 iooi 

d’où l’on déduit les élémens de l’orbite que voici : 

Log. dist. périh 9.869548 

Passage au périh Nov. 18. 15098 

Inclinaison i 5 ° 4 a' 5 o <r 

Lieu du nœud asc 344 56 54 

Lieu du périh. sur l’orbite 149 35 ao 


Ces élémens ne diffèrent pas beaucoup de ceux que nous avons 
trouvés par le calcul précédent, '"mais ils approchent moins 
de la véritable orbite, ainsi qu’on peut s’en assurer en cal- 
culant par leur moyen les lieux de la Comète pour le 33 oc- 
tobre et le 3 novembre. Au reste on devoit s’attendre à une 
approximation assez médiocre dans le second exemple, i“. par- 
ceque les observations embrassent un intervalle de 13 jours, 
qui est un peu grand pour la convergence des séries ; a°. par- 
ceque le lieu moyen, celui du 38, qui devroit être le plus 
exact , est un lieu conclu sur lequel il peut y avoir plusieurs 
minutes d’incertitude. 
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Calcul de P orbite corrigée. 

ao. Ayant indiqué dans le Mémoire précédent diverses 
méthodes pour corriger aussi exactement qu’il est possible l’or- 
bite connue par une première approximation, nous nous propo- 
sons ici de perfectionner ces méthodes, et nous prendrons pour 
exemple la première Comète de i 8 o 5 , dont nous avons rap- 
porté ci-dessus cinq observations. Nous supposerons que les 
élémens corrigés de son orbite sont comme il suit : 

i*. log [1, ou log. dist. périh 9.571598-!- 

mtr est une correction indéterminée dans laquelle m est le 
module o. 434 a 9448 * Cette valeur signifie que la distance , dont 
le logarithme est 9.571598, doit être multipliée par 1 -f -ta-, 
<w étant une petite fraction , pour donner la vraie distance pé- 
rihélie que nous cherchons. 

3°. Instant du passage au périhélie.. .Nov. 18. 18747-+- r; 

r est une fraction de jour qu'il faut ajouter à l’époque trou- 
vée par la première approximation. 

3*. Inclinaison de l’orbite.... I=i5* 58' a4' + *i 

la correction a. est exprimée en parties du rayon , ce qui con- 
tribuera à simplifier les calculs. 

. ’ . < ni' • . .1 . r • • . 

4 °. Lieu du nœud ascendant S = 345 ° 6' 5 i'-f- 6 . 

5 °. Lieu du périhélie , moins le lieu du nœud ascendant 

compté sur l’orbite dans l’ordre des signes r=i 63 ° 53 ' 37 *-f- 7 ; 

les corrections € et y sont encore exprimées en parties dn 
rayon. • . 

31. D’après ces élémens, nous allons calculer le lieu de 
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la Comète pour l’instant de l’observation du 5 novembre, 

Temps donné Nov. 3.7308s 

Passage au périhélie Nov. 18.187474-1 

Différence t =. 14. 46659 4- x 

Le logarithme de 14.46659 est 1.160 366 a; le logarithme de 
a + x, x étant très - petit par rapport à a , est en général 
log a 4- mx ; ainsi on aura 

log ï = 1 . i 6 o 366 a 4- mr 

Retranchant de ce logarithme celui de FI* , qui est 
4 /ogn = 9. 35739704- 77tsr (1 .5), 

il restera 

log T = i. 803969a — mw (1 . 5) 4- mr (i); 

La partie connue de ce logarithme répond au nombre 63 . 5 a 86 : 
or en général log a+mx est le logarithme du nombre a + ax-, 
ainsi on aura 

T= 63 . 5386 — ar(iT)4-r(ï); 

, • . . . . t 

ou en effectuant fcnmWiwsU» 4 . . T et t, dans la partie in- 
déterminée , 

T = 65 > . 5 a 86 — (95 . 3939) 4-r (4. 3 9 i4), 

Or on trouve dans la Table des Comètes que ce nombre 
de jours répond à l’anomalie 6f 46' 4f, et qu’une différence 
d’un jour répond k une différence de a 584 * sur l’anomalie ; 
donc en convertissant ce nombre de secondes en parties durajoD, 

on voit qu’il faut multiplier les corrections par R' étant 


( a 4 ) 

le nombre de seconde* comprises dans le rayon. On obtient 
ainsi l’anomalie de la Comète 

4' = 67° 46' 47 * — <ar (1 . 1014) + t (0.05076) 
i^ = 35 ° 53 ' a 3 'j — <ar(o.55o7)-f-T (o.oa 538 ). 

Pour avoir le log. cosinus de cet arc , il faut suivre la for- 
mule log cos (a -f- x) = log cos a — mx tang a, qui prescrit de 
multiplier les parties indéterminées par — m et par la tangente 
de l’arc, qui est ici 33 ° 53 '. On aura donc par ce calcul 


logcos i-J. = 9.9191363 + mv (o .3699) — mx (0.01705) 
son double = 9.8383734 -f- m <9 (0.7398) — mx (o.o 34 io) 
log n = 9.5716980 + 7TC4T 

donc log r :=ng. 7333 a 56 -f- mtr (o.a6oa) (o.o 34 io) 

log R 9.9960550 

Diff. 0.3637374 — mit (0.3603) — mx (o,o 34 io) 

Cette différence est le log. tangente d’un angle auxiliaire P \ 
et comme on a loglanga-\-mxz=logtang(a-\-xsin a cos a) , 
l’angle dont il s'agit sera 

Ps=6i“ ai'39'.7— txsinVcosV (o.aôoa ) — rsinVcosV (o.o 34 io). 

Retranchant 45 * de cet angle, afin d’avoir un reste Q=P — 45 °, 
le log. tangente de ce reste sera, en vertu de la même for- 
mule ou de son inverse, log tang(a-\-x\s=log tang a-f- — wx — 

^ ^ sm a cotai 


_ iinPcoiP , . . iinPcojP, . 

J, tang Q=9 (o- ( o o3 4 10 ) ! 


jinQcoj<J v 


mais à cause de P = 45 ’ -f- Q , on a == c ol-*Q, ainsi 

en effectuant la substitution dans les parties indéterminées , 
on aura 


L tang Q = 9. 4677a j 9 — mxr (o . 4049) — mx (o . o 53 o 5 ). 

C’est 
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C’est le logarithme qn'il faudra ajouter à ltang\( C+T) 

pour avoir / tang j (C — T), comme l’exige la résolution du 
triangle CST. 

32. Maintenant la distance du périhélie au nœud étant 
i 63 ° 53 ' 3 7 * + y , 


si on en retranche l’anomalie 

^=67° / fi ' 47' — ar (1 . 101 4) + r (0.05076) , 
on aura la distance de la Comète au nœud , 

<r = cfi ’ 6' 5 o" + «-(i . 1014) — r (o.o 5 o 7 6 ) -f- >• 

Avec cette distance et l’inclinaison I, on calculera la la- 
titude héliocentrique X et l’arc x, différence entre la longitude 
liéliocen trique de la Comète et celle du nœud, par les for- 
mules sinh = sinlsina,tangx=cosïtang<T-, on aura d’abord, 
en suivant la formule Lîin(a + x) — Lsina-\-mxcota, 

L sin 1=9 .430707 3+//1 a. ( 3 . 5 7 ao) 

L ji«a , = 9.9975337 — //itv(o . 1 i8o)+7nT(o.oo544) — my(p. 1071) 

L sin X= 9 . 4 a 8 a 3 oo — rrnsrio . 1 1 8o)+//ir(o . oo 544 ) 

+/7ia(3.5 7 ao) — my(p. 1071) 


pliant par tang" X, on tire 

L cos X = 9.9838078 + wtar (0.0091) — T7iT (0.00043) » 
— ma. (o . 3 7 65) + my (o . oo 83 ) / 


L’autre formule tang x = cos I tanga , donnera 


L cos I = 9. 9836 148 — m*( o.a8oo) 

L( — langi) =0.g70i5a3 — m»( 1 o . 4 oo 4 ) 4- mr(o . 47929) — m y( 3 - 44 3 y) 

L(— tang*) = 0.9537671 — nv»(io.4oo4)+mT(o.47929) 1 
— ma . ( 0.2800) — m>(.9-44 2 ,9) > 


4 
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et il en résulte 

* =: 96’ ao' 5 o’ 4 <nr(i.i 4 * 7 ) — r(o.o 5 a 66 ) f 
4 * (o.o 3 o 8 ) 4 > (1.0375) ÿ 

Ajoutant à cette quantité la longitude du nœud 

545 “ & 5 i' + C, 

et retranchant la longitude de la terre 4 '° 33 ' 3 ® > on aura * 
l’angle 

TSK = 4 o» 5 ' rS'-f-'sr (1.1427)— TCo.o 5 a 66 ) + e 
4 * (o.o 3 o 8 ) -f- > (i.o375> 

Calculant ensuite l’angle TSC au moyen de la formule cos TSC 
cos A cor TSK, on aura 

/ogcoiTSK =: g. 8836966 — mv (0.9618) 4 rai (o . c4433) 

— mx (o.oaog) — mC (0.8417) — m‘)(p.'&j' 55 ) 
l,>gcos\=x 9.9838078 4 m* (o.cogi) —mi (0.0004a) 

— mx (0.2765) 4 m>(o.oo83) 

log cos TSC = 9 . 8675044 — m<r (o . 9617) + ">t ( ° • 0409 > ) 

_ mx (o.3ca4) — (0.8417) — ra> (o.865o) 

Ainsi dans le triangle TSC, on a l’anglcr 

g — 4 1 °^ 1 ’ 5 Ti -l^irTi . o 3 oo ) — t(o.o4?88) 

+ *(0.3398)46(0.9180)47,40.9434). 

De là résulte 90’ — i S, ou 

^±Z == 68 “ 44 , a 7 'a—'ar(o. 5 i 95 ) 4 r(o. 02394) 

— *(0.1649)— £(0.4590)— >(0.4717), 

Le log-tangente de cet angle est 

o . 4099814 — msr ( 1 . 5374 ) 4«t (0.07085) 

— m *(0.4880) — jb£( 1 . 3583 ) — tw>(i .39^9) , 


Digitized by Google 


( a 7 ) 

auquel ajoutant L tang Q déjà trouvé; savoir, 

9.4677219 — m'a- ( o . 4049) — tot(o.o 53 o 5 ), 

on aura log tang f (C — T) = 

9.8777033 — mtr( 1 .94a3)-4-WT (0.01780) 

— ma (0.4880) — mt > ( 1 . 3583 ) — my ( 1 .3959); 

donc 

îî^=37* 2' 1 4 ' 5 — ur (o . 9539) 4- t (o . oo 856 ) 

— * (0.2246) — £(o. 653 i ) — y (0.671a). 

Retranchant cette quantité de j(T + C), il restera 

T= 3 i° 4 a/ i 2 ’ 7 + , ® , (o. 4 i 44 ) - i- T ( 0 - ( >i 538 ) 

4 - a (o. 0697)4-^(0. 1941)4-7.(0. «995). 

a 3 . Maintenant la latitude b se calculera par la formule 

. , . .. .. 
sinb — — — - — , de cette manière: 
smb 

Lj//iT =9.7 20592 u-j-mtr(o . 6709)4-11^(0 . oa4go) 

+rai(o. 1 i29)4-/n€(o.3i4 a )4-w>(o.32a9), 
Ls/'ii * = 9 . 4282 5 oo — ma(o. 1 i 8 o) 4 -hit(o.oo 544 ) 

4 -to ai ( 3 . 57 20)— my(o. 1071) 

Somme — 4 s (o ■ 55 29 ) 4 - 1 // t(o . o3o34) 

~\-rna. ( 3 . 6849)4-w€’(o . 3 1 42)4-mj.(o. a i 58 ). 
Z.Æi/jS=9.8a9834o4-m<nr(i . i 33 a)— mr(o.o5a22) 

4-/na(o.5597)4-mC(i .ooia) 4 -wj>(i 0289) 

L sin £=9.5i 8988a — mm(p . 58 o 3 ) 4 -/«t(o . 08 a 56 ) 

4 -ni ai ( 3 . 3 a 5 a) — m£( o . 6870) — my(o. 81 3 1); 

donc la latitude cherchée , ♦ 

b= 12 9 i' 5 a*i — «- (o.i237)4-t(o. 01760) 

* 4-* (0.7087) — £(0.1464) — >(0.1733). 



( a8 ) 

Enfin il reste à calculer la formule cos STK = ££iZ • 

CO JO ' 

L 98166 — m«-(o.a56o) — mr(o.oog5o) 

— m a. (o.o 43 1 ) — 11. 6(0. 1 1 99) — my( o. 1 a 3 a) 
L co s b =9.99o3543+w-®’(o.oa64 ) — //jt^.ooô^S) 

— m <t(o. i5 io)-f-77if(o.o3 1 a)-f-"J>(o.o 3 f> 9 ) 


LcoiSTK=9.93946a4 — m<aip.28a^) — iiit(o.oo5j5) 

-\-m a. (o. 1 079) — mÇÇo. 1 5 1 1 ) — inyÇo. 1 60 1 )j 


donc 


STK=a9°33' 1 G"5 -+-<»• (o. 4979 )-}- t (0.0 1014) 

— a (o. 1904)4-6(0. aG65)4-J'(o.a8a4). 

Cet angle étant retranché du lien du soleil aai“ aa' 36' , il 
restera pour la longitude de la Comète à l’instant donné : 


a= 191*49' 9*5— «-(0.4979) — t(o. 01014) 

4- a (o. 1904 ) — 6(o.a6G5) — j.(o.a8a4). 

Riais la longitude a été observée de igi^G' i 5 ', et la latitude 
de ia“ a' 39', on aura donc les deux erreurs 

E = a' 54 * 5 — «(0.4979) — t(o. 01014) 

4- ot (o. 1904) — ë(o.aG 65 ) — j.(o.a 8 a 4 ), 

£ =— 30*9 — «-(o. 1 ^37 ) 4- t(o. 01 -Go) 

4- *(0.7087) — £(0.1464) — >(0.1733). 

34. On pourrait calculer avec la même généralité les 
lieux de la Comète aux époques des autres observations , et 
ensuite par la méthode des moindres quarrés on détermineroit 
les cinq coefliciens et, ë , y , r , « . Riais pour simplifier les 
calculs , nous égalerons à zéro les deux erreurs précédentes , 
ce qui exprimera que l’orbite satisfait exactement à l’obser- 
vation du 3 novembre ; et de cette manière noug aurons la 
valeur de deux coclliciens indéterminés en fonction des trois 
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autres. Le reste du calcul aura pour objet de déterminer entre 
toutes les paraboles qui satisfont à l’observation du 3 novembre, 
celle qui représente avec le moins d’erreur possible les autres 
observations. 


Faisant donc E = o et e = o, et prenant & et r pour les 
deux inconnues à déterminer par les trois autres , on trouvera , 
après avoir réduit les termes connus de E et e , en parties du 
rayon : 

l 

<sr = o.ooi 5 o 5 -f- a( i .oSaa) — € (0.6164) — y (0.6717.) 
r =0.01934 — a( 3 a. 88 o) -f- £(3.988) -+->.(5.137). 

Ainsi pour les calculs suivans on supposera 

n =9. 57a i 65 -{- nia (1 .o 5 aa) — mS( 0.O1G4) — m>. (0.67 17), 
et l’instant du passage au périhélie • 

Novembre i8.aoG8i — «( 3 a. 88 o)-j- 6 ( 3 .g 88 )+j-( 3 . 137). 

Les autres élémens contenant les corrections a, 6, y, resteront 
les mêmes. 

a 5 . D’après ces élémens qui ne contiennent plus que les 
trois indéterminées a , £ , y , si on calcule les lieux de la 
Comète aux époques des observations des aa et 5 o octobre, on 
trouvera les résultats suivans : 


Iongft. ftn^Üi) — C+o.Haÿi) — y (o.3i53) 

tarit. aa 59 3g. 3 4- * (o.oo+5) — C (o. 1740) — >-(0.0326) 

r longit. 1 83 ° + 5 ' 10* '-f « (o, i 3 + 5 ) — C (o. 1071) — y (o. ioao) 
ù ° c ' l latit. i 5 37 47 -f « (o.oig 5 ) — C (o.o 3 a 3 ) +• >.(0.0113). 

» 

Comparant ces résultats aux observations, on aura ces quatre 
• erreurs, où la partie connue est évaluée en minutes: 

E' — — i 8 '.a 33 -f- a (0.9348) — 

E' = — 3 .567 -f- a (o.i 345 ) — 

e’ = — O .aa8 -|- a (0.0045) — 

e' = o .433 -f- a (o.oig 5 ) — 

tr 


£ (0.3371) — y (o. 5 i 53 ) 
C (o. 1071) — y (o. ioao) 
£ (o. 1740) — y (0.022(1) 
£ (o.o 3 a 3 ) -|- >.(0.011 a). 



t 
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Si l’on sp borne an calcul de ces trois lieux, il Faudra dé- 
terminer x, Ç, y, par ces quatre équations. Or d'après la mé- 
thode des moindres quarrcs, l'équation particulière qui déter- 
mine a est 

o = — !7'.3 o 8 -f- a (0.8758) — 6 (o. 3 i 83 ) — >(o. 3 o 52 ), 
d'où l'on tire 

a = 19'. 808 4- Ç (o . 3 G 43 ) 4- y (o . 34 g 3 ) j 

éliminant donc a au moyen de cette valeur, on anra 

E' = o'.o 85 4 * 6 (0.0098) 4 ~ >(0.0078) 

— E' = o .702 4- £(o.o 58 i) 4- >(o.o 55 o) 

— e’ = o . i 3 g -j - £ (o. 1724) 4 ~ > (0.0210) 
d = o .819 — Ç(o. 025 a) 4~ >(o.or8o). 

On voit immédiatement qu’en faisant €=o et >=o, les 
erreurs restent toutes quatre au-dessous d’une minute , ce qui 
est à-peu-près tout le degré d’approximation qu’on peut dési- 
rer daus la théorie des Comètes. En s’arrêtant donc à ce résul- 
tat, les élémens de l’orbite seront déterminés comme il suit: 

Log. dist. 'périh 9.574798 

Passage au périh Novembre 18.01736 

Inclinaison ; 1 5 ° 58 ' 1 a* 

Lieu du nœud asc....» 345 6 5 i 

Lieu du périhélie - 149 o 28 

Mouvement direct. 

26. Si avec ces élémens on calcule les lieux de la Co- 
mète aux temps des observations données, on trouvera les 
différences suivantes entre le calcul et l’observation : 


* 
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Er. en long. 

Er. en lot. j 

Octobre, ai 

- 33” 

60" 

aa 

— 8 

1 I 

5o 

— 33 

+ 4* 

3i 

+ >5 

— a 

Novembre. 3 

+ >3 

— 8 


On remarquera que le calcul a été dirigé pour satisfaire 
exactement à l’observation du 3 novembre, et que cependant 
il résulte des élémens trouvés, une erreur de lâ'surla lon- 
gitude, et de 8" sur la latitude. Cette différence est due à 
l'omission des termes des corrections qui contiendroicnt les 


quarrés et les produits deux à deux des coelHciens a, 6, y, 
r, /a-, et qui peuvent s’élever à plusieurs secondes, lorsque les 
corrections sont assez considérable* Or dans cet exemple il y 
a eu des changemens assez grands, puisque l’inclinaison est 
augmentée de près de ao'. Lorsque les corrections sont moindres , 
les seconds termes ont beaucoup moins d’influence •, au reste 
ces différences sont toujours très-petites en comparaison des 
erreurs qu’on fait disparoitre en tenant compte seulement des 
premiers termes. 

27. Ayant trouvé un résultat qui satisfait si bien aux 
observations, puisquela plus grande erreur n’excède pas 1', 
on peut se dispenser de toute recherche ultérieure. Cependant 
dans la vue de montrer l’usage de nos méthodes , nous allons 
tâcher de parvenir à un résultat 'encoie plus approché des 
observations. 


Si l’on se bornait à combiner les équations de l’art. a5, 
afin d’en tirer les valeurs de £ et y, par la méthode des moindres 
quarrés, on parvicudroit'facilemcut à représenter à moins d'une 


Or 
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demi-minute les trois observations des aa et 3o octobre, et 5 
novembre; mais alors l’erreur sur le lieu du ai octobre, qui 
est la plus considérable, augmenteroit encore. Pour remédier 
à cet inconvénient, il est donc nécessaire de calculer le lieu 
du ai octobre, en laissant 6 et y indéterminées, et conser- 
vant d’ailleurs tout le reste. On trouvera ainsi sur le lieu du 
a 1 octobre ces deux nouvelles erreurs : 

E" = — o’ .601 -f* Ç (0.0578) -f- y (o.o535) 

— e" = 1 .014 -f- 6 (o. 1943) -f- y (0.0339). 

Formant donc un seul système de ces deux équations et des 
quatre déjà trouvées, on aura pour l’équation du minimum 
par rapport à 6 ; 

0 = ao'.7a3 4-6(7.493) -f-y(i.6ia), 
d’où l’on tire 

6 = ■ — 3^765 — y(o.ai5i). 

Substituant cette valeur dans les six équations, elles de- 
viennent 

E* = — o'.o58 4- y (0.0057) 

E' s= — ■ 0.541 — y (o.o4a5) 

e’ — 0.338 -f- y (0.0161) 

ü = o .888 4 * y (0.0334) 

E" = — o .761 -f - y (0.0411) - 

e" = — o .477 4- y (0.0079). 

Celles-ci donnent pour équation du minimum , 
o = i3.84 4~ J' (4 • 4°) > 

d'où l’on tire y — — 3'.i5 et 6 = — a'.og. Les erreurs s’a- 
baissent de quelque chose par ces valeurs ; mais la plus grando 

montant 
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montant encore à o'.8g, et tombant sur une longitude, le 
premier système paroît préférable. 


Observation sur la méthode précédente. 


a8. La méthode que nous venons de suivre pour corriger 
les élémens de l’orbite , diflère beaucoup de celle que 
nous avons indiquée dans le Mémoire précédent. Comme 
elle nous semble préférable , nous allons en retracer les prin- 
cipaux points. 

i°. On fait varier à-la-fois les cinq élémens, en appliquant 
à chacun une correction indéterminée: cette correction est ex- 
primée en parties du rayon pour les angles, en jours ou frac- 
tions de jour, pour l'instant du passage au périhélie, et celle 
qui s'applique au logarithme de la distance périhélie est dé- 
signée par tnv , m étant le module des tables. 

a". Dans le cours du calcul, les parties des logarithmes dues 
aux corrections, se calculent par des formules directes, sans 
employer les différences des tables; on rapportera ci-après la 
suite de ccs formules. 

S». A l’aide des cinq élémens corrigés, on calcule la longi- 
tude et la latitude de la CorrièfiTpour l’instant d’une obser- 
vation choisie , à laquelle on veut satisfaire exactement. On 
égale donc à zéro les erreurs sur la longitude et la latitude , 
ce qui donne deux équations au moyen desquelles les cinq 
indéterminées se réduisent à trois. 

4°. Avec les élémens qui ne contiennent plus que ces trois 
indéterminées, on calcule la longitude et la latitude de la 
Comète aux temps des diverses observations. Il faut calculer 
au moins deux positions qui, avec celle qu'on a déjà employée, 
comprennent d'aussi grands intervalles qu’on pourra; mais on 

5 
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peut en calculer un plus grand nombre, et les résultats n’en 
seront que plus exacts. 

5 °. En prenantlesdifférencesdes longitudes et latitudes calculées 
avec les longitudeset latitudes observées, onaura les équations des 
erreurs qu'il faudra traiter par la méthode des moindres quar- 
rés , pour en tirer les valeurs des trois indéterminées. Et d'a- 
bord , considérait l’indéterminée qui a le plus d’influence par 
la grandeur de ses coefliciens , on pourra former, par rapport 
à cette indéterminée, l’équation du minimum , (ainsi qu’on 
l'a fait ci-dessus pour a), et éliminer cette quantité, ce qui ne 
laissera plus que deux inconnues dans les équations des erreurs. 
Oq, procédera ensuite de la même manière à l’élimination de 
l’une de ces inconnues, et enfin à la détermination delà troi- 
sième , qui fera connoître toutes les autres. Mais il arrivera 
assez souvent qu’on n’aura pas besoin de pousser l'opération 
jusqu’au bout, et que les erreurs se trouvant suffisamment at- 
ténuées dès la première ou la seconde opération , on pourra 
faire l'inconnue ou les inconnues restantes égales à zéro. 

29. Voici maintenant les formules dont nous avons parlé, 
an moyen desquelles on calcule directement les coefliciens 
de3 corrections. Il suffit de quatre décimales dans ces coeffi- 
cieng, et on doit e mployer pour ce calcul les tables à cinq 
ou six décimales. 

Appelons m le module 0.43429 etc., dont il n’est pas né- 
cessaire de substituer la valeur numérique dans le cours du 
calcul. 

1°. Si l’on demande le logarithme du nombre a-\-bx, dans 
lequel bx est une correction supposée très-petite par rapport 
à s, et t une iuconnue, on aura la formule 

log (cz — bx) = log a mx (jj). 

a*. Réciproquement ayant log y = log a •+- mx (b ) , on en 
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tire le nombre 

y — a -+- x (ah). i 

5*. Si on a l’arc z = a + x(b), il en résulte 

log sin z = log sin a -f- mx {b cota) 
log cos z = log cos a — mx {b tanga) 

log tangz = logtang a-\- mx 

log cot z = log cot a — mx ( -■ — - Y 

\sin a coj al 

4*. Réciproquement si l’on a 

* *T 

log sin z / sin a-f-mx(b), il en résulte z=a-f-z (b tanga) 
log cos z=l cos a-\-mx{b), il en résulte z=a — x (b cot a) 
logtangz=ltanga-\-mx{b), il en résulte z=a-\-r (b sin a cos a) 
log cot z=l cot a-\-mx{b), il en résulte z=a — x (b sin a cos a). 

Ces formules ne diffèrent point essentiellement de celles que 
nous avons données, art. LXIII du Mémoire précédent; mais 
leur forme est plus simple, en ce que les corrections ne sont 
exprimées ni en minutes, ni en décimales d'un ordre déter- 
miné , mais bien en parties du ra/on et de l'unité. Avec ces 
perfectionncmens , la méthode des corrections indéterminées 
pourra être d’un grand secours dans les recherches d’approxi- 
mation; ~ — 

Seconde Comète de i8o5. 


3o. Cette Comète , quoique très-petite , et à peine visible, est 
une de celles qui ont le plus approché de la terre, puisqu'au 
5 décembre sa distance n’étoit que de 7*- de celle du soleil. 
Aussi son mouvement a-t-il été assez rapide, tant en longi- 
tude qu’en latitude, et il a offert l’apparence d’un mouvement 
rétrograde continué pendant 17 jours, tandis que le mouve- 
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ment véritable étoit direct. A raison de cette proximité, le cal* 
cul de l’orbite exigeant une assez grande précision, nous avons 
cru devoir en faire l'objet d'une seconde application de nos 
méthodes. 

Voici quatre observations réduites de cette Comète , qui 
nous ont été communiquées par M. Bouvard, et auxquelles nous 
avons joint les lieux correspondans du soleil. Les époques sont 
fixées en temps moyen à l’Observatoire de Paris. 


Epoques. 

Longitude. 

Latit. boréale. 

Lieu du ©. 

Log. R. 

Nûv. l 8 . 3 qq-.J 7 

a3.3W,i 
3o.5iog5 
Déc. .S.aqôiii 

a 8 ° 14 ' i3" 

24 41 5 
>5 38 3G 
a 6 33 

29 * 3a' 17 ” 
27 a5 ig 
■ 9 u5 6 

3 19 sa 

a56“ g' 3 i*5 
241 8 23.3 
□48 a5 47.5 
a53 17 5.6 

9-9946‘4' 
9 -994^4 3 
q.nq56655 
9.9933784 


Par les lieux des 18 , a3 et 5o novembre, on en a calculé 
un nouveau qui fût à égales distances avec ceux des 3o no- _ 
vembre et 5 décembre. On a formé ainsi les trois lieux suivans, 
qui serviront à déterminer l'orbite. 


Epoque*.. 

l-onsinde. 

Latit. boréale. 

Lieu du ©. 

Log. R. 

Nov. 25.72809 
3o . 5 1 oq5 
Déc. 5 .s;|58t 

22° IO' l5' 

i5 38 36 
a G 33 

a5« a3' 4 a' 
19 a5 6 
.4 li) 2 -2 

a48’ a5' 47*5 

g.g 9 3G655 


La méthode du mémoire précédent n’étant, dans ce cas, 
sujette à aucune difficulté, si on en fait usage, et que pour 
simplifier le calcul , on retranche i5° 58' 56* de chacune des 
longitudes a*, a, a'. A, afin de rendre a nulle, on trouvera 
les Videurs logarithmiques qui suivent : 
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fl 8 .gi 545 i 

f sin Ÿ . . 7 . 938844 


(— C)... 9. 594187 


( 3 7 ) 

. . . .8.938715 


0 


Q 9. i 5 i 53 o 

H 9.09393a 

""*)•• .1.464937 


«C 


M 9.775268 

N 9.894790 

(— W')... 9. 910213 

N' 9.78295a 


Les équations ( 35 ) et ( 36 ) qu’il faut résoudre seront : 

a .... o . o 5 o 863 
a'. . . .0.415938 

a*. ... 1 . 262676 


r*=o.97ia5o+/> [a ] •+•/>* 
i =o.5i46y3+/>[a']+7>* [ a' ] 


et leur solution donne les valeurs logarithmiques 
f ■ . . 8 . 998666 


r. . .0.019365 


( — m '). . .0. ioi58i 

«'. . . .9.634745 

{-p')... 9 . 5575.5 


m .. .9.839953 
«... 9.894790 
p . . .8.528208 


(mn — ni n ). . .0.1 io 5 y 3 
— (tnp — ni p ) . . .9.31 6 iG 3 
— {np—n'p ). . .9.474058 


D’où résultent les élémens suivans , dans lesquels on a eu 
soin d’ajouter i 5 ° 58 ' 36 ", au lieu du nœud, et à celui du 
périhélie. 

Log dist. périh..' 9.g532i6 

Passage nwoiirili ...... Déc embre 5 o. 45687 

Inclinaison i5° 42 ' 3o* 

Lieu du noeud asc a 5 o 5 o 13 

Lieu du périh 168 9 12 

Mouvement * direct. 

Ces élémens, qui doivent toujours satisfaire à l’observation 
mojenne , ne représentent qu’à 20 ou 3 o* près les lieux des 35 
novembre et 5 décembre; ce qui n’est pas étonnant, vu la 
rapidité du mouvement de la Comète , et l’incertitude assee 
grande qu’il y a sur le lieu du a 5 , connu seulement par in- 
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ferpolafion. Nous allons maintenant nous occuper de la correc- 
tion de ces élémeus. 

Calcul de P Orbite corrigée. 

Si. Nous supposerons 


Log. dist. périh =9.953216 -f-mOT 

Passage au périh Décembre 30.43687 -f- 1008 (*) 

Inclinaison i 5 ° 4 a ° 5 o 


et d'après ces trois élémens , nous allons calculer les deux 
autres, de manière qu’ils satisfassent exactement à l'observa- 
tion du 5 décembre, la plus importante, soit par son exacti- 
tude , soit parceque la Comète étoit alors dans sa plus grande 
proximité à la terre. 

3 a. On trouve d’abord l'anomalie 4 et le rajon vecteur r 
comme il suit : 

4=38° 8' i4'6-f-0 (2.3790) — -ar(o. 8596) 

Log. r= o.ooaa 35 a-f-m 0 (0.7880) -f-mar (0.7039) 

Considérons le triangle STC formé par le soleil , la ferre et 
la Comète , soj t cnqim e dans la fig. I, du mémoire précédent,' 
K le lieu de la Comète projeté sur l'écliptique, et du point K 
imaginons qu’on abaisse une perpendiculaire KZ sur le rayon ST 
prolongé : cette figure très-simple suffira pour diriger le calcul. 

Puisqu’on a la longitude de la Comète a = a° 6' 33 *, et celle 
de la terre A=73* 17' 5 * 6 , il en résulte l’angle 


(*) On prend ici pour correction de l’époque 1 00 9 au lieu de 9 ou t , afin 
d'éviter les fractions trop petites qui se rencontreraient dans les multiplica- 
teurs de 9. 


Digitized by Google 


( 3 9 ) 

KT.Z= 7 i a 10' 5a*6. D’ailleurs la latitude b= 3° 19' aa\ De là, 
au moyen de l’équation cos CTZ = cos b cos KTZ, on ob- 
tient CTZ = 7 i°ia' îo^. Ensuite l’angle C du triangle SCT 

se calcule par la formule sin C = ! isin ^ 1 1 ‘ , qui donne 

Log. sin C = 0.0673544 — infl ( 0.7880) — nw ( 0.7030) 

C = 68*3'4a'— 8 (i. 9 564)— -st(i. 7 45i) 

Tl n'y a point d’ambiguité, parceque l’angle STC est obtus; 
on aura donc le troisième angle TSC ou 

S = 3» 8' 48'8 + 0 ( 1 . 9 5G4 ) + «T ( 1 . 7 45i ). 

La latitude héliocentrique X est donnée par la formule 

sin X= sin b , d’où résulte 

log sin X = 7 . 5a64644 ~i~ m9 ( 35.585) + war ( 3i . 7 4a). 
Mais en appelant <7 la distance de la Comète au nœud des- 
cendant, on a sin 7; donc 

log sin <r= 8 . 09391 1 a+;wfl(35 . 585 )+otw(3 1 . 7 4a) — ma(3 . 5550) 

o - = 0° 4a' 4°'6 — a (0.0441 ) +8 (0.44*8) + «•(0.3941 ). 

Retranchant <r de l’anomalie il restera la distance du pé- 
rihélie au nœud «Ucnamln if . nn l'élém ent 

T = 3 7 ° a5' 34* + a (o . o44 1 ) + 8 ( 1 .85 7 a ) — 1 ? ( 1 . a53 7 ). 

Cette quantité , ajoutée au lieu du nœud descendant , seroit ce 
qu’on appelle le lien du périhélie; mais il est plus simple de 
l’employer sous celte forme , et dans le calcul des autres lieux 
on aura toujours s , =4/ — T. 

Pour calculer le lieu du nœud , on cherchera d’abord l’angle 
KSZ par la formule long KSZ = cos b tang S . » qui 

donne 


( 4 ° ) 

L tang KSZ=8. 7393803 -)- m 9 ( 35.69a) 3 i .837 ) 

KSZ = 3 “ 8 ' 27*6 + 9(1.9528)-)- '»•( 1.7419) 

Reiranchant l’angle KSZ de la longitude A=-j$° 17' 5 " 6 , il 
restera la longitude héliocentrique de la Comète 

<p = 70” 8' 38 * — 9 ( 1 .9528 ) — ®- ( 1 .7419). 

-Enfin la distance x sur l’écliptique entre la Comète et le nœud , 
se calcule par la formule tang x = cos I tanga-, qui donne 

x= 4 »' 5*+9 (o. 4 a 53 ) -f-®- (0.3794) — <x (o.o 458 ); 

1 

et cette quantité étant ajoutée à la longitude <p , on a le lieu 
du nœud descendant 

t 3 =70° 4 ^ 45 * — * (o.o 458 ) — 9 ( 1 .5275) — iw ( 1 . 36 a 5 ). 

Tels sont les deux élémens T et ü qui» avec les trois sup- 
posés , satisfont généralement à l’observation du 5 décembre. 
Onauroit trouvé le même résultat , en faisant varier à la-fois 
les cinq élémens, comme nous l’avons fait pour la i rt Comète; 
et déterminant les deux élémens r et ü » de manière que l’ob- 
servation du 5 décembre fût satisfaite ; mais le calcul tel que 
nous venons de le faire , est plus simple» et pourra servir déréglé 



33 . Cela posé , à l’aide de ces cinq élémens qui ne con- 
tiennent que les trois indéterminées a, B, <rs , on calculera le 
lieu géocentrique de la Comète , pour l’instant de l’observa- 
tion du 5 o novembre. 

Voici le détail de ce calcul qui, à cause des petits angles qui 
s’y rencontrent , exige des attentions particulières. 

On trouvera à l’ordinaire l’anomalie et le rayon vecteur : 

4' = 44’9' 20 *+ 6(3,1107)— or (0.9475) 

log r = 0.0193616-)- m8 (o. 856 i )-) -m® (0.6157). 

Retranchant 
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Retranchant l'élément T de 4. , on aura la distance au nœud 
< 7 — 6 ° 43 ' 46' — « (0.044* )+ 9 (o. 3735 )+<®-(o. 3 o 6 a) j 

ensuite le calcul des deux formules sin X—sin I sin a , 

tang x = cos I tanga , donne 

log sin X— 8 . 5 o i 4 1 o 8 +wa(" . 1 8 ig)-f-rnfl(2 . 3 1 . 5 g 5 o) 

L tang%= 9 . 0553295 — rn«(o . G6oa)-t-;«8(2 . 35 oa)-|-w w(a . 63 1 a) 
x=6° 28' 49 * — a( 0.0740)-!- 6(0. aG 35 )-f- asip . ay 5 1 ) 

Ajoutant à x la quantité U , on a la longitude héliocentrique 

tp = 64* ao' 54 * + * (0.0283) — 8( 1 .7910) — «• ( 1 .6576) 

qui étant retranchée de la longitude de la terre 68* a 5 ' 47*5 , 
donne 

KST= 4 ‘’ 4 / 53 * 5 — a(o.oa 8 a) 4 - 0 (* . 79io)+'w(i .G576). 

Maintenant de log sin X on déduit, en multipliant les parties 
indéterminées par tang'X : 

L cos X=9 . 99978 1 3 — ma (o . oo 3 a) — mÜ(o . 0034) — o . 0026). 
Ajoutant logr, on a 

log SK =0 . o 191429 — ma (o . oo 3 a)-J-w 8 (o . 855 j)-i-mmr(o . 61 3 i). 
Ajoutant encore à celui-ci 

log cos KST = 9 . 998897s -f ma (o.ooao) — m!(o. 1078) — mst (o. u 83 ) 
on aura 

log SZ =0 . o 1 8040 1 — ma(o .0012) + m8(o . 73 59)+m-sr(o . 4948) . 

De là il faut déduire TZ=SZ — ST. Or en général si l'on a 
log y— log a -f- mx , il en résulte 

log{y — — i^ + logb+mx^^)-. 

6 


( 4 * ) 

par cette formule on trouvera 

log TZ=8. 755o6io — ra*(o .0220) +«6(i 3 . 3 oo)+m«( 9 . 066) ; 
mais on a KZ=SKs//iKST , ce qui donne : 

/g# KZ = 8.87 1 4756 — 7714(0.3984) + mfl (a 5 .g 53 ) -f- m-ar (33.843). 
Donc log KZ — logTZ , ou 

Lr<xng-.KTZ=o.i 164146 - mn(o.3764H-mS(i3 653)4-wi^(i4.777) 

KTZ= 5 a° 35 ' 30'— a(o. i 8 i 6 )+ 9 ( 6 . io 58 )+ «-(7 i 3 o 8 ). 

Cet angle étant retranché de 68’ a 5 ' 47 ^ * ^ restera la lon- 
gitude géocentrique de la Comète 

a=i 5 ° 5 o' 37 ' 5 +a(o. 1816)— 6(6. xo58)—«r(7.i3o8). 


Cette longitude a été observée de i 5 “ 38 ' 36 ", ainsi on aura 
l’erreur 

E'=s u'5i*5+« ( 0 . 1816 )— 9 (6.io58)— w(7-i3o8). 

La latitude b se calculera par la formule 

long £=j^.= sin KTZ » 

qui donne 


L tan g Z> =9 . 549a 796+7714(3 . 44 1 4 ) — w, ®( 1 ® ■ 1 ° 9 ) — mu ( 1 5 • 1 87 ) 
■ 4(1 .o85i) — 6(5.700)— «(4.780); 

mais la latitude observée est de 19-35' 6*, ainsi on aura cette 
seconde erreur 


c'=5' i4'+a(i .o 83 i)— 6(5.700)— «( 4 - 780). 

34 . Calculant semblablement le lieu de la Comète pour le 18 
novembre , on aura ces deut résultats: 

|E'= n'.858+*(o.i8a)— 6(6.106) — ar( 7 .i3i ) 

5oN °v*\ e '— 5 .a33+*(i.o83)— 8(5.700)— «(4.780) 

j E*=— 39 .437 + 3(0.393)— 6(5.388) — «(7.468) 
.8 Nov. | e ,_ — ,880+4(1 .6i5) — 6(3.385) — «(3.160) 
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auxquels il faut maintenant appliquer la méthode des moindres 
quarrés , pour rendre aussi petites qu’il est possible , les quatre 
erreurs. En vertu de cette méthode , l'équation du minimum 
par rapport à a , est 

o= — 4 f • «a* +*(3.968) — 0( i 4 - 863 ) — «-( 12.890) 
d’où résulte 

a=i 1 . 876+6(3. 746)+<sr(3. 249). 

Substituant cette valeur dans les quatre équations, on obtient 

E' = 14.024— 0(5-424)— ar(6. 540) 

e'= 18.095 — 8(1.643) — <ar( 1.263) 

E'= — 34.772 — 9(5.9 a o) — 'ar( 6 . 194) 
e'=. — 7.7 Oo+0(2.665) + w(5.o87). 

Maintenant l’équation du minimum par rapport à 0, est 

0=29.21 1 — 0( 54.%) — 3r (70.05), 

d’où résulte 

0=0.535 — ^(i .283). 

Eliminant 0 par cette valeur , on aura de nouveau 

E'= 1 1 . 121 +<sr(o.42i ) 

?= ) ..*^0.346) 

E'= — 26.870 — , ar ( i . i63) 
e'=— 6.274— ar(o. 333). 

Enfin l’équation du minimum par rapport à <o- , sera • 

0=52. 5g+ tr ( 2.356) , d’où résulte 'V— — 22'.3i 6 , et les 
quatre erreurs deviennent 

E'=i' 73, e'= — 1'66, E'=o'93 , c'= — i'i6} 

desorte qu’elles sont réduites au dessous de 2'; mais ce résultat, 
où l’on n'a eu égard qu'aux termes du premier ordre , ne sau- 
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roit être bien exact , attendu que les correction* sont assez 
considérables pour que les termes omis puissent s'élever à 
plusieurs minutes. 

<nr étant connu , on en déduit les deux autres corrections 
0 = 29' 1 74 > * = 48'66 i = 48' 4 °*> L’inclinaison corrigée est 
donc 

I=i6° 3 i' 10'. 

Quant aux valeurs de 8 et de <ar , comme elles sont expri- 
mées en minutes , il faut les diviser par le nombre de mi- 
nutes comprises dans le rayon. Ou trouvera ainsi 

mn — — 0.0028192, 1006=0.84864, ce qui donnera les élé- 
mens corrigés: 

Log. dist. périh 4 — 9.9502968 

Passage au périh Décembre 3 i.a 855 i. 

Les deux autres élémens F et ü pourroient se déduire des for- 
mules qui les représentent ; mais comme ces formules ne con- 
tiennent que les termes du premier ordre , il sera plus exact 
de les déterminer directement , de manière qu’avec les trois 
autres élémens ils satisfassent exactement à l’observation du 5 
décembre. On obtient ainsi 

r= 38 a 48' 5 o*" 
t3 = 70 3 a 16 

35 . Pour mieux comparer les résultats du calcul à ceux de 
l’observation , nous joignons ici un tableau complet des obser- 
vations de la seconde Comète de i 8 o 5 , corrigées de l’aber- 
ration et delà parallaxe, avec les lieux du soleil correspondant, 
calculés par les nouvelles tables de Delambre. 
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♦ ■ « 


Epoques. 

Longitude. 

Laf. boréale. 

Lieu du®. 

Log. R. 

Nov. i6.45ao3 
17-36386 
18.33937 
a3.3aa4i 
3o.5>oq5 
Déc. 5.ag58i 

39 ° a3' 37 " 
a 8 5i a3 
28 14 18 
34 41 4 
■5 3g 4° 
a 7 11 

3o° a' 40 ' 
ag 48 5a * 
39 3a 44 
37 a5 35 
îg a5 a 8 
3 ao 45 

a34° 1 1 ' 36" 
a35 6 5a 

a36 g 3g 
a4> 8 3o 

348 a5 4 1 
a 53 17 4 

9- 9947 8 47 
9- 9947°" 

9.9946166 

9 - 994 2 ° 7 ‘ 

9.9936673 

9-993364» 


Si d'après nos cinq élémens corrigés , on calcule les lieux de 
la Comèle pour le 18, le a 3 et le 5 o novembre, on trouvera 
les erreurs suivantes 


18 


/ E"=- 6 ' 41" , E'=- 6 ' 3 .* r E'=-3' 

l e* = — a' i3* 33 ( e*=— 3' a 7 * 30 l e' = — 5' 


E'= — 3 ' 4 a* 


où l’on voit que celles du 18 et du 3 o sont loin d'être au- 
dessous de 3', comme nous les avions trouvées ci-dessus. C’est 
sans doute à l'omission des termes du second ordre , qu’on doit 

attribuer ces différences. Eu effet , puisqu'on a trouvé 

33' 3 iG ; la seule variation de la distance périhélie 
produit dans li* uue différence = — -ar ( 7 . 468 ) = a* 46' 39". 
Cette différence est compensée, pour la majeure partie, par 
une tlifférMucc.senibljablc due à la correction Q , tuais il u’en est 
pas moins vrai que si les termes~ 3 irpremier ordre dont on tient 
compte sont si grands , ceux du second qu’on a négligés , 
peuvent être encore de quelques minutes. Il n’est donc pas 
étonnant que la plus grande des erreurs s’élève à 6' \i“. 

36 . Il importe maintenant de faire disparaître ccs nouvelles 
erreurs. Mais quelle méthode doitron employer pour cela? Re- 
commencera-t-on un nouveau calcul d’après les élémens trou- 
vés, joints à de nouvelles iudélerminées? Ce calcul devien- 
droit très-long et très-fastidieux , mais heureusement il est fa- 
cile de l’éviter. En effet les élémens corrigés ne diffèrent pas 
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assez des élétnens non corrigés , pour qu'on ne puisse regar- 
der les coefficieus des corrections , comme étant à-peu-près 
les mêmes dans les deux cas (*). 

Si donc par une seconde correction on prend 

log n =9 . 95 o3<j 68 -f- nvrt', 

l’instant du passage au périb.. . Décembre 3i ,a855t-f-roo9' 

et l’inclinaison 1= 16’ 3i' i o’+a' , 

les deux élémens r et 13 qui continueront de satisfaire à l’ob- 
•ervation du 5 décembre , serout à-peu-près 

r=38» 48- 5o*4-*'Co.o4)+9'C * -84)— « • a5) 

Î3=70 3a 16 — ai'(o.o5) — 6'(i.53) — <®-'(t.36), 

et les erreurs sur les observations des 3o et 18 novembre , se- 
ront aussi à-peu-près 

E'=s — 3.6974-*'(o. i 8) — &'(6. n ) — >zr' ( 7. i 3 ) 
e'=— 5 .oi 3 -f-*'( 1 .08)— ô'( 5 . 7 o) — «*'(4.78) 
E'=— 6 . 685 -t-*'Co. 3 9 )-fl'( 5 . 5 9 )— «-'(7.47; 
c’=. — a.222-i-a'( t .61 ) — Ô'( 5 . 58 ) — <ar'(a. 16). 

On y peut joindre les erreurs calculées semblablement pour 
l'observaTion OTTîH'rto'rernbre , lesquelles sont * 

E"= — 6.5a3 + *'(o.34 ) — l Y (6.07) — «■'(7.92) 
e*= — 3.45o+«.'(t .5a) — 8'(4.33) — <s-'(3, 17). 


(*) On se Fonde ici sur ce principe d’analyse : soient 9 et f’ deux fonctions 
semblables, l’une des variables x, y, z , etc . , l’autre des variables J , ÿ, s/, etc., 
peu différentes de celles-là. Si l’on a , 

ç{x+*, y+G, i-f-7-etc.) = ? + A<t + BC-f-Cj.+«tc.,j . 
on aura aussi à-très-peu-près 

f (xf -f-«e,y -1-C, *'+>» etc. ) =s»'-f-A*-f-BS-f-Cy-f-etc. 
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5>7> En appliquant la méthode des moindres qnarrés à ces 
équations, il seroit facile d’obtenir le résultat qui donne ab- 
solument les plus petites erreurs ; mais sans entrer dans au- 
cun calcul pénible , on voit immédiatement que les erreurs 
seront beaucoup diminuées en faisant a'=o,Ô' = o, et 
Alors on conserve l'inclinaison de 16’ 5i' 10', et le passage 
au périhélie Déc. 3i.a855i , et on ne change que le log. de 
la distance périhélie , qui devient 9.9503700. D’après ces trois 
élemens , les valeurs de r et U , qui satisfont à l’observation 
du 5.décembre. seront 

r=38°5o' 5* 

Î3 =70 33 34 

Si avec ces cinq élémens on calcule les lieux de la Comète 
aux temps des observations rapportées dans le tableau pré- 
cédent , on trouvera les différences suivantes entre l’observa- 
tion et le calcul. 



16 IVovem. 

î/Novem. 

i8Novem. 

o 3 Noveni . 

3o ^forera . 

5 Déc. 

E 

e 

— »' 5 o* 
— 0 16 

— 0' 3 ' 
+0 35 

+ o' 17* 
— 0 i 3 

+ * 4 ? r 

—0 ag 

•f a' 3 G* 
-o 44 

0 

0 


La plus grande erreur a' 36' tombe sur l’observation du 3o 
novembre ; on auroit pu aisément la diminuer par une solu- 
tion plus recherchée de nos dernières équations ; mais peut-être 
n auroit-on pas réussi aussi bien à atténuer les autres erreurs. 
En général ectte solution paraît très-satisfaisante , puisqu’on 
peut supposer qu une observation est moins précise que les 
autres, et nous pouvons établir avec confiance les élémens 
de 1 orbite de la seconde Comète de i8o5, comme il suit : 
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I.og- dist. périh : 9 • güoajoo 

Passage au périh Décembre 3i.a855i 

Inclinaison i6*5i'io* 

Lieu du nœud asc a5o 53 34 

Lieu du périh. sur l'orbite 109 a3 3g 

Mouvement direct. 

Il est facile d’en conclure que le 8 décembre , à a‘ 5a' , temps 
moyen,] la Comète s’est trouvée dans sa plus grande proxi- 
mité à l’orbite terrestre, dont elle n’étoit éloignée que de 
0.01488 de la distance du soleil, ce qui ne fait pas 6 foi» 
la distance de la lune. 

Calcul de la même orbite par trois observations 
inégalement distantes, 

38. Si la condition de l’égalité de distance entre les ob- 
servations, sert à faciliter le calcul de l’orbite, il faut avouer 
que cette facilité est balancée par l’embarras d’un calcul pré- 
liminaire d’interpolation, et par l’incertitude qu il y a néces- 
sairement sur son résultat, surtout lorsque le lieu calculé est 
à u a-pcu grande des lieux observés. Ainsi, quoique 

la solution du problème entraîne des calculs plus compliqués , 
quand on ne suppose pas égaux les intervalles de temps entre 
les observations ; cependant il peut être plus avantageux à 
d’autres égards, d’appliquer directement lo calcul aux trois 
observations données sans faire usage d’aucune interpolation. 
Pour répandre quelque lumière sur cet objet , et vérifier èn 
même temps les formules trouvées ci-dessus pour la solution 
générale du problème, nous allons donner , avec tout le détail 
nécessaire , l’application de ces formules aux observations des 
a3 novembre, 5o novembre et 5 décembre. 

59. 
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39. Voici donc les données d’après lesquelles nous effec- 
tuerons le calcul : on a diminué les longitudes d’une même quan- 
tité pour rendre nulle celle du 3o novembre. 


9 =7 . 18854 

a°= 9° a ' a 9* 

b‘= 2 f 

35 ' 19 ' 

A=5a*4 7 ' 1 1*5 

6' =4. 78486 

a — 000 

b =19 

35 6 



a'=346 37 5 j 

b'= 3 

19 33 

R=[ 9 .995Ô655](‘) 


Avec les temps 8, ô\ multipliés parle nombre [8.a3558ai ], on 

g/ 1 Q g' g 

forme les quantités r = — - — , r=— ^ — . Avec les longitudes 
a’, a, a et les latitudes b", b, b', on forme les quantités 

- . _ r . , . tang b° tang b . tang b' 

f = tang a" , f — tang a , = g=^, = 

On détermine ensuite les coefficiens F, F', F*, G, G', G*, 
par les formules 

p_. f + f * p' f — t?» _ f — F-t-i'F' 

G G'=i=£,G-= 

a ar 9 wt> 

4o. Enfin les quantités M, N, M', N' , relatives à l’orbite de la 
terre , se calculent par les formules 

IVTt — R cos A. , M'= — ^ — J sin A •— 1 sin'V cos A 

N=Rj»«A, N'=-|-( ■ ^ cos A— -r sin’lr sin A 


(*) Nous désignons par [ g . gg36655 ] le nombre dont le logarithme est 
g.9956655; notation commode et très-propre à prévenir les erreurs désignés 
dans les calculs tels que celui qu’on va développer, et où l'on considère plus 
souvent les logarithmes des quantités, que les quantités elles-mêmes. On prend 
la somme ou la différence de deux expressions de cette sorte , en suivant la 


formule log (a±5) ss log a-f- log^ 1 


7 


C 5 o ) 


oh l’on peut supposer é=[8.335o5i], i — i i’= [9.999959] , 
et l’anomalie vraie Ÿ=long.© — long, apogée ( La longitude 
de l'apogée est maintenant de 99’ 35 ', et elle augmente d’en- 
viron l' chaque aunée ). Ces valeurs donneront le terme 1 sin ’L 
avec une exactitude suffisante; on pourroit même n’emprunter 
des tables que la seule longitude du soleil , et calculer log R 
par la formule 7 og'R=mod. ( écoj'F— t*+[ Scos’^). 

4 i. D'après ocs diverses formules, on trouvera les valeurs 
logarithmiques suivantes 


= [9-°« 376963'/“ 


T 

66 '= 

N = 
M'= 
N'= 


-[ 8 . 3 154387 
[8.00767 1 8] 
[9. 7,75268]! 

[ 9 - 8< J / t 79 ° ] 
-[9. 9 1031 3 ] 

7 8a 9 5a ] 


[9. 20 17 584 ]! 


V « 

f =■ 

F =■ 

F' = 

F 

N —FM =[9. 908040]] 
N'~FM=[g. 7 5853 i]| 


-[<>38 149" 3] 
-[8. 6 io 5 7 8] 
: — [o. 288087] 
= [o. 900416] 


s ’ = 

[9.7304607] 

s ~ 

[9.5471780] 

/ 

8 = 

[8.7760966] 

G = 

[9. 466186] 

G=- 

-[0. 354 a 58 ] 

G*= 

[1. 020553] 

GM = 

[9. 341454] 

GM=- 

-[9. 376399] 


Soit pour abréger /*' — £r'=û.', 
on aura les deux équations 


F>-F'fe *' ( N — FM ) -Ki «r' C NV-FM' ) 

G>— G' 5 '= 3 .'C + —GM') 

qui en substituant les valeurs des coefficiens, deviennent 


H [0.9004 1 6] -f-Ç [0.288087]= «'(9.908040] — ai [7.596839] 

/r [ 1 . o3o55a ]+? [o. 354 a 58 ] = — «.' [9.341 454] — a [7". 2 14707] 
Pour éliminer a> de ces deux équations, je multiplie la pre- 
mière par [9. 333 / t i 4 ], et j’ajoute le produit à la seconde, j’ai 

/ 3 [o. 25583 o] +$[9-621495] J j — a [6! 920253] 

+^[1 . 02055 a] +$[o. 354 a 5 F] \ ( — w[ 7. '21 4707] 
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A^C i-o86a85J+Ç[o. 427983]=— »[ 7 .5y6594] 
ce qui donne 

£= — /*[ o. 65830a ] — u [6. 96841 1 ] 

Mais en substituant les valeurs des coefficiens, on a 

^ =ft' + pr [8 • » 3 g 53 7 ] + a [ 7 . 789666]- 

Donc enfin 

/u'= — /*[ o.6583o2] — [8". 1 3q 33 7 ] — a>[ 7 .85o 7 J2]. 

43. C’est l’équation principale du problème , laquelle est 
représentée généralement par la formule (/-'•), donnée ci-dessus, 
et qui doit être combinée avec les autres équations , au nombre 
de six , savoir : 

r* = R*-H*(2RccsA)-Bf ( ^5) 

1 1 

“ — r* IV ■ : 

m' — M' 

n'=N' + ( F + r'F')At' + F> + a.T , (N — FM) 

p‘~ GM ) 

î = ,»'* + «'*+;/*. 


Or je remarque que ces équations admettront toujours une 
solution - , car si l’on fait /*= o , on aura r = R, &>=o,/u.'=o, 
»'=N', p'=o, et la dernière équation donnant... 

y =M'*+N'*=-^- — 1, ou 7 = A— on vo jt q Ue cette se- 
conde valeur de y est plus petite que la première , qui est Au 
contraire, si l’on fait ji=infini, la valeur de y, déduite de la 


première équation sera 


( 5 a) 


l cns b 



et la râleur déduite de la 


dernière équation, sera de la forme i=A^*} desorte que la 
seconde valeur de £ sera infiniment plus grande que la pre- 
mière. Donc entre les deux valeurs p=o , */ = infini, il jr en 
aura nécessairement une qui rendra identiques les deux valeurs 
de r , données par la première et par la dernière des équations 
précédentes. On parviendra donc toujours à une solution. 

43 . Par la substitution des valeurs numériques , les sept 
équations précédentes deviendront 

r*= 0.971250 -f-ft [o. 07639 d]-f-/**[o.o 5 o 869 ] 

a* = pp — 1.044729 

p-' — — p. [o.G 583 o 2 ] — [8. 1 39337] — m [7.8507 1 2] 

m'— — [9.910213] 

ri— [9.783953] — p [6.818592] — /*[o. 32808 1 ] — t» [8.323469] 
p' = /*' [9.530559 J — p [o. 554 25>8]-f- «[7. 556883 ] 
- r =nï'+ri‘-\-p'‘ 

et on trouvera facilement , après quelques essais , la solution 
suivante : 

, n*'= s — - £ 0 . 7ÜÔ83TÏ ] 

r= [0.0205995] , ri— [9.632016] 

P 1 — — [9.574554] 

au moyen de laquelle les deux valeurs de r s’accordent jusques 
dans la sixième décimale. 

44. Il faut ensuite calculer les coefficiens m, n,p, par les 

formules 4 . 

7n=M-{- ^=[9.841661] 

N 4-//*= [9 . 894790 ] 

P— [8- 540406] 
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mn — m'n = [0.107505] 

mp — mp— [g.356i52] 
np' — n'p =. — [9.490064] 


Puisque la quantité mn’ — m’n est positive, il s’ensuit que le 
mouvement de la Comité est direct ; comme au contraire , 
ainsi qu’il est facile de le voir, la quantité mm' -\-nri -\-pp‘ est 
négative, on en conclura qu'au 5o novembre la Comète n'avoit 
pas encore passé au périhélie. Le lieu du nœud S se calcule 
par la formule 




d’où résulte S=54*56'47*7, ou S = 234°56'47 l '7. L’incertitude 
sera fixée par la valeur de tang I, qui doit être positive. Or on a 


Tang I 


mp' — m' p 
( mn — m n) eos S 


[g.aa8b'47 ] 
coj S ’ 


donc il faut que cos S soit négatif, et qu’on ait S= 334’56'47*7. 
Le nœud ainsi déterminé est le nœud ascendant ( art. XXV ) , 
puisque le mouvement est direct et la latitude boréale. On aura 
ensuite tang I=[ 9.46947^ ] > et parconséquent l’inclinaison 
ï== 16” 55 ^ 5*5 ; 


de là résulte la distance périhélie 

n — = [9-9501644 ]. 

Ensuite l’anomalie se déduit de la formule 


cos*i4=y = [9. 9395649] 
cos i4= [9-9 6 47® a 4l 

4 = 45* 3 1 ' 5 7 ' 7 . 


% 
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Cette anomalie correspond dans la fable des Comités au temps 
T = 36> 5a 1 7a , et par cclui-ci on trouve le temps employé à 
parcourir l’anomalie 4 : 

r=rrT=[ 1.4877979] = 50.74666 

Epoque de l’obs. moy... Novemb.. . 5o.5iog5 

Passage au périhélie... Décemb... 5i. 35761 

La longitude béliocentriquc <p de la Comète se calcule par la 
formule 

Tang p = £ = [o.oSSiago] , 

qui donne (p=48° 39' 45*5, ou $ = 328’ 39' 45*5. La première 
est celle qu’il faut choisir , «pareeque sin <p doit être du même 
sigue que n. Enfin on calculera la distance or du périhélie au 
nœud par la formule 

tang <r= ^ j L r s ) — — [9.07 .4146], 

d’où résulte ( <r devant être dans le même quadraus que p — S) 
<r= 173° i6'38"6, et le lfeu du périhélie sur l’orbite 

•S + c -f- 4 — 93’ 45' 4’* 

R assem b la n t mi<li»«ri résultats, et ajoctiin* rS**W' 56' aux lon- 
gitudes fictives pour avoir les longitudes vraies, on aura les 
élémens suivans : 

Log. dist. périh 9.9501644 

Passage au périh Décembre 31.35761 

Inclinaison 16° 35' 35’ 

Lieu du nœud asc s5o 35 24 

Lieu du pér. sur l’orbite 109 25 40 

En comparant ce résultat à celui des calculs précédées , on 
est étonné de voir combien il approche de la vérité. Une solu- 
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tion si satisfaisante, obtenue de premier abord, sans aucune 
préparation des lieux observés , prouve que nos formules ont 
toute l’exactitude qu'on peut désirer, et que leur application peut 
se faire immédiatement à trois observations inégalement dis- 
tantes, sans donner lieu à des calculs trop compliqués. 

Au reste nous avons eu soin de rassembler dans cet exemple 
toutes les formules nécessaires pour calculer les élémens d'une 
orbite, d’après trois observations données; desorte qu’en prenant 
cet exemple pour guide, on n’aura pas besoin d’avoir recours à 
aucun autre article de cet ouvrage. Lorsque ces Méthodes seront 
connues, j’ose croire qu'on les jugera préférables à tontes les 
autres, et surtout à celles qui, fondées seulement sur des cons- 
tructions approximatives, n’ont pas été soumises à l’examen 
d’une analyse exacte et approfondie. 

Paris , le l 5 janvier 180S. 

*. . * 


ERRA TA du Mémoire précédent. 

„ .. ~ 77 3Kt ~p~ 3 ht ZKl 

Page a, ligne dermere, — -jrr . lues .... 

> 43 s- a5 . le signe h , lisez le signe de //. 

ao, lig. 17, l'équation ( 34 ), Usez l’équation ( 35 ). 

Z fr $ \ 

an , lig. 5 , au lieu du second coefficient H sei — J.; 

a A oa g 

Ibid. lig. 6 et 7 , même correction. 

Page 3 S , lig. ai , cos C, lisez cos. c. 

Page 41 , lig. 1 , 4i" , lù»* 44* • 

Ib. lig. a4, 1781, lisez 1780. 

Page 54 , au bas, { log. Il, lisez J- log. n. 
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